Implementation#
We want to solve the PDE-constrained shape optimization problem
subject to \((\Omega,u_\Omega)\) satisfying
where \(\Omega \subset \mathbb{R}^2\) and \(u_d, f \in C^1(\mathbb{R}^2)\).
from ngsolve import *
from netgen.occ import *
from ngsolve.webgui import Draw
import matplotlib.pyplot as plt
from netgen.geom2d import SplineGeometry
import numpy as np
from ngsolve.nonlinearsolvers import Newton
maxH = 0.05
geo = OCCGeometry(Circle((0, 0), 0.9).Face(), dim=2).GenerateMesh(maxh=maxH)
mesh = Mesh(geo)
vol = 3e-2
a = 0.6
b = 0.4
ud = 1 - ((x-0.01)**2/a**2 + (y-0.01)**2/b**2)
f = (2/a**2 + 2/b**2) + ud
grad_f = CoefficientFunction( (f.Diff(x), f.Diff(y) ) )
grad_ud = CoefficientFunction( (ud.Diff(x), ud.Diff(y) ) )
#Draw(ud, mesh, "target_function")
fes = H1(mesh, order=2, dirichlet=".*")
u, v = fes.TnT()
gfu = GridFunction(fes)
scene_gfu = Draw (gfu, mesh, "state")
a = BilinearForm(fes, symmetric=True)
a += grad(u)*grad(v)*dx
a += u*v*dx
fstate = LinearForm(fes)
fstate += f*v*dx
def SolveStateEquation():
rhs = gfu.vec.CreateVector()
rhs.data = fstate.vec - a.mat * gfu.vec
update = gfu.vec.CreateVector()
update.data = a.mat.Inverse(fes.FreeDofs()) * rhs
gfu.vec.data += update
a.Assemble()
fstate.Assemble()
SolveStateEquation()
scene_gfu.Redraw()
Adjoint equation#
The adjoint p satisfies
def Cost(u):
cost = InnerProduct(grad(u) - grad_ud, grad(u) - grad_ud) * dx
cost += InnerProduct(CF(vol), CF(1)) * dx
return cost
p, w = fes.TnT()
gfp = GridFunction(fes)
scene_gfp = Draw (gfp, mesh, "adjoint")
fadjoint = LinearForm(fes)
#fadjoint += -1*Cost(gfu).Diff(gfu,w)
fadjoint += -2* InnerProduct(grad(gfu) - grad_ud, grad(w)) * dx
def SolveAdjointEquation():
rhs = gfp.vec.CreateVector()
rhs.data = fadjoint.vec - a.mat.T * gfp.vec
update = gfp.vec.CreateVector()
update.data = a.mat.Inverse(fes.FreeDofs()).T * rhs
gfp.vec.data += update
fadjoint.Assemble()
SolveAdjointEquation()
scene_gfp.Redraw()
The shape derivative of the shape function
at \(\Omega\) in direction \(X\) is given by
We now assemble this shape derivative in NGSolve as follows:
def hesse(u):
return CoefficientFunction(
(u.Diff(x).Diff(x), u.Diff(x).Diff(y),
u.Diff(y).Diff(x), u.Diff(y).Diff(y)),
dims=(2,2)
)
#VEC = H1(mesh, order=3, dirichlet=".*", dim=2)
VEC = H1(mesh, order=1, dim = 2)
X, PHI = VEC.TnT()
dJOmega = LinearForm(VEC)
dJOmega += InnerProduct((div(PHI)*(grad(gfu)- grad_ud)), grad(gfu)- grad_ud) * dx
dJOmega += -2* InnerProduct(hesse(ud)*PHI, grad(gfu)- grad_ud) * dx
dJOmega += div(PHI)* gfu *gfp *dx
dJOmega += InnerProduct(div(PHI)*grad(gfu), grad(gfp)) * dx
dJOmega += -InnerProduct(PHI.Deriv()*grad(gfu), grad(gfp)) * dx
dJOmega += -InnerProduct(grad(gfu), PHI.Deriv()*grad(gfp)) * dx
dJOmega += -div(PHI)* f * gfp * dx
dJOmega += -grad_f * PHI*gfp * dx
# added volumen constraint
dJOmega += -vol * div(PHI)* dx
Next, we want to find a vector field \(X\) which yields a decrease of the objective function, i.e. we want to find a vector field \(X\) such that
\[ d\mathcal{J}(\Omega; X) < 0. \]This can be achieved by solving an auxiliary boundary value problem of the type
\[ \text{Find } X \in H : \qquad b(X, \Phi) = d\mathcal{J}(\Omega; \Phi) \quad \text{for all } \Phi \in H, \]for a suitable Hilbert space \(H\) (e.g. \(H = H^1(\Omega)^2\)). Here, \(b(\cdot,\cdot) : H \times H \to \mathbb{R}\) is a positive definite bilinear form which we are free to choose. Then \(-X\) is a descent direction since
\[ d\mathcal{J}(\Omega; -X) = - d\mathcal{J}(\Omega; X) = - b(X, X) < 0. \]First we start with the bilinear form
\[ b(X, Y) = \int_\Omega \nabla X : \nabla Y + X \cdot Y \, dx . \]or
\[ \tilde{b}(X,Y) = b(X, Y) + \gamma \int BX \cdot BY dx \]with $\( B := \begin{pmatrix} -\partial_x & \partial_y \\ \partial_y & \partial_x \end{pmatrix}.\)$
we can also choose X as $\( \min _{X \in H^1(\Omega)} \frac{1}{p} \int _\Omega |\nabla X|^p - DJ(\Omega)(X) \)$
since this is a continous functional on \(H^1(\Omega)'\), it is equivalent to finding \(X \in H^1(\Omega)\) sucht that,
\[ \int _\Omega |\nabla X | ^{p-2} \nabla X \cdot \nabla \Phi = DJ(\Omega)(\Phi) \qquad \forall \Phi \in H^1(\Omega) \]
def pose_deformation_equation(gamma=None, p = None):
b = BilinearForm(VEC)
# standard case
if not p:
b += InnerProduct(grad(X),grad(PHI))*dx + InnerProduct(X,PHI)*dx
if gamma:
# standard case with gamma regularization
print('Using gamma =', gamma)
b = BilinearForm(VEC)
b += (1.0/gamma)*InnerProduct(grad(X),grad(PHI))*dx + (1.0/gamma)*InnerProduct(X,PHI)*dx
b += gamma * (PHI.Deriv()[0,0] - PHI.Deriv()[1,1])*(X.Deriv()[0,0] - X.Deriv()[1,1]) *dx
b += gamma * (PHI.Deriv()[1,0] + PHI.Deriv()[0,1])*(X.Deriv()[1,0] + X.Deriv()[0,1]) *dx
#print("expression tree\n")
#print(b)
return b
elif p:
# solving minimization problem
print('Using p =', p)
eps = 1E-5
grad_norm = sqrt(InnerProduct(grad(X), grad(X)) + eps**2)
c = BilinearForm(VEC)
c += (grad_norm**(p-2)) * InnerProduct(grad(X), grad(PHI)) * dx
c += InnerProduct(X,PHI)*dx
# -DJ(Omega)(PHI)
c += -InnerProduct((div(PHI)*(grad(gfu)- grad_ud)), grad(gfu)- grad_ud) * dx
c += 2* InnerProduct(grad(gfu)- grad_ud, hesse(ud)*PHI) * dx
c += -div(PHI)* gfu *gfp *dx
c += -InnerProduct(div(PHI)*grad(gfu), grad(gfp)) * dx
c += InnerProduct(PHI.Deriv()*grad(gfu), grad(gfp)) * dx
c += InnerProduct(grad(gfu), PHI.Deriv()*grad(gfp)) * dx
c += div(PHI)* f * gfp * dx
c += grad_f*PHI*gfp * dx
c += vol * div(PHI)*dx
# c += -div(PHI)* f * gfp * dx
# c += -grad_f*PHI*gfp * dx
return c
#print("expression tree\n")
#print(b)
return b
def SolveDeformationEquation(linear=True):
if linear:
#print("solving linear deformation equation")
rhs = gfX.vec.CreateVector()
rhs.data = dJOmega.vec - b.mat * gfX.vec
update = gfX.vec.CreateVector()
update.data = b.mat.Inverse(VEC.FreeDofs()) * rhs
gfX.vec.data += update
#scene_gfset.Redraw()
elif not linear:
#print("solving nonlinear deformation equation with Newton")
Newton(c, gfX, maxerr = 1E-8, freedofs=VEC.FreeDofs(), maxit = 1000, printing=True)
#print("Deformation norm:", Norm(gfX.vec)
#scene_gfset.Redraw()
checking the mesh quality#
We are estimating the meshquality by taking the elementwise (approximated) diameter
def check_mesh_quality():
h = specialcf.mesh_size
gf_h = GridFunction(L2(mesh, order=0))
gf_h.Set(h)
qual = gf_h.vec.FV().NumPy()
qual /= maxH
print("=================Mesh Quality Check=================")
print("max quality:", np.max(qual), " min quality:", np.min(qual), " mean quality:", np.mean(qual))
print("===================================================")
#return Integrate(h, mesh, element_wise=True)
def optimize_shape(linear = True, LineSearch = True):
gfset.Set((0,0))
mesh.SetDeformation(gfset)
scene.Redraw()
a.Assemble()
fstate.Assemble()
SolveStateEquation()
iter_max = 600
Jold = Integrate(Cost(gfu), mesh)
converged = False
# input("Press enter to start optimization")
print('Cost at initial design', Integrate (Cost(gfu), mesh))
for k in range(iter_max):
mesh.SetDeformation(gfset)
scene.Redraw()
print('cost at iteration', k, ': ', Jold)
a.Assemble()
fstate.Assemble()
SolveStateEquation()
fadjoint.Assemble()
SolveAdjointEquation()
if linear:
b.Assemble()
#plt.spy(b.mat.ToDense())
dJOmega.Assemble()
SolveDeformationEquation(linear=linear)
#scale = 0.01 / Norm(gfX.vec)
scale = 1E-4
gfsetOld.vec[:] = gfset.vec[:]
gfset.vec.data -= scale * gfX.vec
Jnew = Integrate(Cost(gfu), mesh)
if not LineSearch:
if Jnew >= Jold:
print("No descent found for current step size")
converged = True
#print("gradient norm", np.linalg.norm(dJOmega.vec.FV().NumPy()))
if LineSearch:
while Jnew > Jold and scale > 1e-12:
#input('a')
scale = min(1E-3, scale/2)
print("scale = ", scale)
#print("scale without normalization factor = ", scale* Norm(gfX.vec))
if scale <= 1e-12:
converged = True
break
gfset.vec.data = gfsetOld.vec - scale * gfX.vec # scale times -X is descent direction
mesh.SetDeformation(gfset) # applying deformation and solve new state
a.Assemble() # note that during line search we only have to solve state to check if update yields lower costs
fstate.Assemble()
SolveStateEquation()
Jnew = Integrate(Cost(gfu), mesh)
print("New cost during line search: ", Jnew)
#print("scaling parameter:", scale)
#print("descet direction value", -1*np.dot(dJOmega.vec.FV().NumPy(), gfX.vec.FV().NumPy()))
if converged==True:
print("No more descent can be found")
break
Jold = Jnew
Redraw(blocking=True)
#print("Norm of gset", Integrate(gfset**2, mesh))
linear#
a1 = check_mesh_quality()
gfX = GridFunction(VEC)
b = pose_deformation_equation()
gfset = GridFunction(VEC)
gfset.Set((0,0))
gfsetOld = GridFunction(VEC)
scene = Draw(gfset,mesh,"gfset")
SetVisualization (deformation=True)
optimize_shape(linear = True)
mesh.SetDeformation(gfset)
a2 = check_mesh_quality()
# gfX2 = GridFunction(VEC)
# gfX2.vec.data =
# print("=====")
# print(np.dot(dJOmega.vec.FV().NumPy(), gfX.vec.FV().NumPy()))
# print(Integrate(Cost(gfu), mesh))
=================Mesh Quality Check=================
max quality: 1.452321790073649 min quality: 0.5768477774464785 mean quality: 0.8831324760843785
===================================================
Cost at initial design 13.96238810275433
cost at iteration 0 : 13.96238810275433
cost at iteration 1 : 13.735550739621843
cost at iteration 2 : 13.526642840331323
cost at iteration 3 : 13.322392417508466
cost at iteration 4 : 13.12265849375249
cost at iteration 5 : 12.927305689331154
cost at iteration 6 : 12.736203943744316
cost at iteration 7 : 12.549228253853594
cost at iteration 8 : 12.36625842744546
cost at iteration 9 : 12.187178851181772
cost at iteration 10 : 12.011878271971462
cost at iteration 11 : 11.840249590867716
cost at iteration 12 : 11.67218966866377
cost at iteration 13 : 11.507599142420847
cost at iteration 14 : 11.34638225221715
cost at iteration 15 : 11.188446677459389
cost at iteration 16 : 11.03370338214665
cost at iteration 17 : 10.88206646851788
cost at iteration 18 : 10.733453038557334
cost at iteration 19 : 10.587783062867315
cost at iteration 20 : 10.44497925645389
cost at iteration 21 : 10.304966961000465
cost at iteration 22 : 10.167674033236612
cost at iteration 23 : 10.03303073903283
cost at iteration 24 : 9.900969652880633
cost at iteration 25 : 9.771425562437086
cost at iteration 26 : 9.64433537783695
cost at iteration 27 : 9.519638045494702
cost at iteration 28 : 9.397274466134974
cost at iteration 29 : 9.277187416810861
cost at iteration 30 : 9.159321476681768
cost at iteration 31 : 9.043622956339174
cost at iteration 32 : 8.930039830481821
cost at iteration 33 : 8.818521673753528
cost at iteration 34 : 8.709019599570668
cost at iteration 35 : 8.601486201774575
cost at iteration 36 : 8.495875498956943
cost at iteration 37 : 8.392142881313857
cost at iteration 38 : 8.290245059893905
cost at iteration 39 : 8.190140018113214
cost at iteration 40 : 8.091786965419237
cost at iteration 41 : 7.995146292990532
cost at iteration 42 : 7.900179531368054
cost at iteration 43 : 7.806849309918279
cost at iteration 44 : 7.715119318036037
cost at iteration 45 : 7.624954267997505
cost at iteration 46 : 7.53631985938323
cost at iteration 47 : 7.449182744990578
cost at iteration 48 : 7.363510498164074
cost at iteration 49 : 7.279271581473583
cost at iteration 50 : 7.1964353166745845
cost at iteration 51 : 7.114971855889361
cost at iteration 52 : 7.034852153950404
cost at iteration 53 : 6.956047941851222
cost at iteration 54 : 6.878531701251429
cost at iteration 55 : 6.802276639988421
cost at iteration 56 : 6.727256668547503
cost at iteration 57 : 6.653446377446716
cost at iteration 58 : 6.580821015495286
cost at iteration 59 : 6.509356468885111
cost at iteration 60 : 6.4390292410776775
cost at iteration 61 : 6.369816433451991
cost at iteration 62 : 6.30169572667845
cost at iteration 63 : 6.2346453627871155
cost at iteration 64 : 6.168644127900511
cost at iteration 65 : 6.1036713356005805
cost at iteration 66 : 6.039706810904678
cost at iteration 67 : 5.97673087482172
cost at iteration 68 : 5.91472432946657
cost at iteration 69 : 5.853668443706879
cost at iteration 70 : 5.793544939320945
cost at iteration 71 : 5.734335977645499
cost at iteration 72 : 5.67602414669258
cost at iteration 73 : 5.618592448716386
cost at iteration 74 : 5.562024288211479
cost at iteration 75 : 5.506303460325568
cost at iteration 76 : 5.45141413966889
cost at iteration 77 : 5.3973408695056415
cost at iteration 78 : 5.344068551311047
cost at iteration 79 : 5.291582434680406
cost at iteration 80 : 5.239868107575666
cost at iteration 81 : 5.188911486897126
cost at iteration 82 : 5.138698809366406
cost at iteration 83 : 5.089216622710051
cost at iteration 84 : 5.040451777131078
cost at iteration 85 : 4.992391417058067
cost at iteration 86 : 4.9450229731612945
cost at iteration 87 : 4.898334154625466
cost at iteration 88 : 4.852312941669703
cost at iteration 89 : 4.806947578305664
cost at iteration 90 : 4.762226565324642
cost at iteration 91 : 4.718138653505423
cost at iteration 92 : 4.6746728370347395
cost at iteration 93 : 4.631818347132968
cost at iteration 94 : 4.589564645876579
cost at iteration 95 : 4.547901420211695
cost at iteration 96 : 4.506818576150781
cost at iteration 97 : 4.466306233146395
cost at iteration 98 : 4.426354718635686
cost at iteration 99 : 4.386954562749661
cost at iteration 100 : 4.348096493181397
cost at iteration 101 : 4.309771430207505
cost at iteration 102 : 4.271970481857837
cost at iteration 103 : 4.234684939228223
cost at iteration 104 : 4.197906271931019
cost at iteration 105 : 4.161626123679323
cost at iteration 106 : 4.125836307999768
cost at iteration 107 : 4.090528804069933
cost at iteration 108 : 4.055695752676047
cost at iteration 109 : 4.02132945228679
cost at iteration 110 : 3.98742235523977
cost at iteration 111 : 3.953967064036446
cost at iteration 112 : 3.9209563277422292
cost at iteration 113 : 3.8883830384882834
cost at iteration 114 : 3.856240228071364
cost at iteration 115 : 3.8245210646493395
cost at iteration 116 : 3.793218849528267
cost at iteration 117 : 3.7623270140387124
cost at iteration 118 : 3.7318391164983775
cost at iteration 119 : 3.7017488392583173
cost at iteration 120 : 3.672049985829624
cost at iteration 121 : 3.64273647808877
cost at iteration 122 : 3.6138023535586794
cost at iteration 123 : 3.5852417627631556
cost at iteration 124 : 3.557048966652588
cost at iteration 125 : 3.52921833409841
cost at iteration 126 : 3.501744339454925
cost at iteration 127 : 3.474621560184901
cost at iteration 128 : 3.44784467454924
cost at iteration 129 : 3.421408459356327
cost at iteration 130 : 3.3953077877712023
cost at iteration 131 : 3.369537627181565
cost at iteration 132 : 3.344093037119673
cost at iteration 133 : 3.3189691672374315
cost at iteration 134 : 3.294161255334572
cost at iteration 135 : 3.2696646254369646
cost at iteration 136 : 3.245474685924152
cost at iteration 137 : 3.221586927704939
cost at iteration 138 : 3.1979969224390046
cost at iteration 139 : 3.174700320803726
cost at iteration 140 : 3.151692850804505
cost at iteration 141 : 3.1289703161277425
cost at iteration 142 : 3.106528594534672
cost at iteration 143 : 3.0843636362954583
cost at iteration 144 : 3.062471462662174
cost at iteration 145 : 3.040848164379228
cost at iteration 146 : 3.019489900230858
cost at iteration 147 : 2.998392895623941
cost at iteration 148 : 2.9775534412056017
cost at iteration 149 : 2.9569678915144437
cost at iteration 150 : 2.9366326636644047
cost at iteration 151 : 2.916544236060737
cost at iteration 152 : 2.8966991471463954
cost at iteration 153 : 2.877093994178978
cost at iteration 154 : 2.8577254320364656
cost at iteration 155 : 2.838590172051993
cost at iteration 156 : 2.819684980875429
cost at iteration 157 : 2.8010066793627577
cost at iteration 158 : 2.7825521414908816
cost at iteration 159 : 2.7643182932984156
cost at iteration 160 : 2.746302111851069
cost at iteration 161 : 2.7285006242307546
cost at iteration 162 : 2.710910906549088
cost at iteration 163 : 2.693530082982694
cost at iteration 164 : 2.6763553248313885
cost at iteration 165 : 2.6593838495978726
cost at iteration 166 : 2.642612920088303
cost at iteration 167 : 2.626039843533671
cost at iteration 168 : 2.6096619707309516
cost at iteration 169 : 2.593476695203691
cost at iteration 170 : 2.577481452381752
cost at iteration 171 : 2.561673718799414
cost at iteration 172 : 2.546051011311209
cost at iteration 173 : 2.5306108863257855
cost at iteration 174 : 2.515350939056333
cost at iteration 175 : 2.5002688027880033
cost at iteration 176 : 2.4853621481613435
cost at iteration 177 : 2.470628682471653
cost at iteration 178 : 2.4560661489834916
cost at iteration 179 : 2.4416723262605333
cost at iteration 180 : 2.427445027509923
cost at iteration 181 : 2.413382099940702
cost at iteration 182 : 2.3994814241364444
cost at iteration 183 : 2.3857409134412375
cost at iteration 184 : 2.3721585133589196
cost at iteration 185 : 2.358732200965482
cost at iteration 186 : 2.3454599843335524
cost at iteration 187 : 2.3323399019696147
cost at iteration 188 : 2.3193700222628615
cost at iteration 189 : 2.3065484429460783
cost at iteration 190 : 2.293873290567446
cost at iteration 191 : 2.2813427199739977
cost at iteration 192 : 2.2689549138054415
cost at iteration 193 : 2.256708081998858
cost at iteration 194 : 2.244600461303798
cost at iteration 195 : 2.2326303148069733
cost at iteration 196 : 2.2207959314676313
cost at iteration 197 : 2.2090956256614747
cost at iteration 198 : 2.1975277367350303
cost at iteration 199 : 2.1860906285684303
cost at iteration 200 : 2.1747826891474222
cost at iteration 201 : 2.1636023301441427
cost at iteration 202 : 2.1525479865062165
cost at iteration 203 : 2.141618116054566
cost at iteration 204 : 2.130811199088807
cost at iteration 205 : 2.1201257380012137
cost at iteration 206 : 2.109560256897829
cost at iteration 207 : 2.099113301227706
cost at iteration 208 : 2.0887834374192007
cost at iteration 209 : 2.078569252523587
cost at iteration 210 : 2.0684693538658676
cost at iteration 211 : 2.0584823687022507
cost at iteration 212 : 2.0486069438847796
cost at iteration 213 : 2.038841745532014
cost at iteration 214 : 2.029185458706767
cost at iteration 215 : 2.0196367870996013
cost at iteration 216 : 2.010194452718936
cost at iteration 217 : 2.0008571955867405
cost at iteration 218 : 1.9916237734406657
cost at iteration 219 : 1.9824929614414906
cost at iteration 220 : 1.9734635518862256
cost at iteration 221 : 1.9645343539271067
cost at iteration 222 : 1.9557041932955332
cost at iteration 223 : 1.9469719120314786
cost at iteration 224 : 1.9383363682180776
cost at iteration 225 : 1.929796435721272
cost at iteration 226 : 1.9213510039341817
cost at iteration 227 : 1.9129989775265157
cost at iteration 228 : 1.9047392761987394
cost at iteration 229 : 1.8965708344405479
cost at iteration 230 : 1.8884926012941636
cost at iteration 231 : 1.8805035401220322
cost at iteration 232 : 1.8726026283785673
cost at iteration 233 : 1.8647888573865248
cost at iteration 234 : 1.8570612321171416
cost at iteration 235 : 1.8494187709746175
cost at iteration 236 : 1.841860505584373
cost at iteration 237 : 1.8343854805853856
cost at iteration 238 : 1.8269927534260753
cost at iteration 239 : 1.8196813941641619
cost at iteration 240 : 1.8124504852699752
cost at iteration 241 : 1.8052991214334695
cost at iteration 242 : 1.7982264093745732
cost at iteration 243 : 1.7912314676572485
cost at iteration 244 : 1.7843134265065905
cost at iteration 245 : 1.7774714276292956
cost at iteration 246 : 1.7707046240374826
cost at iteration 247 : 1.7640121798754804
cost at iteration 248 : 1.757393270249858
cost at iteration 249 : 1.7508470810623302
cost at iteration 250 : 1.7443728088456039
cost at iteration 251 : 1.737969660602488
cost at iteration 252 : 1.731636853647192
cost at iteration 253 : 1.725373615450051
cost at iteration 254 : 1.7191791834845536
cost at iteration 255 : 1.7130528050772313
cost at iteration 256 : 1.7069937372601294
cost at iteration 257 : 1.70100124662571
cost at iteration 258 : 1.6950746091846793
cost at iteration 259 : 1.6892131102256227
cost at iteration 260 : 1.683416044177652
cost at iteration 261 : 1.6776827144752418
cost at iteration 262 : 1.6720124334249806
cost at iteration 263 : 1.6664045220754165
cost at iteration 264 : 1.6608583100883452
cost at iteration 265 : 1.6553731356126518
cost at iteration 266 : 1.6499483451603318
cost at iteration 267 : 1.6445832934844677
cost at iteration 268 : 1.639277343459288
cost at iteration 269 : 1.6340298659624126
cost at iteration 270 : 1.6288402397588018
cost at iteration 271 : 1.6237078513869072
cost at iteration 272 : 1.6186320950464848
cost at iteration 273 : 1.6136123724886295
cost at iteration 274 : 1.6086480929071705
cost at iteration 275 : 1.6037386728322989
cost at iteration 276 : 1.5988835360257612
cost at iteration 277 : 1.5940821133777066
cost at iteration 278 : 1.5893338428053423
cost at iteration 279 : 1.5846381691534075
cost at iteration 280 : 1.5799945440958447
cost at iteration 281 : 1.5754024260395718
cost at iteration 282 : 1.5708612800294617
cost at iteration 283 : 1.5663705776550487
cost at iteration 284 : 1.56192979695875
cost at iteration 285 : 1.557538422345481
cost at iteration 286 : 1.5531959444936536
cost at iteration 287 : 1.548901860268023
cost at iteration 288 : 1.5446556726333833
cost at iteration 289 : 1.5404568905700111
cost at iteration 290 : 1.5363050289903868
cost at iteration 291 : 1.5321996086571132
cost at iteration 292 : 1.5281401561023675
cost at iteration 293 : 1.5241262035484744
cost at iteration 294 : 1.5201572888296888
cost at iteration 295 : 1.5162329553153822
cost at iteration 296 : 1.5123527518343562
cost at iteration 297 : 1.508516232600388
cost at iteration 298 : 1.5047229571387146
cost at iteration 299 : 1.5009724902141657
cost at iteration 300 : 1.4972644017598626
cost at iteration 301 : 1.493598266807315
cost at iteration 302 : 1.4899736654177438
cost at iteration 303 : 1.4863901826139896
cost at iteration 304 : 1.4828474083142
cost at iteration 305 : 1.4793449372656573
cost at iteration 306 : 1.475882368980472
cost at iteration 307 : 1.472459307671675
cost at iteration 308 : 1.4690753621905175
cost at iteration 309 : 1.4657301459648182
cost at iteration 310 : 1.4624232769379895
cost at iteration 311 : 1.459154377509365
cost at iteration 312 : 1.4559230744750364
cost at iteration 313 : 1.4527289989697263
cost at iteration 314 : 1.4495717864098823
cost at iteration 315 : 1.4464510764369465
cost at iteration 316 : 1.4433665128621622
cost at iteration 317 : 1.4403177436117856
cost at iteration 318 : 1.4373044206732672
cost at iteration 319 : 1.4343262000422319
cost at iteration 320 : 1.431382741670177
cost at iteration 321 : 1.428473709413187
cost at iteration 322 : 1.4255987709810212
cost at iteration 323 : 1.4227575978872453
cost at iteration 324 : 1.4199498654001614
cost at iteration 325 : 1.417175252494063
cost at iteration 326 : 1.4144334418017905
cost at iteration 327 : 1.4117241195673917
cost at iteration 328 : 1.409046975599941
cost at iteration 329 : 1.4064017032277107
cost at iteration 330 : 1.4037879992532847
cost at iteration 331 : 1.4012055639090877
cost at iteration 332 : 1.3986541008136941
cost at iteration 333 : 1.396133316928895
cost at iteration 334 : 1.3936429225169267
cost at iteration 335 : 1.3911826310988555
cost at iteration 336 : 1.3887521594132843
cost at iteration 337 : 1.3863512273758267
cost at iteration 338 : 1.3839795580386887
cost at iteration 339 : 1.3816368775516255
cost at iteration 340 : 1.3793229151226805
cost at iteration 341 : 1.3770374029799521
cost at iteration 342 : 1.3747800763339746
cost at iteration 343 : 1.372550673339937
cost at iteration 344 : 1.3703489350613844
cost at iteration 345 : 1.3681746054339314
cost at iteration 346 : 1.3660274312292202
cost at iteration 347 : 1.3639071620201206
cost at iteration 348 : 1.3618135501456547
cost at iteration 349 : 1.3597463506771514
cost at iteration 350 : 1.3577053213840993
cost at iteration 351 : 1.3556902227012104
cost at iteration 352 : 1.3537008176953111
cost at iteration 353 : 1.3517368720332115
cost at iteration 354 : 1.3497981539496846
cost at iteration 355 : 1.3478844342161382
cost at iteration 356 : 1.3459954861093983
cost at iteration 357 : 1.344131085381208
cost at iteration 358 : 1.342291010228171
cost at iteration 359 : 1.3404750412617914
cost at iteration 360 : 1.3386829614792632
cost at iteration 361 : 1.3369145562345277
cost at iteration 362 : 1.3351696132096162
cost at iteration 363 : 1.333447922386621
cost at iteration 364 : 1.3317492760198009
cost at iteration 365 : 1.3300734686084696
cost at iteration 366 : 1.3284202968695196
cost at iteration 367 : 1.3267895597111736
cost at iteration 368 : 1.325181058206404
cost at iteration 369 : 1.3235945955673651
cost at iteration 370 : 1.322029977119296
cost at iteration 371 : 1.3204870102757758
cost at iteration 372 : 1.3189655045134994
cost at iteration 373 : 1.3174652713478447
cost at iteration 374 : 1.3159861243086943
cost at iteration 375 : 1.314527878916308
cost at iteration 376 : 1.3130903526578683
cost at iteration 377 : 1.3116733649644474
cost at iteration 378 : 1.310276737187427
cost at iteration 379 : 1.3089002925764317
cost at iteration 380 : 1.3075438562566144
cost at iteration 381 : 1.306207255206912
cost at iteration 382 : 1.3048903182377662
cost at iteration 383 : 1.3035928759701794
cost at iteration 384 : 1.3023147608141397
cost at iteration 385 : 1.3010558069479445
cost at iteration 386 : 1.2998158502972734
cost at iteration 387 : 1.2985947285150612
cost at iteration 388 : 1.29739228096135
cost at iteration 389 : 1.2962083486832132
cost at iteration 390 : 1.2950427743954231
cost at iteration 391 : 1.293895402461041
cost at iteration 392 : 1.2927660788722422
cost at iteration 393 : 1.2916546512315608
cost at iteration 394 : 1.2905609687333333
cost at iteration 395 : 1.2894848821451101
cost at iteration 396 : 1.2884262437899847
cost at iteration 397 : 1.2873849075281925
cost at iteration 398 : 1.2863607287398187
cost at iteration 399 : 1.2853535643072802
cost at iteration 400 : 1.2843632725979808
cost at iteration 401 : 1.2833897134475774
cost at iteration 402 : 1.2824327481429731
cost at iteration 403 : 1.2814922394058883
cost at iteration 404 : 1.2805680513765503
cost at iteration 405 : 1.2796600495974786
cost at iteration 406 : 1.2787681009975265
cost at iteration 407 : 1.2778920738762698
cost at iteration 408 : 1.2770318378884116
cost at iteration 409 : 1.2761872640284655
cost at iteration 410 : 1.2753582246155069
cost at iteration 411 : 1.2745445932784343
cost at iteration 412 : 1.273746244941086
cost at iteration 413 : 1.2729630558076461
cost at iteration 414 : 1.272194903348229
cost at iteration 415 : 1.271441666284778
cost at iteration 416 : 1.2707032245769816
cost at iteration 417 : 1.2699794594083245
cost at iteration 418 : 1.2692702531725586
cost at iteration 419 : 1.268575489460036
cost at iteration 420 : 1.2678950530445068
cost at iteration 421 : 1.2672288298697183
cost at iteration 422 : 1.266576707036594
cost at iteration 423 : 1.2659385727904116
cost at iteration 424 : 1.2653143165078875
cost at iteration 425 : 1.264703828684697
cost at iteration 426 : 1.2641070009232351
cost at iteration 427 : 1.2635237259201735
cost at iteration 428 : 1.2629538974546575
cost at iteration 429 : 1.2623974103758737
cost at iteration 430 : 1.2618541605918403
cost at iteration 431 : 1.2613240450573697
cost at iteration 432 : 1.2608069617625788
cost at iteration 433 : 1.260302809721812
cost at iteration 434 : 1.259811488962035
cost at iteration 435 : 1.259332900511973
cost at iteration 436 : 1.258866946391148
cost at iteration 437 : 1.2584135295989214
cost at iteration 438 : 1.2579725541039437
cost at iteration 439 : 1.2575439248334177
cost at iteration 440 : 1.2571275476629469
cost at iteration 441 : 1.2567233294058464
cost at iteration 442 : 1.256331177803315
cost at iteration 443 : 1.255951001514014
cost at iteration 444 : 1.2555827101045347
cost at iteration 445 : 1.2552262140391546
cost at iteration 446 : 1.2548814246703741
cost at iteration 447 : 1.254548254229384
cost at iteration 448 : 1.254226615816268
cost at iteration 449 : 1.2539164233909819
cost at iteration 450 : 1.2536175917640013
cost at iteration 451 : 1.2533300365868727
cost at iteration 452 : 1.2530536743438068
cost at iteration 453 : 1.25278842234207
cost at iteration 454 : 1.2525341987036143
cost at iteration 455 : 1.252290922356204
cost at iteration 456 : 1.2520585130247126
cost at iteration 457 : 1.251836891222852
cost at iteration 458 : 1.2516259782445036
cost at iteration 459 : 1.2514256961555208
cost at iteration 460 : 1.251235967785526
cost at iteration 461 : 1.2510567167196875
cost at iteration 462 : 1.2508878672909112
cost at iteration 463 : 1.2507293445715988
cost at iteration 464 : 1.2505810743660208
cost at iteration 465 : 1.2504429832025268
cost at iteration 466 : 1.2503149983259212
cost at iteration 467 : 1.2501970476897117
cost at iteration 468 : 1.250089059948902
cost at iteration 469 : 1.2499909644523188
cost at iteration 470 : 1.249902691235647
cost at iteration 471 : 1.2498241710136928
cost at iteration 472 : 1.2497553351738397
cost at iteration 473 : 1.249696115768505
cost at iteration 474 : 1.249646445508458
cost at iteration 475 : 1.2496062577556994
cost at iteration 476 : 1.2495754865168993
cost at iteration 477 : 1.2495540664364273
cost at iteration 478 : 1.2495419327897916
cost at iteration 479 : 1.2495390214770865
scale = 5e-05
New cost during line search: 1.242616032021647
cost at iteration 480 : 1.242616032021647
scale = 5e-05
New cost during line search: 1.2426156656099823
cost at iteration 481 : 1.2426156656099823
scale = 5e-05
New cost during line search: 1.242617568717261
scale = 2.5e-05
New cost during line search: 1.2426165396036053
scale = 1.25e-05
New cost during line search: 1.2426160832181463
scale = 6.25e-06
New cost during line search: 1.2426158695670388
scale = 3.125e-06
New cost during line search: 1.24261576637678
scale = 1.5625e-06
New cost during line search: 1.2426157156904207
scale = 7.8125e-07
New cost during line search: 1.2426156905745076
scale = 3.90625e-07
New cost during line search: 1.242615678073288
scale = 1.953125e-07
New cost during line search: 1.2426156718368953
scale = 9.765625e-08
New cost during line search: 1.2426156687222547
scale = 4.8828125e-08
New cost during line search: 1.2426156671658437
scale = 2.44140625e-08
New cost during line search: 1.242615666387826
scale = 1.220703125e-08
New cost during line search: 1.2426156659988907
scale = 6.103515625e-09
New cost during line search: 1.2426156658044298
scale = 3.0517578125e-09
New cost during line search: 1.2426156657071838
scale = 1.52587890625e-09
New cost during line search: 1.2426156656586058
scale = 7.62939453125e-10
New cost during line search: 1.242615665634289
scale = 3.814697265625e-10
New cost during line search: 1.242615665622116
scale = 1.9073486328125e-10
New cost during line search: 1.2426156656160388
scale = 9.5367431640625e-11
New cost during line search: 1.2426156656130096
scale = 4.76837158203125e-11
New cost during line search: 1.2426156656114715
scale = 2.384185791015625e-11
New cost during line search: 1.2426156656107425
scale = 1.1920928955078126e-11
New cost during line search: 1.2426156656103344
scale = 5.960464477539063e-12
New cost during line search: 1.2426156656101475
scale = 2.9802322387695314e-12
New cost during line search: 1.2426156656100888
scale = 1.4901161193847657e-12
New cost during line search: 1.2426156656099985
scale = 7.450580596923828e-13
No more descent can be found
=================Mesh Quality Check=================
max quality: 1.4420554859725245 min quality: 0.5166402075400832 mean quality: 0.8593656411107816
===================================================
check_mesh_quality()
gfX = GridFunction(VEC)
b = pose_deformation_equation(gamma=1)
gfset = GridFunction(VEC)
gfset.Set((0,0))
gfsetOld = GridFunction(VEC)
scene = Draw(gfset,mesh,"gfset")
SetVisualization (deformation=True)
optimize_shape(linear = True)
mesh.SetDeformation(gfset)
check_mesh_quality()
=================Mesh Quality Check=================
max quality: 1.4420554859725245 min quality: 0.5166402075400832 mean quality: 0.8593656411107816
===================================================
Using gamma = 1
Cost at initial design 13.962388102754325
cost at iteration 0 : 13.962388102754325
cost at iteration 1 : 13.825772613698412
cost at iteration 2 : 13.708401589313848
cost at iteration 3 : 13.592734781168076
cost at iteration 4 : 13.478731015689526
cost at iteration 5 : 13.366350749245594
cost at iteration 6 : 13.255555977256009
cost at iteration 7 : 13.146310149336143
cost at iteration 8 : 13.038578090036056
cost at iteration 9 : 12.932325924773831
cost at iteration 10 : 12.827521010591616
cost at iteration 11 : 12.724131871390405
cost at iteration 12 : 12.62212813732685
cost at iteration 13 : 12.521480488076158
cost at iteration 14 : 12.422160599690764
cost at iteration 15 : 12.324141094801222
cost at iteration 16 : 12.227395495926604
cost at iteration 17 : 12.131898181677869
cost at iteration 18 : 12.03762434565384
cost at iteration 19 : 11.944549957844472
cost at iteration 20 : 11.85265172836919
cost at iteration 21 : 11.761907073390157
cost at iteration 22 : 11.672294083053641
cost at iteration 23 : 11.583791491321508
cost at iteration 24 : 11.496378647565159
cost at iteration 25 : 11.410035489805589
cost at iteration 26 : 11.324742519486739
cost at iteration 27 : 11.240480777683652
cost at iteration 28 : 11.15723182264898
cost at iteration 29 : 11.074977708610167
cost at iteration 30 : 10.993700965735849
cost at iteration 31 : 10.913384581196372
cost at iteration 32 : 10.834011981246448
cost at iteration 33 : 10.755567014265564
cost at iteration 34 : 10.678033934694788
cost at iteration 35 : 10.60139738781274
cost at iteration 36 : 10.525642395299371
cost at iteration 37 : 10.4507543415359
cost at iteration 38 : 10.376718960597987
cost at iteration 39 : 10.30352232389744
cost at iteration 40 : 10.231150828433448
cost at iteration 41 : 10.159591185616444
cost at iteration 42 : 10.088830410629775
cost at iteration 43 : 10.018855812297367
cost at iteration 44 : 9.94965498342675
cost at iteration 45 : 9.881215791600304
cost at iteration 46 : 9.813526370387274
cost at iteration 47 : 9.746575110953774
cost at iteration 48 : 9.680350654046608
cost at iteration 49 : 9.614841882330133
cost at iteration 50 : 9.550037913055997
cost at iteration 51 : 9.485928091047798
cost at iteration 52 : 9.42250198198223
cost at iteration 53 : 9.359749365951204
cost at iteration 54 : 9.29766023128913
cost at iteration 55 : 9.236224768651516
cost at iteration 56 : 9.175433365331248
cost at iteration 57 : 9.115276599799587
cost at iteration 58 : 9.055745236461473
cost at iteration 59 : 8.996830220611923
cost at iteration 60 : 8.938522673585503
cost at iteration 61 : 8.880813888087435
cost at iteration 62 : 8.82369532369782
cost at iteration 63 : 8.76715860254074
cost at iteration 64 : 8.711195505109323
cost at iteration 65 : 8.655797966239959
cost at iteration 66 : 8.600958071228272
cost at iteration 67 : 8.546668052079939
cost at iteration 68 : 8.492920283890067
cost at iteration 69 : 8.439707281345619
cost at iteration 70 : 8.387021695344725
cost at iteration 71 : 8.33485630972766
cost at iteration 72 : 8.283204038114569
cost at iteration 73 : 8.232057920845817
cost at iteration 74 : 8.181411122019073
cost at iteration 75 : 8.131256926620479
cost at iteration 76 : 8.081588737744742
cost at iteration 77 : 8.032400073901128
cost at iteration 78 : 7.983684566401613
cost at iteration 79 : 7.9354359568272175
cost at iteration 80 : 7.887648094570333
cost at iteration 81 : 7.840314934449432
cost at iteration 82 : 7.7934305343928365
cost at iteration 83 : 7.746989053190266
cost at iteration 84 : 7.700984748307666
cost at iteration 85 : 7.655411973764373
cost at iteration 86 : 7.610265178069866
cost at iteration 87 : 7.565538902217125
cost at iteration 88 : 7.521227777731821
cost at iteration 89 : 7.4773265247739875
cost at iteration 90 : 7.433829950291335
cost at iteration 91 : 7.390732946221673
cost at iteration 92 : 7.348030487743153
cost at iteration 93 : 7.30571763157005
cost at iteration 94 : 7.263789514293292
cost at iteration 95 : 7.222241350763824
cost at iteration 96 : 7.181068432517166
cost at iteration 97 : 7.140266126237953
cost at iteration 98 : 7.099829872263646
cost at iteration 99 : 7.0597551831250165
cost at iteration 100 : 7.020037642123312
cost at iteration 101 : 6.980672901942279
cost at iteration 102 : 6.9416566832942195
cost at iteration 103 : 6.902984773598968
cost at iteration 104 : 6.864653025694407
cost at iteration 105 : 6.8266573565784245
cost at iteration 106 : 6.788993746180299
cost at iteration 107 : 6.751658236161375
cost at iteration 108 : 6.7146469287437895
cost at iteration 109 : 6.677955985566658
cost at iteration 110 : 6.641581626568691
cost at iteration 111 : 6.605520128896418
cost at iteration 112 : 6.569767825837985
cost at iteration 113 : 6.534321105780406
cost at iteration 114 : 6.499176411191092
cost at iteration 115 : 6.4643302376220655
cost at iteration 116 : 6.42977913273647
cost at iteration 117 : 6.395519695356891
cost at iteration 118 : 6.361548574534466
cost at iteration 119 : 6.327862468638919
cost at iteration 120 : 6.294458124468221
cost at iteration 121 : 6.261332336377701
cost at iteration 122 : 6.228481945428376
cost at iteration 123 : 6.195903838553306
cost at iteration 124 : 6.163594947741946
cost at iteration 125 : 6.131552249242174
cost at iteration 126 : 6.099772762778906
cost at iteration 127 : 6.068253550789795
cost at iteration 128 : 6.03699171767664
cost at iteration 129 : 6.005984409072565
cost at iteration 130 : 5.975228811124798
cost at iteration 131 : 5.944722149792106
cost at iteration 132 : 5.914461690156875
cost at iteration 133 : 5.88444473575141
cost at iteration 134 : 5.854668627898181
cost at iteration 135 : 5.825130745063481
cost at iteration 136 : 5.79582850222414
cost at iteration 137 : 5.766759350247391
cost at iteration 138 : 5.737920775283032
cost at iteration 139 : 5.70931029816814
cost at iteration 140 : 5.680925473843277
cost at iteration 141 : 5.652763890780821
cost at iteration 142 : 5.624823170424442
cost at iteration 143 : 5.59710096663976
cost at iteration 144 : 5.569594965175941
cost at iteration 145 : 5.542302883137909
cost at iteration 146 : 5.515222468468557
cost at iteration 147 : 5.488351499441873
cost at iteration 148 : 5.461687784165084
cost at iteration 149 : 5.435229160090971
cost at iteration 150 : 5.408973493539715
cost at iteration 151 : 5.3829186792295065
cost at iteration 152 : 5.357062639816603
cost at iteration 153 : 5.331403325444119
cost at iteration 154 : 5.305938713299228
cost at iteration 155 : 5.280666807179087
cost at iteration 156 : 5.255585637064755
cost at iteration 157 : 5.230693258703373
cost at iteration 158 : 5.205987753198138
cost at iteration 159 : 5.181467226605926
cost at iteration 160 : 5.157129809542473
cost at iteration 161 : 5.132973656795099
cost at iteration 162 : 5.1089969469424465
cost at iteration 163 : 5.085197881981353
cost at iteration 164 : 5.0615746869606255
cost at iteration 165 : 5.038125609621544
cost at iteration 166 : 5.014848920045193
cost at iteration 167 : 4.991742910305761
cost at iteration 168 : 4.968805894130814
cost at iteration 169 : 4.946036206567214
cost at iteration 170 : 4.923432203653479
cost at iteration 171 : 4.90099226209791
cost at iteration 172 : 4.878714778962587
cost at iteration 173 : 4.8565981713529744
cost at iteration 174 : 4.834640876113276
cost at iteration 175 : 4.812841349527101
cost at iteration 176 : 4.791198067023474
cost at iteration 177 : 4.769709522888272
cost at iteration 178 : 4.748374229980345
cost at iteration 179 : 4.727190719453457
cost at iteration 180 : 4.706157540482098
cost at iteration 181 : 4.685273259993144
cost at iteration 182 : 4.664536462401425
cost at iteration 183 : 4.643945749350507
cost at iteration 184 : 4.6234997394575155
cost at iteration 185 : 4.603197068062952
cost at iteration 186 : 4.583036386984256
cost at iteration 187 : 4.563016364274103
cost at iteration 188 : 4.543135683982712
cost at iteration 189 : 4.523393045924095
cost at iteration 190 : 4.50378716544673
cost at iteration 191 : 4.484316773207632
cost at iteration 192 : 4.464980614950813
cost at iteration 193 : 4.445777451289177
cost at iteration 194 : 4.426706057490039
cost at iteration 195 : 4.407765223264819
cost at iteration 196 : 4.388953752561654
cost at iteration 197 : 4.370270463361901
cost at iteration 198 : 4.351714187479962
cost at iteration 199 : 4.333283770366416
cost at iteration 200 : 4.3149780709145515
cost at iteration 201 : 4.296795961270009
cost at iteration 202 : 4.27873632664377
cost at iteration 203 : 4.260798065127959
cost at iteration 204 : 4.242980087514956
cost at iteration 205 : 4.225281317119616
cost at iteration 206 : 4.207700689603904
cost at iteration 207 : 4.190237152804947
cost at iteration 208 : 4.1728896665655855
cost at iteration 209 : 4.155657202567927
cost at iteration 210 : 4.138538744169364
cost at iteration 211 : 4.1215332862414815
cost at iteration 212 : 4.10463983501145
cost at iteration 213 : 4.087857407906072
cost at iteration 214 : 4.071185033398243
cost at iteration 215 : 4.054621750856027
cost at iteration 216 : 4.0381666103940494
cost at iteration 217 : 4.021818672727305
cost at iteration 218 : 4.005577009027234
cost at iteration 219 : 3.989440700780418
cost at iteration 220 : 3.973408839648983
cost at iteration 221 : 3.957480527333594
cost at iteration 222 : 3.941654875438743
cost at iteration 223 : 3.925931005339714
cost at iteration 224 : 3.9103080480520895
cost at iteration 225 : 3.894785144102904
cost at iteration 226 : 3.8793614434042225
cost at iteration 227 : 3.864036105128496
cost at iteration 228 : 3.848808297585798
cost at iteration 229 : 3.8336771981030195
cost at iteration 230 : 3.8186419929052002
cost at iteration 231 : 3.8037018769982414
cost at iteration 232 : 3.788856054053679
cost at iteration 233 : 3.7741037362954715
cost at iteration 234 : 3.759444144388007
cost at iteration 235 : 3.7448765073262873
cost at iteration 236 : 3.7304000623275106
cost at iteration 237 : 3.7160140547244733
cost at iteration 238 : 3.7017177378606054
cost at iteration 239 : 3.6875103729864156
cost at iteration 240 : 3.6733912291577493
cost at iteration 241 : 3.659359583135402
cost at iteration 242 : 3.6454147192863853
cost at iteration 243 : 3.631555929486598
cost at iteration 244 : 3.617782513024946
cost at iteration 245 : 3.6040937765090435
cost at iteration 246 : 3.590489033772125
cost at iteration 247 : 3.5769676057815403
cost at iteration 248 : 3.5635288205483513
cost at iteration 249 : 3.5501720130387353
cost at iteration 250 : 3.5368965250861537
cost at iteration 251 : 3.5237017053051995
cost at iteration 252 : 3.510586909006592
cost at iteration 253 : 3.4975514981134235
cost at iteration 254 : 3.4845948410786276
cost at iteration 255 : 3.4717163128037187
cost at iteration 256 : 3.458915294558712
cost at iteration 257 : 3.4461911739031206
cost at iteration 258 : 3.4335433446081187
cost at iteration 259 : 3.420971206580122
cost at iteration 260 : 3.4084741657849307
cost at iteration 261 : 3.396051634173521
cost at iteration 262 : 3.383703029608659
cost at iteration 263 : 3.3714277757923736
cost at iteration 264 : 3.359225302194955
cost at iteration 265 : 3.347095043984623
cost at iteration 266 : 3.3350364419580947
cost at iteration 267 : 3.32304894247264
cost at iteration 268 : 3.3111319973786593
cost at iteration 269 : 3.299285063953235
cost at iteration 270 : 3.287507604834982
cost at iteration 271 : 3.2757990879595726
cost at iteration 272 : 3.2641589864959597
cost at iteration 273 : 3.2525867787840115
cost at iteration 274 : 3.2410819482725657
cost at iteration 275 : 3.2296439834587516
cost at iteration 276 : 3.218272377827663
cost at iteration 277 : 3.206966629793401
cost at iteration 278 : 3.1957262426406157
cost at iteration 279 : 3.1845507244668863
cost at iteration 280 : 3.1734395881261572
cost at iteration 281 : 3.162392351172511
cost at iteration 282 : 3.1514085358052246
cost at iteration 283 : 3.1404876688140906
cost at iteration 284 : 3.129629281526015
cost at iteration 285 : 3.1188329097518332
cost at iteration 286 : 3.108098093734223
cost at iteration 287 : 3.0974243780961825
cost at iteration 288 : 3.086811311790229
cost at iteration 289 : 3.0762584480482973
cost at iteration 290 : 3.065765344332298
cost at iteration 291 : 3.0553315622854633
cost at iteration 292 : 3.044956667684129
cost at iteration 293 : 3.034640230390382
cost at iteration 294 : 3.024381824305271
cost at iteration 295 : 3.014181027322641
cost at iteration 296 : 3.004037421283487
cost at iteration 297 : 2.993950591931224
cost at iteration 298 : 2.9839201288671626
cost at iteration 299 : 2.9739456255068797
cost at iteration 300 : 2.9640266790370684
cost at iteration 301 : 2.95416289037284
cost at iteration 302 : 2.9443538641159943
cost at iteration 303 : 2.93459920851317
cost at iteration 304 : 2.9248985354152612
cost at iteration 305 : 2.9152514602369415
cost at iteration 306 : 2.905657601916772
cost at iteration 307 : 2.8961165828780073
cost at iteration 308 : 2.8866280289897066
cost at iteration 309 : 2.877191569528519
cost at iteration 310 : 2.867806837140803
cost at iteration 311 : 2.8584734678055432
cost at iteration 312 : 2.849191100797221
cost at iteration 313 : 2.839959378649823
cost at iteration 314 : 2.830777947120737
cost at iteration 315 : 2.8216464551555402
cost at iteration 316 : 2.81256455485293
cost at iteration 317 : 2.803531901430348
cost at iteration 318 : 2.794548153189873
cost at iteration 319 : 2.785612971484631
cost at iteration 320 : 2.7767260206856648
cost at iteration 321 : 2.7678869681491616
cost at iteration 322 : 2.7590954841840993
cost at iteration 323 : 2.750351242020376
cost at iteration 324 : 2.7416539177772665
cost at iteration 325 : 2.733003190432423
cost at iteration 326 : 2.7243987417909237
cost at iteration 327 : 2.715840256455236
cost at iteration 328 : 2.7073274217950845
cost at iteration 329 : 2.6988599279179004
cost at iteration 330 : 2.690437467639729
cost at iteration 331 : 2.6820597364562824
cost at iteration 332 : 2.673726432514631
cost at iteration 333 : 2.6654372565848696
cost at iteration 334 : 2.6571919120325873
cost at iteration 335 : 2.6489901047912983
cost at iteration 336 : 2.6408315433355103
cost at iteration 337 : 2.6327159386538064
cost at iteration 338 : 2.624643004222619
cost at iteration 339 : 2.6166124559800505
cost at iteration 340 : 2.608624012300176
cost at iteration 341 : 2.6006773939674965
cost at iteration 342 : 2.5927723241520804
cost at iteration 343 : 2.584908528384302
cost at iteration 344 : 2.5770857345308076
cost at iteration 345 : 2.56930367276997
cost at iteration 346 : 2.5615620755681605
cost at iteration 347 : 2.5538606776560253
cost at iteration 348 : 2.546199216005103
cost at iteration 349 : 2.5385774298049624
cost at iteration 350 : 2.5309950604402665
cost at iteration 351 : 2.5234518514682605
cost at iteration 352 : 2.515947548596616
cost at iteration 353 : 2.5084818996615112
cost at iteration 354 : 2.501054654605828
cost at iteration 355 : 2.4936655654577473
cost at iteration 356 : 2.4863143863096298
cost at iteration 357 : 2.4790008732970406
cost at iteration 358 : 2.471724784578061
cost at iteration 359 : 2.4644858803129277
cost at iteration 360 : 2.4572839226437084
cost at iteration 361 : 2.450118675674617
cost at iteration 362 : 2.4429899054519395
cost at iteration 363 : 2.435897379944866
cost at iteration 364 : 2.4288408690261813
cost at iteration 365 : 2.421820144453088
cost at iteration 366 : 2.414834979848614
cost at iteration 367 : 2.407885150682844
cost at iteration 368 : 2.4009704342548415
cost at iteration 369 : 2.3940906096741497
cost at iteration 370 : 2.387245457843142
cost at iteration 371 : 2.3804347614392123
cost at iteration 372 : 2.3736583048972406
cost at iteration 373 : 2.3669158743922702
cost at iteration 374 : 2.360207257822471
cost at iteration 375 : 2.353532244792174
cost at iteration 376 : 2.346890626595135
cost at iteration 377 : 2.3402821961981233
cost at iteration 378 : 2.3337067482244263
cost at iteration 379 : 2.3271640789378374
cost at iteration 380 : 2.3206539862266475
cost at iteration 381 : 2.3141762695878505
cost at iteration 382 : 2.307730730111632
cost at iteration 383 : 2.301317170465829
cost at iteration 384 : 2.294935394880756
cost at iteration 385 : 2.2885852091342054
cost at iteration 386 : 2.28226642053642
cost at iteration 387 : 2.2759788379154244
cost at iteration 388 : 2.269722271602483
cost at iteration 389 : 2.2634965334176997
cost at iteration 390 : 2.2573014366557596
cost at iteration 391 : 2.2511367960719033
cost at iteration 392 : 2.245002427867939
cost at iteration 393 : 2.2388981496785156
cost at iteration 394 : 2.2328237805576
cost at iteration 395 : 2.226779140964867
cost at iteration 396 : 2.22076405275258
cost at iteration 397 : 2.2147783391522573
cost at iteration 398 : 2.2088218247617903
cost at iteration 399 : 2.2028943355326214
cost at iteration 400 : 2.1969956987567536
cost at iteration 401 : 2.1911257430545645
cost at iteration 402 : 2.185284298362089
cost at iteration 403 : 2.179471195918811
cost at iteration 404 : 2.1736862682554667
cost at iteration 405 : 2.167929349182026
cost at iteration 406 : 2.162200273775836
cost at iteration 407 : 2.1564988783697827
cost at iteration 408 : 2.150825000540728
cost at iteration 409 : 2.1451784790978556
cost at iteration 410 : 2.1395591540714687
cost at iteration 411 : 2.133966866701642
cost at iteration 412 : 2.1284014594269918
cost at iteration 413 : 2.122862775873796
cost at iteration 414 : 2.1173506608450197
cost at iteration 415 : 2.1118649603095836
cost at iteration 416 : 2.106405521391615
cost at iteration 417 : 2.1009721923599995
cost at iteration 418 : 2.0955648226178654
cost at iteration 419 : 2.0901832626922348
cost at iteration 420 : 2.084827364223976
cost at iteration 421 : 2.07949697995747
cost at iteration 422 : 2.074191963730797
cost at iteration 423 : 2.0689121704657145
cost at iteration 424 : 2.063657456158043
cost at iteration 425 : 2.0584276778677353
cost at iteration 426 : 2.053222693709529
cost at iteration 427 : 2.0480423628433977
cost at iteration 428 : 2.0428865454651546
cost at iteration 429 : 2.037755102797126
cost at iteration 430 : 2.032647897079085
cost at iteration 431 : 2.0275647915591453
cost at iteration 432 : 2.0225056504846646
cost at iteration 433 : 2.017470339093534
cost at iteration 434 : 2.0124587236052056
cost at iteration 435 : 2.0074706712121477
cost at iteration 436 : 2.0025060500710534
cost at iteration 437 : 1.9975647292945102
cost at iteration 438 : 1.9926465789423775
cost at iteration 439 : 1.9877514700136314
cost at iteration 440 : 1.982879274437885
cost at iteration 441 : 1.9780298650674404
cost at iteration 442 : 1.9732031156690253
cost at iteration 443 : 1.9683989009158784
cost at iteration 444 : 1.9636170963797466
cost at iteration 445 : 1.9588575785232047
cost at iteration 446 : 1.9541202246916343
cost at iteration 447 : 1.9494049131058178
cost at iteration 448 : 1.9447115228541376
cost at iteration 449 : 1.9400399338851813
cost at iteration 450 : 1.9353900270002058
cost at iteration 451 : 1.930761683845777
cost at iteration 452 : 1.9261547869065354
cost at iteration 453 : 1.921569219497869
cost at iteration 454 : 1.9170048657589447
cost at iteration 455 : 1.9124616106453542
cost at iteration 456 : 1.9079393399224565
cost at iteration 457 : 1.9034379401581238
cost at iteration 458 : 1.898957298716071
cost at iteration 459 : 1.8944973037489519
cost at iteration 460 : 1.8900578441918245
cost at iteration 461 : 1.88563880975517
cost at iteration 462 : 1.8812400909186653
cost at iteration 463 : 1.8768615789243626
cost at iteration 464 : 1.8725031657704054
cost at iteration 465 : 1.8681647442045735
cost at iteration 466 : 1.8638462077179465
cost at iteration 467 : 1.8595474505386596
cost at iteration 468 : 1.8552683676256985
cost at iteration 469 : 1.8510088546627164
cost at iteration 470 : 1.846768808051967
cost at iteration 471 : 1.8425481249083493
cost at iteration 472 : 1.8383467030533396
cost at iteration 473 : 1.8341644410092341
cost at iteration 474 : 1.8300012379930357
cost at iteration 475 : 1.8258569939110172
cost at iteration 476 : 1.8217316093526965
cost at iteration 477 : 1.8176249855853683
cost at iteration 478 : 1.81353702454827
cost at iteration 479 : 1.8094676288472926
cost at iteration 480 : 1.8054167017491978
cost at iteration 481 : 1.8013841471763428
cost at iteration 482 : 1.7973698697012255
cost at iteration 483 : 1.7933737745411453
cost at iteration 484 : 1.7893957675528274
cost at iteration 485 : 1.785435755227312
cost at iteration 486 : 1.781493644684605
cost at iteration 487 : 1.7775693436687017
cost at iteration 488 : 1.7736627605423319
cost at iteration 489 : 1.7697738042820614
cost at iteration 490 : 1.76590238447314
cost at iteration 491 : 1.7620484113046821
cost at iteration 492 : 1.7582117955646677
cost at iteration 493 : 1.7543924486351927
cost at iteration 494 : 1.7505902824875663
cost at iteration 495 : 1.7468052096775428
cost at iteration 496 : 1.7430371433407568
cost at iteration 497 : 1.7392859971877888
cost at iteration 498 : 1.735551685499818
cost at iteration 499 : 1.7318341231238346
cost at iteration 500 : 1.7281332254681958
cost at iteration 501 : 1.7244489084981307
cost at iteration 502 : 1.7207810887312072
cost at iteration 503 : 1.7171296832329783
cost at iteration 504 : 1.7134946096127064
cost at iteration 505 : 1.7098757860188099
cost at iteration 506 : 1.7062731311347472
cost at iteration 507 : 1.7026865641747577
cost at iteration 508 : 1.6991160048795564
cost at iteration 509 : 1.6955613735122526
cost at iteration 510 : 1.6920225908542317
cost at iteration 511 : 1.688499578200949
cost at iteration 512 : 1.6849922573579899
cost at iteration 513 : 1.6815005506370084
cost at iteration 514 : 1.678024380851817
cost at iteration 515 : 1.6745636713143124
cost at iteration 516 : 1.6711183458306498
cost at iteration 517 : 1.667688328697402
cost at iteration 518 : 1.6642735446976848
cost at iteration 519 : 1.6608739190973945
cost at iteration 520 : 1.6574893776413469
cost at iteration 521 : 1.6541198465497196
cost at iteration 522 : 1.6507652525142187
cost at iteration 523 : 1.6474255226944374
cost at iteration 524 : 1.6441005847143233
cost at iteration 525 : 1.6407903666584605
cost at iteration 526 : 1.6374947970686191
cost at iteration 527 : 1.6342138049401824
cost at iteration 528 : 1.6309473197185778
cost at iteration 529 : 1.6276952712959678
cost at iteration 530 : 1.6244575900077578
cost at iteration 531 : 1.6212342066290346
cost at iteration 532 : 1.6180250523715307
cost at iteration 533 : 1.6148300588799223
cost at iteration 534 : 1.61164915822883
cost at iteration 535 : 1.6084822829193206
cost at iteration 536 : 1.6053293658757424
cost at iteration 537 : 1.6021903404425513
cost at iteration 538 : 1.5990651403810356
cost at iteration 539 : 1.595953699866223
cost at iteration 540 : 1.5928559534836677
cost at iteration 541 : 1.5897718362264412
cost at iteration 542 : 1.5867012834919412
cost at iteration 543 : 1.5836442310790193
cost at iteration 544 : 1.5806006151847636
cost at iteration 545 : 1.5775703724015537
cost at iteration 546 : 1.5745534397141807
cost at iteration 547 : 1.5715497544968733
cost at iteration 548 : 1.5685592545102658
cost at iteration 549 : 1.565581877898658
cost at iteration 550 : 1.562617563187081
cost at iteration 551 : 1.5596662492784263
cost at iteration 552 : 1.5567278754507556
cost at iteration 553 : 1.5538023813542918
cost at iteration 554 : 1.5508897070089045
cost at iteration 555 : 1.5479897928011683
cost at iteration 556 : 1.5451025794818032
cost at iteration 557 : 1.5422280081628288
cost at iteration 558 : 1.5393660203150048
cost at iteration 559 : 1.5365165577651203
cost at iteration 560 : 1.5336795626935251
cost at iteration 561 : 1.530854977631256
cost at iteration 562 : 1.528042745457633
cost at iteration 563 : 1.5252428093977621
cost at iteration 564 : 1.522455113019884
cost at iteration 565 : 1.519679600232826
cost at iteration 566 : 1.516916215283718
cost at iteration 567 : 1.514164902755255
cost at iteration 568 : 1.5114256075635313
cost at iteration 569 : 1.5086982749553746
cost at iteration 570 : 1.5059828505060466
cost at iteration 571 : 1.5032792801168928
cost at iteration 572 : 1.500587510012929
cost at iteration 573 : 1.4979074867404847
cost at iteration 574 : 1.4952391571649495
cost at iteration 575 : 1.4925824684683713
cost at iteration 576 : 1.489937368147271
cost at iteration 577 : 1.4873038040103368
cost at iteration 578 : 1.4846817241761814
cost at iteration 579 : 1.4820710770710444
cost at iteration 580 : 1.4794718114267855
cost at iteration 581 : 1.4768838762784446
cost at iteration 582 : 1.474307220962263
cost at iteration 583 : 1.4717417951134668
cost at iteration 584 : 1.4691875486641068
cost at iteration 585 : 1.4666444318409806
cost at iteration 586 : 1.4641123951635793
cost at iteration 587 : 1.461591389441834
cost at iteration 588 : 1.4590813657743456
cost at iteration 589 : 1.4565822755460305
cost at iteration 590 : 1.45409407042634
cost at iteration 591 : 1.4516167023670545
cost at iteration 592 : 1.4491501236004325
cost at iteration 593 : 1.4466942866372101
cost at iteration 594 : 1.4442491442645735
cost at iteration 595 : 1.4418146495442976
cost at iteration 596 : 1.4393907558107133
cost at iteration 597 : 1.436977416668975
cost at iteration 598 : 1.4345745859929717
cost at iteration 599 : 1.4321822179235317
=================Mesh Quality Check=================
max quality: 1.1832106186704068 min quality: 0.5059467464043369 mean quality: 0.8036611401632984
===================================================
non linear with p#
check_mesh_quality()
c = pose_deformation_equation(p=2)
gfX = GridFunction(VEC)
gfX.Set((x *(1-x)*y*(1-y),0))
gfsetOld = GridFunction(VEC)
gfset = GridFunction(VEC)
gfset.Set((0,0))
scene = Draw(gfset,mesh,"gfset")
SetVisualization (deformation=True)
optimize_shape(linear=False)
mesh.SetDeformation(gfset)
check_mesh_quality()
=================Mesh Quality Check=================
max quality: 1.1832106186704068 min quality: 0.5059467464043369 mean quality: 0.8036611401632984
===================================================
Using p = 2
Cost at initial design 13.962388102754339
cost at iteration 0 : 13.962388102754339
Newton iteration 0
err = 37.410916722694864
Newton iteration 1
err = 1.8869214813442775e-13
cost at iteration 1 : 13.735607075299352
Newton iteration 0
err = 0.5046122613881028
Newton iteration 1
err = 4.4254487315379895e-14
cost at iteration 2 : 13.526734084141825
Newton iteration 0
err = 0.4939055405164836
Newton iteration 1
err = 4.363456293335879e-14
cost at iteration 3 : 13.3225169520256
Newton iteration 0
err = 0.48353678713810977
Newton iteration 1
err = 4.3915394551139703e-14
cost at iteration 4 : 13.122814774368264
Newton iteration 0
err = 0.473492208258463
Newton iteration 1
err = 4.218799185213006e-14
cost at iteration 5 : 12.9274922406713
Newton iteration 0
err = 0.46375867035191304
Newton iteration 1
err = 4.150894460032659e-14
cost at iteration 6 : 12.736419356284364
Newton iteration 0
err = 0.4543236642196425
Newton iteration 1
err = 4.182856003489478e-14
cost at iteration 7 : 12.54947118072345
Newton iteration 0
err = 0.4451752719111788
Newton iteration 1
err = 4.132104602605138e-14
cost at iteration 8 : 12.366527581410036
Newton iteration 0
err = 0.4363021355767377
Newton iteration 1
err = 4.0760320972295e-14
cost at iteration 9 : 12.187473001786444
Newton iteration 0
err = 0.42769342812216765
Newton iteration 1
err = 3.829579513147409e-14
cost at iteration 10 : 12.0121962428425
Newton iteration 0
err = 0.4193388255528737
Newton iteration 1
err = 3.9107369256031866e-14
cost at iteration 11 : 11.840590257157793
Newton iteration 0
err = 0.4112284808975276
Newton iteration 1
err = 3.829863816750436e-14
cost at iteration 12 : 11.672551954633821
Newton iteration 0
err = 0.4033529996112698
Newton iteration 1
err = 3.8760547806858766e-14
cost at iteration 13 : 11.50798201915008
Newton iteration 0
err = 0.3957034163676085
Newton iteration 1
err = 3.735585281847013e-14
cost at iteration 14 : 11.34678473543316
Newton iteration 0
err = 0.3882711731502286
Newton iteration 1
err = 3.679358962415122e-14
cost at iteration 15 : 11.188867825481541
Newton iteration 0
err = 0.3810480985658445
Newton iteration 1
err = 3.6355290237819265e-14
cost at iteration 16 : 11.034142293936283
Newton iteration 0
err = 0.3740263883041259
Newton iteration 1
err = 3.6423255470679594e-14
cost at iteration 17 : 10.882522281828033
Newton iteration 0
err = 0.3671985866726246
Newton iteration 1
err = 3.5168145324956277e-14
cost at iteration 18 : 10.733924928177329
Newton iteration 0
err = 0.3605575691446813
Newton iteration 1
err = 3.618354054421274e-14
cost at iteration 19 : 10.588270238956433
Newton iteration 0
err = 0.3540965258578649
Newton iteration 1
err = 3.466003726418036e-14
cost at iteration 20 : 10.445480962958058
Newton iteration 0
err = 0.3478089460072846
Newton iteration 1
err = 3.417259683553717e-14
cost at iteration 21 : 10.305482474149086
Newton iteration 0
err = 0.3416886030808486
Newton iteration 1
err = 3.363225970136936e-14
cost at iteration 22 : 10.1682026601132
Newton iteration 0
err = 0.3357295408871752
Newton iteration 1
err = 3.3657680969529875e-14
cost at iteration 23 : 10.03357181621658
Newton iteration 0
err = 0.32992606033011274
Newton iteration 1
err = 3.306700354588007e-14
cost at iteration 24 : 9.901522545154027
Newton iteration 0
err = 0.3242727068864402
Newton iteration 1
err = 3.3252305797305016e-14
cost at iteration 25 : 9.77198966155669
Newton iteration 0
err = 0.3187642587465984
Newton iteration 1
err = 3.226146710550975e-14
cost at iteration 26 : 9.644910101363886
Newton iteration 0
err = 0.3133957155809668
Newton iteration 1
err = 3.1820854747058744e-14
cost at iteration 27 : 9.520222835680451
Newton iteration 0
err = 0.3081622878946944
Newton iteration 1
err = 3.0166350286928206e-14
cost at iteration 28 : 9.397868788861103
Newton iteration 0
err = 0.30305938694005086
Newton iteration 1
err = 3.069483372461462e-14
cost at iteration 29 : 9.277790760578783
Newton iteration 0
err = 0.29808261515309536
Newton iteration 1
err = 3.1562027199591346e-14
cost at iteration 30 : 9.159933351650135
Newton iteration 0
err = 0.2932277570861597
Newton iteration 1
err = 3.1170655957994946e-14
cost at iteration 31 : 9.04424289340649
Newton iteration 0
err = 0.28849077080878344
Newton iteration 1
err = 3.098536667357062e-14
cost at iteration 32 : 8.930667380411426
Newton iteration 0
err = 0.28386777975081473
Newton iteration 1
err = 3.0822062436365796e-14
cost at iteration 33 : 8.819156406339548
Newton iteration 0
err = 0.27935506496283785
Newton iteration 1
err = 3.05640745024032e-14
cost at iteration 34 : 8.709661102842102
Newton iteration 0
err = 0.2749490577729587
Newton iteration 1
err = 2.9802068828511664e-14
cost at iteration 35 : 8.602134081235981
Newton iteration 0
err = 0.27064633281544587
Newton iteration 1
err = 2.9897951777283047e-14
cost at iteration 36 : 8.496529376863588
Newton iteration 0
err = 0.26644360141381834
Newton iteration 1
err = 2.979879115117831e-14
cost at iteration 37 : 8.392802395979752
Newton iteration 0
err = 0.2623377052972754
Newton iteration 1
err = 2.974981838288552e-14
cost at iteration 38 : 8.290909865030997
Newton iteration 0
err = 0.2583256106333532
Newton iteration 1
err = 2.916544314229156e-14
cost at iteration 39 : 8.190809782200082
Newton iteration 0
err = 0.2544044023596154
Newton iteration 1
err = 2.897880130158268e-14
cost at iteration 40 : 8.092461371098164
Newton iteration 0
err = 0.250571278798772
Newton iteration 1
err = 2.9149188223265397e-14
cost at iteration 41 : 7.995825036490887
Newton iteration 0
err = 0.24682354654094163
Newton iteration 1
err = 2.8388198077377854e-14
cost at iteration 42 : 7.900862321955695
Newton iteration 0
err = 0.243158615580995
Newton iteration 1
err = 2.927281966432458e-14
cost at iteration 43 : 7.8075358693684205
Newton iteration 0
err = 0.23957399469454663
Newton iteration 1
err = 2.8065210185213697e-14
cost at iteration 44 : 7.715809380129594
Newton iteration 0
err = 0.23606728704329613
Newton iteration 1
err = 2.7923198285731088e-14
cost at iteration 45 : 7.62564757803949
Newton iteration 0
err = 0.2326361859949322
Newton iteration 1
err = 2.770632745826125e-14
cost at iteration 46 : 7.5370161737411125
Newton iteration 0
err = 0.22927847114741998
Newton iteration 1
err = 2.7424721947950033e-14
cost at iteration 47 : 7.449881830653638
Newton iteration 0
err = 0.22599200454819124
Newton iteration 1
err = 2.6656721630787922e-14
cost at iteration 48 : 7.364212132320758
Newton iteration 0
err = 0.222774727095234
Newton iteration 1
err = 2.781901017996575e-14
cost at iteration 49 : 7.279975551107125
Newton iteration 0
err = 0.21962465511355989
Newton iteration 1
err = 2.5895490588667782e-14
cost at iteration 50 : 7.197141418175964
Newton iteration 0
err = 0.21653987709637457
Newton iteration 1
err = 2.697390467807852e-14
cost at iteration 51 : 7.115679894686386
Newton iteration 0
err = 0.2135185506024936
Newton iteration 1
err = 2.6024578110089555e-14
cost at iteration 52 : 7.0355619441525
Newton iteration 0
err = 0.21055889930284774
Newton iteration 1
err = 2.5894851503522644e-14
cost at iteration 53 : 6.956759305908962
Newton iteration 0
err = 0.20765921016731087
Newton iteration 1
err = 2.6162624102862953e-14
cost at iteration 54 : 6.879244469630437
Newton iteration 0
err = 0.20481783078554475
Newton iteration 1
err = 2.472819737989273e-14
cost at iteration 55 : 6.802990650856605
Newton iteration 0
err = 0.2020331668146619
Newton iteration 1
err = 2.4954669047009874e-14
cost at iteration 56 : 6.727971767475329
Newton iteration 0
err = 0.19930367954724418
Newton iteration 1
err = 2.5392850327199726e-14
cost at iteration 57 : 6.654162417119806
Newton iteration 0
err = 0.19662788359360964
Newton iteration 1
err = 2.4883511309256957e-14
cost at iteration 58 : 6.58153785543873
Newton iteration 0
err = 0.19400434467260577
Newton iteration 1
err = 2.495309677577878e-14
cost at iteration 59 : 6.510073975198968
Newton iteration 0
err = 0.19143167750582743
Newton iteration 1
err = 2.417306722744785e-14
cost at iteration 60 : 6.439747286183503
Newton iteration 0
err = 0.18890854380886424
Newton iteration 1
err = 2.395625253251977e-14
cost at iteration 61 : 6.370534895849263
Newton iteration 0
err = 0.1864336503764679
Newton iteration 1
err = 2.4179326833341218e-14
cost at iteration 62 : 6.302414490711011
Newton iteration 0
err = 0.18400574725527524
Newton iteration 1
err = 2.344351739420432e-14
cost at iteration 63 : 6.235364318419066
Newton iteration 0
err = 0.18162362600094725
Newton iteration 1
err = 2.2853774338649303e-14
cost at iteration 64 : 6.169363170500706
Newton iteration 0
err = 0.17928611801448188
Newton iteration 1
err = 2.3180929241673046e-14
cost at iteration 65 : 6.104390365736333
Newton iteration 0
err = 0.1769920929551006
Newton iteration 1
err = 2.29884149808342e-14
cost at iteration 66 : 6.040425734142831
Newton iteration 0
err = 0.1747404572242098
Newton iteration 1
err = 2.2407926063751413e-14
cost at iteration 67 : 5.977449601538489
Newton iteration 0
err = 0.1725301525187306
Newton iteration 1
err = 2.206037414049852e-14
cost at iteration 68 : 5.915442774664087
Newton iteration 0
err = 0.17036015444878166
Newton iteration 1
err = 2.2592987813551075e-14
cost at iteration 69 : 5.854386526837393
Newton iteration 0
err = 0.16822947121793563
Newton iteration 1
err = 2.2332218979478733e-14
cost at iteration 70 : 5.7942625841178
Newton iteration 0
err = 0.16613714236173774
Newton iteration 1
err = 2.3162700265990732e-14
cost at iteration 71 : 5.735053111960767
Newton iteration 0
err = 0.16408223754294976
Newton iteration 1
err = 2.0775037205834536e-14
cost at iteration 72 : 5.676740702340764
Newton iteration 0
err = 0.16206385539922852
Newton iteration 1
err = 2.0964340905126647e-14
cost at iteration 73 : 5.619308361324378
Newton iteration 0
err = 0.16008112244253891
Newton iteration 1
err = 2.0825835946906057e-14
cost at iteration 74 : 5.562739497074374
Newton iteration 0
err = 0.15813319200569742
Newton iteration 1
err = 2.0811534547971275e-14
cost at iteration 75 : 5.5070179082677875
Newton iteration 0
err = 0.15621924323612166
Newton iteration 1
err = 2.0992040522753272e-14
cost at iteration 76 : 5.452127772910877
Newton iteration 0
err = 0.15433848013196255
Newton iteration 1
err = 2.0233013255686562e-14
cost at iteration 77 : 5.3980536375354795
Newton iteration 0
err = 0.15249013062120986
Newton iteration 1
err = 2.1085335798391568e-14
cost at iteration 78 : 5.344780406760929
Newton iteration 0
err = 0.15067344567927646
Newton iteration 1
err = 2.031073308090293e-14
cost at iteration 79 : 5.292293333208101
Newton iteration 0
err = 0.14888769848541278
Newton iteration 1
err = 2.01752988900952e-14
cost at iteration 80 : 5.240578007750145
Newton iteration 0
err = 0.14713218361408853
Newton iteration 1
err = 2.0531484264300172e-14
cost at iteration 81 : 5.189620350088669
Newton iteration 0
err = 0.14540621626164266
Newton iteration 1
err = 2.0550330125789433e-14
cost at iteration 82 : 5.139406599640823
Newton iteration 0
err = 0.14370913150480172
Newton iteration 1
err = 1.9492766330431437e-14
cost at iteration 83 : 5.089923306726742
Newton iteration 0
err = 0.14204028359084603
Newton iteration 1
err = 1.988326561839075e-14
cost at iteration 84 : 5.041157324045196
Newton iteration 0
err = 0.14039904525719915
Newton iteration 1
err = 2.0092245142869654e-14
cost at iteration 85 : 4.993095798426067
Newton iteration 0
err = 0.1387848070794409
Newton iteration 1
err = 1.9595886835325316e-14
cost at iteration 86 : 4.945726162850187
Newton iteration 0
err = 0.13719697684644117
Newton iteration 1
err = 1.9221262179063385e-14
cost at iteration 87 : 4.899036128725188
Newton iteration 0
err = 0.13563497896062468
Newton iteration 1
err = 1.9633910333095527e-14
cost at iteration 88 : 4.853013678409138
Newton iteration 0
err = 0.1340982538638475
Newton iteration 1
err = 2.0035138301752925e-14
cost at iteration 89 : 4.80764705797135
Newton iteration 0
err = 0.13258625748560676
Newton iteration 1
err = 1.89933144326837e-14
cost at iteration 90 : 4.762924770182638
Newton iteration 0
err = 0.13109846071423814
Newton iteration 1
err = 1.9084821237066645e-14
cost at iteration 91 : 4.7188355677261935
Newton iteration 0
err = 0.12963434888950617
Newton iteration 1
err = 1.908089189428801e-14
cost at iteration 92 : 4.6753684466205785
Newton iteration 0
err = 0.1281934213149521
Newton iteration 1
err = 1.9183504421377437e-14
cost at iteration 93 : 4.632512639848179
Newton iteration 0
err = 0.12677519079062152
Newton iteration 1
err = 1.9424698294404588e-14
cost at iteration 94 : 4.590257611180182
Newton iteration 0
err = 0.12537918316337374
Newton iteration 1
err = 1.8321507594471017e-14
cost at iteration 95 : 4.548593049192574
Newton iteration 0
err = 0.12400493689564202
Newton iteration 1
err = 1.8534113542856063e-14
cost at iteration 96 : 4.507508861465317
Newton iteration 0
err = 0.12265200265086788
Newton iteration 1
err = 1.8902169458220728e-14
cost at iteration 97 : 4.466995168958002
Newton iteration 0
err = 0.12131994289514317
Newton iteration 1
err = 1.8254690088613713e-14
cost at iteration 98 : 4.427042300557031
Newton iteration 0
err = 0.12000833151442203
Newton iteration 1
err = 1.8367255808737127e-14
cost at iteration 99 : 4.3876407877865615
Newton iteration 0
err = 0.11871675344618546
Newton iteration 1
err = 1.7992305711889887e-14
cost at iteration 100 : 4.348781359678971
Newton iteration 0
err = 0.11744480432582347
Newton iteration 1
err = 1.884271770104458e-14
cost at iteration 101 : 4.310454937798219
Newton iteration 0
err = 0.11619209014621458
Newton iteration 1
err = 1.7319665579270483e-14
cost at iteration 102 : 4.2726526314112485
Newton iteration 0
err = 0.11495822692982222
Newton iteration 1
err = 1.774913269200734e-14
cost at iteration 103 : 4.235365732802734
Newton iteration 0
err = 0.1137428404140657
Newton iteration 1
err = 1.7829726332219462e-14
cost at iteration 104 : 4.1985857127274455
Newton iteration 0
err = 0.11254556574790923
Newton iteration 1
err = 1.7405712286700443e-14
cost at iteration 105 : 4.162304215995716
Newton iteration 0
err = 0.11136604719987281
Newton iteration 1
err = 1.756172201034928e-14
cost at iteration 106 : 4.126513057188112
Newton iteration 0
err = 0.1102039378767274
Newton iteration 1
err = 1.717652900010677e-14
cost at iteration 107 : 4.09120421649457
Newton iteration 0
err = 0.109058899453365
Newton iteration 1
err = 1.7852769246321188e-14
cost at iteration 108 : 4.0563698356731175
Newton iteration 0
err = 0.10793060191130138
Newton iteration 1
err = 1.746576537520936e-14
cost at iteration 109 : 4.0220022141255996
Newton iteration 0
err = 0.10681872328794062
Newton iteration 1
err = 1.742867705322332e-14
cost at iteration 110 : 3.9880938050851227
Newton iteration 0
err = 0.10572294943417872
Newton iteration 1
err = 1.676228308777403e-14
cost at iteration 111 : 3.954637211912643
Newton iteration 0
err = 0.10464297378118169
Newton iteration 1
err = 1.669486191878093e-14
cost at iteration 112 : 3.9216251844979313
Newton iteration 0
err = 0.10357849711504515
Newton iteration 1
err = 1.689193527853468e-14
cost at iteration 113 : 3.889050615763077
Newton iteration 0
err = 0.10252922736021403
Newton iteration 1
err = 1.7766500066388993e-14
cost at iteration 114 : 3.856906538263273
Newton iteration 0
err = 0.10149487936977013
Newton iteration 1
err = 1.6725253937344247e-14
cost at iteration 115 : 3.8251861208834375
Newton iteration 0
err = 0.10047517472378296
Newton iteration 1
err = 1.6238910010123773e-14
cost at iteration 116 : 3.793882665626482
Newton iteration 0
err = 0.09946984153432696
Newton iteration 1
err = 1.6395255617308234e-14
cost at iteration 117 : 3.7629896044906963
Newton iteration 0
err = 0.09847861425743942
Newton iteration 1
err = 1.5727642125647482e-14
cost at iteration 118 : 3.732500496433299
Newton iteration 0
err = 0.09750123351158986
Newton iteration 1
err = 1.568034726830412e-14
cost at iteration 119 : 3.7024090244177703
Newton iteration 0
err = 0.09653744590232963
Newton iteration 1
err = 1.6037676393262092e-14
cost at iteration 120 : 3.672708992541458
Newton iteration 0
err = 0.09558700385309694
Newton iteration 1
err = 1.6033422571830485e-14
cost at iteration 121 : 3.643394323241782
Newton iteration 0
err = 0.0946496654415284
Newton iteration 1
err = 1.6751463785325726e-14
cost at iteration 122 : 3.614459054578219
Newton iteration 0
err = 0.09372519424160372
Newton iteration 1
err = 1.5615148123479615e-14
cost at iteration 123 : 3.5858973375876624
Newton iteration 0
err = 0.09281335917098991
Newton iteration 1
err = 1.6064218745597007e-14
cost at iteration 124 : 3.557703433710705
Newton iteration 0
err = 0.0919139343430389
Newton iteration 1
err = 1.60652103742555e-14
cost at iteration 125 : 3.5298717122874064
Newton iteration 0
err = 0.09102669892492495
Newton iteration 1
err = 1.5693592792509927e-14
cost at iteration 126 : 3.502396648119054
Newton iteration 0
err = 0.09015143699881792
Newton iteration 1
err = 1.5608244602409452e-14
cost at iteration 127 : 3.475272819095457
Newton iteration 0
err = 0.08928793742898121
Newton iteration 1
err = 1.578886090981116e-14
cost at iteration 128 : 3.4484949038845554
Newton iteration 0
err = 0.08843599373253572
Newton iteration 1
err = 1.5148270305905005e-14
cost at iteration 129 : 3.4220576796829736
Newton iteration 0
err = 0.08759540395469989
Newton iteration 1
err = 1.5070675741909975e-14
cost at iteration 130 : 3.3959560200255274
Newton iteration 0
err = 0.08676597054772177
Newton iteration 1
err = 1.5696967849700388e-14
cost at iteration 131 : 3.3701848926520497
Newton iteration 0
err = 0.08594750025410273
Newton iteration 1
err = 1.500713706640934e-14
cost at iteration 132 : 3.3447393574297646
Newton iteration 0
err = 0.08513980399342592
Newton iteration 1
err = 1.516075667593954e-14
cost at iteration 133 : 3.319614564329166
Newton iteration 0
err = 0.08434269675235208
Newton iteration 1
err = 1.509716673750858e-14
cost at iteration 134 : 3.2948057514526643
Newton iteration 0
err = 0.0835559974789482
Newton iteration 1
err = 1.515976880432024e-14
cost at iteration 135 : 3.2703082431133517
Newton iteration 0
err = 0.08277952897943444
Newton iteration 1
err = 1.5225933734585967e-14
cost at iteration 136 : 3.2461174479632677
Newton iteration 0
err = 0.08201311781874558
Newton iteration 1
err = 1.4631550580765774e-14
cost at iteration 137 : 3.2222288571694198
Newton iteration 0
err = 0.08125659422386022
Newton iteration 1
err = 1.442773747625835e-14
cost at iteration 138 : 3.198638042635913
Newton iteration 0
err = 0.08050979199041783
Newton iteration 1
err = 1.3838707331107502e-14
cost at iteration 139 : 3.175340655271209
Newton iteration 0
err = 0.07977254839172224
Newton iteration 1
err = 1.4623982618863782e-14
cost at iteration 140 : 3.1523324232990286
Newton iteration 0
err = 0.07904470409115459
Newton iteration 1
err = 1.453715757656939e-14
cost at iteration 141 : 3.1296091506115546
Newton iteration 0
err = 0.07832610305660725
Newton iteration 1
err = 1.4349225516809257e-14
cost at iteration 142 : 3.107166715164028
Newton iteration 0
err = 0.07761659247786462
Newton iteration 1
err = 1.431283348533617e-14
cost at iteration 143 : 3.0850010674090296
Newton iteration 0
err = 0.07691602268648064
Newton iteration 1
err = 1.3941155891119381e-14
cost at iteration 144 : 3.063108228769916
Newton iteration 0
err = 0.07622424707793228
Newton iteration 1
err = 1.4540268373360353e-14
cost at iteration 145 : 3.041484290151681
Newton iteration 0
err = 0.07554112203586151
Newton iteration 1
err = 1.37695543053199e-14
cost at iteration 146 : 3.0201254104887747
Newton iteration 0
err = 0.07486650685916968
Newton iteration 1
err = 1.4098390029979167e-14
cost at iteration 147 : 2.9990278153282937
Newton iteration 0
err = 0.07420026369051555
Newton iteration 1
err = 1.352876507212318e-14
cost at iteration 148 : 2.9781877954480436
Newton iteration 0
err = 0.07354225744750793
Newton iteration 1
err = 1.3943105355031476e-14
cost at iteration 149 : 2.9576017055079973
Newton iteration 0
err = 0.07289235575546647
Newton iteration 1
err = 1.3870780890945365e-14
cost at iteration 150 : 2.9372659627345326
Newton iteration 0
err = 0.0722504288823758
Newton iteration 1
err = 1.443055284328891e-14
cost at iteration 151 : 2.917177045636685
Newton iteration 0
err = 0.07161634967566578
Newton iteration 1
err = 1.3677476582190203e-14
cost at iteration 152 : 2.897331492752833
Newton iteration 0
err = 0.07098999350057675
Newton iteration 1
err = 1.3224128636860656e-14
cost at iteration 153 : 2.8777259014281062
Newton iteration 0
err = 0.07037123818071718
Newton iteration 1
err = 1.3239281864153681e-14
cost at iteration 154 : 2.8583569266202615
Newton iteration 0
err = 0.06975996393977388
Newton iteration 1
err = 1.3366927330370637e-14
cost at iteration 155 : 2.839221279734652
Newton iteration 0
err = 0.06915605334540897
Newton iteration 1
err = 1.3027927935681592e-14
cost at iteration 156 : 2.8203157274863218
Newton iteration 0
err = 0.06855939125407003
Newton iteration 1
err = 1.3285580588119147e-14
cost at iteration 157 : 2.8016370907894563
Newton iteration 0
err = 0.06796986475810478
Newton iteration 1
err = 1.2929030979614705e-14
cost at iteration 158 : 2.783182243672615
Newton iteration 0
err = 0.06738736313372283
Newton iteration 1
err = 1.3232730305129686e-14
cost at iteration 159 : 2.7649481122195065
Newton iteration 0
err = 0.066811777790548
Newton iteration 1
err = 1.3319236468666535e-14
cost at iteration 160 : 2.746931673534717
Newton iteration 0
err = 0.06624300222282171
Newton iteration 1
err = 1.2468109900809142e-14
cost at iteration 161 : 2.7291299547333683
Newton iteration 0
err = 0.06568093196158645
Newton iteration 1
err = 1.2489869650503688e-14
cost at iteration 162 : 2.711540031954288
Newton iteration 0
err = 0.0651254645283338
Newton iteration 1
err = 1.208742928932281e-14
cost at iteration 163 : 2.694159029395877
Newton iteration 0
err = 0.06457649938984061
Newton iteration 1
err = 1.2145135797743977e-14
cost at iteration 164 : 2.6769841183745475
Newton iteration 0
err = 0.06403393791446017
Newton iteration 1
err = 1.2379827987839742e-14
cost at iteration 165 : 2.660012516404277
Newton iteration 0
err = 0.06349768332898342
Newton iteration 1
err = 1.1496994039170966e-14
cost at iteration 166 : 2.6432414862978044
Newton iteration 0
err = 0.06296764067752605
Newton iteration 1
err = 1.1795401191088332e-14
cost at iteration 167 : 2.6266683352878166
Newton iteration 0
err = 0.06244371678062058
Newton iteration 1
err = 1.1963292218033105e-14
cost at iteration 168 : 2.610290414168519
Newton iteration 0
err = 0.061925820195903813
Newton iteration 1
err = 1.19608416839995e-14
cost at iteration 169 : 2.5941051164563564
Newton iteration 0
err = 0.06141386117996197
Newton iteration 1
err = 1.127006141057721e-14
cost at iteration 170 : 2.578109877569801
Newton iteration 0
err = 0.060907751650451314
Newton iteration 1
err = 1.1685477546714063e-14
cost at iteration 171 : 2.562302174027734
Newton iteration 0
err = 0.060407405150140196
Newton iteration 1
err = 1.1870829775855695e-14
cost at iteration 172 : 2.5466795226655083
Newton iteration 0
err = 0.05991273681108933
Newton iteration 1
err = 1.1201846772749181e-14
cost at iteration 173 : 2.5312394798687916
Newton iteration 0
err = 0.05942366332028936
Newton iteration 1
err = 1.1982155855904087e-14
cost at iteration 174 : 2.515979640824251
Newton iteration 0
err = 0.058940102885943535
Newton iteration 1
err = 1.1523888841153917e-14
cost at iteration 175 : 2.5008976387870554
Newton iteration 0
err = 0.05846197520441013
Newton iteration 1
err = 1.1056230714241185e-14
cost at iteration 176 : 2.48599114436449
Newton iteration 0
err = 0.05798920142857383
Newton iteration 1
err = 1.1365928150840587e-14
cost at iteration 177 : 2.471257864815464
Newton iteration 0
err = 0.0575217041364129
Newton iteration 1
err = 1.1287123139500858e-14
cost at iteration 178 : 2.4566955433651088
Newton iteration 0
err = 0.05705940730062341
Newton iteration 1
err = 1.1209242436060726e-14
cost at iteration 179 : 2.4423019585348325
Newton iteration 0
err = 0.056602236259090166
Newton iteration 1
err = 1.1447131095083275e-14
cost at iteration 180 : 2.428074923486592
Newton iteration 0
err = 0.05615011768587242
Newton iteration 1
err = 1.1720226700249744e-14
cost at iteration 181 : 2.4140122853817636
Newton iteration 0
err = 0.055702979563093016
Newton iteration 1
err = 1.150235688536667e-14
cost at iteration 182 : 2.400111924753508
Newton iteration 0
err = 0.05526075115335434
Newton iteration 1
err = 1.1764797866550177e-14
cost at iteration 183 : 2.3863717548932826
Newton iteration 0
err = 0.05482336297315851
Newton iteration 1
err = 1.0966104965182652e-14
cost at iteration 184 : 2.372789721249856
Newton iteration 0
err = 0.054390746766426724
Newton iteration 1
err = 1.3886464952224779e-14
cost at iteration 185 : 2.359363800841905
Newton iteration 0
err = 0.053962835479271824
Newton iteration 1
err = 1.1007439335482281e-14
cost at iteration 186 : 2.346092001682705
Newton iteration 0
err = 0.0535395632349775
Newton iteration 1
err = 1.0910820269684659e-14
cost at iteration 187 : 2.3329723622172556
Newton iteration 0
err = 0.05312086530969652
Newton iteration 1
err = 1.1712142769008469e-14
cost at iteration 188 : 2.320002950771466
Newton iteration 0
err = 0.05270667810877874
Newton iteration 1
err = 1.0992817642057991e-14
cost at iteration 189 : 2.3071818650128186
Newton iteration 0
err = 0.052296939143755224
Newton iteration 1
err = 1.0568932073922597e-14
cost at iteration 190 : 2.2945072314226214
Newton iteration 0
err = 0.05189158700946137
Newton iteration 1
err = 1.07806256627934e-14
cost at iteration 191 : 2.281977204779218
Newton iteration 0
err = 0.051490561362246716
Newton iteration 1
err = 1.0602168922705796e-14
cost at iteration 192 : 2.269589967652175
Newton iteration 0
err = 0.05109380289845302
Newton iteration 1
err = 1.092523542843782e-14
cost at iteration 193 : 2.2573437299068657
Newton iteration 0
err = 0.05070125333329303
Newton iteration 1
err = 1.1203151189071521e-14
cost at iteration 194 : 2.2452367282195502
Newton iteration 0
err = 0.050312855380289456
Newton iteration 1
err = 1.1161432876548931e-14
cost at iteration 195 : 2.233267225602631
Newton iteration 0
err = 0.04992855273148209
Newton iteration 1
err = 1.079607730955441e-14
cost at iteration 196 : 2.2214335109393466
Newton iteration 0
err = 0.04954829003755455
Newton iteration 1
err = 1.0625111162586442e-14
cost at iteration 197 : 2.2097338985286035
Newton iteration 0
err = 0.049172012889007605
Newton iteration 1
err = 1.0826366877414574e-14
cost at iteration 198 : 2.1981667276386285
Newton iteration 0
err = 0.048799667797245996
Newton iteration 1
err = 1.0645821841132533e-14
cost at iteration 199 : 2.186730362070373
Newton iteration 0
err = 0.048431202176540865
Newton iteration 1
err = 1.0287885543556046e-14
cost at iteration 200 : 2.175423189729347
Newton iteration 0
err = 0.04806656432609454
Newton iteration 1
err = 1.0135520412293773e-14
cost at iteration 201 : 2.164243622206382
Newton iteration 0
err = 0.04770570341269443
Newton iteration 1
err = 1.0571375978387038e-14
cost at iteration 202 : 2.1531900943670617
Newton iteration 0
err = 0.04734856945355085
Newton iteration 1
err = 1.055526034513479e-14
cost at iteration 203 : 2.142261063949316
Newton iteration 0
err = 0.04699511329996871
Newton iteration 1
err = 1.052770201069458e-14
cost at iteration 204 : 2.1314550111690376
Newton iteration 0
err = 0.04664528662084675
Newton iteration 1
err = 1.027684725761626e-14
cost at iteration 205 : 2.120770438333968
Newton iteration 0
err = 0.04629904188682957
Newton iteration 1
err = 1.0186749844140971e-14
cost at iteration 206 : 2.110205869465021
Newton iteration 0
err = 0.045956332354955345
Newton iteration 1
err = 1.0220591737986611e-14
cost at iteration 207 : 2.0997598499254235
Newton iteration 0
err = 0.04561711205321229
Newton iteration 1
err = 1.0480088207610617e-14
cost at iteration 208 : 2.0894309460571536
Newton iteration 0
err = 0.04528133576600432
Newton iteration 1
err = 1.0308991655980906e-14
cost at iteration 209 : 2.079217744824544
Newton iteration 0
err = 0.04494895901937785
Newton iteration 1
err = 1.0055049983706952e-14
cost at iteration 210 : 2.0691188534652327
Newton iteration 0
err = 0.0446199380670377
Newton iteration 1
err = 1.0267110724762877e-14
cost at iteration 211 : 2.059132899147489
Newton iteration 0
err = 0.04429422987630728
Newton iteration 1
err = 1.04209815132137e-14
cost at iteration 212 : 2.0492585286351055
Newton iteration 0
err = 0.043971792114818345
Newton iteration 1
err = 9.799951514141768e-15
cost at iteration 213 : 2.0394944079579895
Newton iteration 0
err = 0.043652583136896204
Newton iteration 1
err = 1.0173706033227935e-14
cost at iteration 214 : 2.0298392220899695
Newton iteration 0
err = 0.04333656197090666
Newton iteration 1
err = 9.923679075918152e-15
cost at iteration 215 : 2.020291674632348
Newton iteration 0
err = 0.043023688306315984
Newton iteration 1
err = 1.0310929777184594e-14
cost at iteration 216 : 2.010850487504017
Newton iteration 0
err = 0.042713922481594425
Newton iteration 1
err = 9.86972975455831e-15
cost at iteration 217 : 2.001514400637221
Newton iteration 0
err = 0.042407225471618556
Newton iteration 1
err = 1.0228614191292706e-14
cost at iteration 218 : 1.9922821716796657
Newton iteration 0
err = 0.04210355887614734
Newton iteration 1
err = 9.571072201271871e-15
cost at iteration 219 : 1.9831525757021429
Newton iteration 0
err = 0.04180288490822253
Newton iteration 1
err = 9.961330197063449e-15
cost at iteration 220 : 1.974124404911329
Newton iteration 0
err = 0.04150516638236573
Newton iteration 1
err = 9.264800212373807e-15
cost at iteration 221 : 1.9651964683693186
Newton iteration 0
err = 0.04121036670397088
Newton iteration 1
err = 9.611506619708045e-15
cost at iteration 222 : 1.9563675917170427
Newton iteration 0
err = 0.04091844985799478
Newton iteration 1
err = 9.758634754387624e-15
cost at iteration 223 : 1.9476366169041894
Newton iteration 0
err = 0.040629380398669875
Newton iteration 1
err = 9.581598069818273e-15
cost at iteration 224 : 1.939002401923543
Newton iteration 0
err = 0.04034312343856659
Newton iteration 1
err = 9.84477477369681e-15
cost at iteration 225 : 1.9304638205506306
Newton iteration 0
err = 0.040059644638886754
Newton iteration 1
err = 9.56323526775078e-15
cost at iteration 226 : 1.9220197620883508
Newton iteration 0
err = 0.03977891019917632
Newton iteration 1
err = 9.406645109163497e-15
cost at iteration 227 : 1.9136691311161722
Newton iteration 0
err = 0.03950088684755654
Newton iteration 1
err = 9.550944053726653e-15
cost at iteration 228 : 1.9054108472444122
Newton iteration 0
err = 0.0392255418311759
Newton iteration 1
err = 9.509808278021026e-15
cost at iteration 229 : 1.8972438448727744
Newton iteration 0
err = 0.03895284290677868
Newton iteration 1
err = 9.543664633182348e-15
cost at iteration 230 : 1.8891670729536754
Newton iteration 0
err = 0.03868275833162445
Newton iteration 1
err = 9.392355158156096e-15
cost at iteration 231 : 1.8811794947598897
Newton iteration 0
err = 0.038415256854325106
Newton iteration 1
err = 9.19469262035993e-15
cost at iteration 232 : 1.8732800876563072
Newton iteration 0
err = 0.038150307706107354
Newton iteration 1
err = 9.353927558094818e-15
cost at iteration 233 : 1.8654678428763722
Newton iteration 0
err = 0.03788788059219876
Newton iteration 1
err = 9.364163939199157e-15
cost at iteration 234 : 1.857741765302336
Newton iteration 0
err = 0.037627945683314824
Newton iteration 1
err = 9.000083050109269e-15
cost at iteration 235 : 1.850100873249579
Newton iteration 0
err = 0.03737047360740218
Newton iteration 1
err = 9.161811975728316e-15
cost at iteration 236 : 1.8425441982549575
Newton iteration 0
err = 0.037115435441485496
Newton iteration 1
err = 9.054763167763328e-15
cost at iteration 237 : 1.83507078486923
Newton iteration 0
err = 0.03686280270382594
Newton iteration 1
err = 9.091233634866697e-15
cost at iteration 238 : 1.8276796904528105
Newton iteration 0
err = 0.0366125473459175
Newton iteration 1
err = 9.358502731421925e-15
cost at iteration 239 : 1.820369984975784
Newton iteration 0
err = 0.036364641745107965
Newton iteration 1
err = 9.02674909724339e-15
cost at iteration 240 : 1.8131407508211574
Newton iteration 0
err = 0.03611905869692779
Newton iteration 1
err = 9.243068807294554e-15
cost at iteration 241 : 1.8059910825919052
Newton iteration 0
err = 0.03587577140784971
Newton iteration 1
err = 9.120282639027886e-15
cost at iteration 242 : 1.798920086921396
Newton iteration 0
err = 0.0356347534880395
Newton iteration 1
err = 8.907400695484304e-15
cost at iteration 243 : 1.791926882287321
Newton iteration 0
err = 0.035395978944427274
Newton iteration 1
err = 9.116842223798918e-15
cost at iteration 244 : 1.7850105988289253
Newton iteration 0
err = 0.035159422173700335
Newton iteration 1
err = 9.776157539729783e-15
cost at iteration 245 : 1.7781703781674647
Newton iteration 0
err = 0.03492505795557075
Newton iteration 1
err = 9.040293483995416e-15
cost at iteration 246 : 1.7714053732299668
Newton iteration 0
err = 0.03469286144619063
Newton iteration 1
err = 8.854989880527818e-15
cost at iteration 247 : 1.7647147480762109
Newton iteration 0
err = 0.034462808171624636
Newton iteration 1
err = 8.698419961648744e-15
cost at iteration 248 : 1.75809767772856
Newton iteration 0
err = 0.03423487402141801
Newton iteration 1
err = 8.821326787631025e-15
cost at iteration 249 : 1.7515533480049628
Newton iteration 0
err = 0.03400903524246598
Newton iteration 1
err = 8.79791680765444e-15
cost at iteration 250 : 1.7450809553549127
Newton iteration 0
err = 0.03378526843285553
Newton iteration 1
err = 8.754063350856767e-15
cost at iteration 251 : 1.738679706698288
Newton iteration 0
err = 0.03356355053579524
Newton iteration 1
err = 9.600195281082473e-15
cost at iteration 252 : 1.73234881926702
Newton iteration 0
err = 0.033343858833751294
Newton iteration 1
err = 8.62249460571416e-15
cost at iteration 253 : 1.7260875204496053
Newton iteration 0
err = 0.03312617094281825
Newton iteration 1
err = 8.913762178379869e-15
cost at iteration 254 : 1.7198950476381134
Newton iteration 0
err = 0.03291046480679048
Newton iteration 1
err = 8.2545569559591e-15
cost at iteration 255 : 1.7137706480782255
Newton iteration 0
err = 0.03269671869188567
Newton iteration 1
err = 9.309215277754185e-15
cost at iteration 256 : 1.707713578721636
Newton iteration 0
err = 0.03248491118110761
Newton iteration 1
err = 8.742127662734057e-15
cost at iteration 257 : 1.7017231060809876
Newton iteration 0
err = 0.03227502116901751
Newton iteration 1
err = 8.869286219743448e-15
cost at iteration 258 : 1.6957985060875984
Newton iteration 0
err = 0.03206702785637755
Newton iteration 1
err = 8.875557449390921e-15
cost at iteration 259 : 1.689939063951304
Newton iteration 0
err = 0.03186091074517047
Newton iteration 1
err = 8.533144650340584e-15
cost at iteration 260 : 1.684144074022993
Newton iteration 0
err = 0.03165664963344586
Newton iteration 1
err = 8.558013309963941e-15
cost at iteration 261 : 1.678412839659345
Newton iteration 0
err = 0.03145422461037338
Newton iteration 1
err = 8.676350916376211e-15
cost at iteration 262 : 1.6727446730899125
Newton iteration 0
err = 0.03125361605145948
Newton iteration 1
err = 8.390347964326866e-15
cost at iteration 263 : 1.667138895286503
Newton iteration 0
err = 0.031054804613713714
Newton iteration 1
err = 8.5360704219807e-15
cost at iteration 264 : 1.6615948358349548
Newton iteration 0
err = 0.03085777123106561
Newton iteration 1
err = 8.265887461544966e-15
cost at iteration 265 : 1.656111832808618
Newton iteration 0
err = 0.030662497109662193
Newton iteration 1
err = 8.331200950328902e-15
cost at iteration 266 : 1.6506892326446538
Newton iteration 0
err = 0.030468963723464165
Newton iteration 1
err = 9.095647199825274e-15
cost at iteration 267 : 1.6453263900217099
Newton iteration 0
err = 0.030277152809807112
Newton iteration 1
err = 8.14405257501335e-15
cost at iteration 268 : 1.640022667740313
Newton iteration 0
err = 0.03008704636505082
Newton iteration 1
err = 8.296304556449326e-15
cost at iteration 269 : 1.6347774366049042
Newton iteration 0
err = 0.029898626640313064
Newton iteration 1
err = 8.458986810645183e-15
cost at iteration 270 : 1.6295900753078367
Newton iteration 0
err = 0.029711876137422703
Newton iteration 1
err = 8.351842636506941e-15
cost at iteration 271 : 1.6244599703155504
Newton iteration 0
err = 0.029526777604571237
Newton iteration 1
err = 8.139835503796986e-15
cost at iteration 272 : 1.6193865157564487
Newton iteration 0
err = 0.029343314032563468
Newton iteration 1
err = 7.885385115291997e-15
cost at iteration 273 : 1.6143691133107332
Newton iteration 0
err = 0.029161468650673324
Newton iteration 1
err = 8.33638732154356e-15
cost at iteration 274 : 1.6094071721021148
Newton iteration 0
err = 0.02898122492285875
Newton iteration 1
err = 8.269826229598571e-15
cost at iteration 275 : 1.604500108591133
Newton iteration 0
err = 0.028802566543947545
Newton iteration 1
err = 8.058224588326022e-15
cost at iteration 276 : 1.5996473464705097
Newton iteration 0
err = 0.02862547743589727
Newton iteration 1
err = 9.045948371734885e-15
cost at iteration 277 : 1.5948483165619856
Newton iteration 0
err = 0.02844994174405817
Newton iteration 1
err = 8.176254215301458e-15
cost at iteration 278 : 1.5901024567150657
Newton iteration 0
err = 0.02827594383370081
Newton iteration 1
err = 8.240728358040853e-15
cost at iteration 279 : 1.5854092117072285
Newton iteration 0
err = 0.028103468286348143
Newton iteration 1
err = 8.391086223953496e-15
cost at iteration 280 : 1.580768033145884
Newton iteration 0
err = 0.027932499896412267
Newton iteration 1
err = 8.557067579853645e-15
cost at iteration 281 : 1.5761783793719932
Newton iteration 0
err = 0.02776302366763554
Newton iteration 1
err = 7.737330875101282e-15
cost at iteration 282 : 1.571639715365117
Newton iteration 0
err = 0.027595024809915766
Newton iteration 1
err = 8.125209069618717e-15
cost at iteration 283 : 1.5671515126500843
Newton iteration 0
err = 0.027428488735891616
Newton iteration 1
err = 8.026627691704212e-15
cost at iteration 284 : 1.5627132492051914
Newton iteration 0
err = 0.027263401057694883
Newton iteration 1
err = 7.78343667748406e-15
cost at iteration 285 : 1.55832440937178
Newton iteration 0
err = 0.027099747583913612
Newton iteration 1
err = 7.982513934726638e-15
cost at iteration 286 : 1.5539844837655796
Newton iteration 0
err = 0.026937514316254164
Newton iteration 1
err = 7.973795284697961e-15
cost at iteration 287 : 1.549692969189044
Newton iteration 0
err = 0.02677668744670462
Newton iteration 1
err = 7.48007636062614e-15
cost at iteration 288 : 1.5454493685452957
Newton iteration 0
err = 0.026617253354342545
Newton iteration 1
err = 7.500205247459698e-15
cost at iteration 289 : 1.5412531907537554
Newton iteration 0
err = 0.02645919860240618
Newton iteration 1
err = 7.650441019754136e-15
cost at iteration 290 : 1.5371039506663529
Newton iteration 0
err = 0.02630250993554459
Newton iteration 1
err = 7.583058943236667e-15
cost at iteration 291 : 1.5330011689861016
Newton iteration 0
err = 0.026147174276699287
Newton iteration 1
err = 7.664425478913346e-15
cost at iteration 292 : 1.528944372185966
Newton iteration 0
err = 0.025993178724475798
Newton iteration 1
err = 7.575901710533628e-15
cost at iteration 293 : 1.52493309242979
Newton iteration 0
err = 0.025840510550363294
Newton iteration 1
err = 8.050068489905822e-15
cost at iteration 294 : 1.5209668674939323
Newton iteration 0
err = 0.025689157195950975
Newton iteration 1
err = 7.882893167940142e-15
cost at iteration 295 : 1.5170452406905635
Newton iteration 0
err = 0.025539106270275315
Newton iteration 1
err = 7.366113700053655e-15
cost at iteration 296 : 1.5131677607918528
Newton iteration 0
err = 0.025390345547250854
Newton iteration 1
err = 7.540186047190117e-15
cost at iteration 297 : 1.5093339819555065
Newton iteration 0
err = 0.025242862963100743
Newton iteration 1
err = 7.20292447477048e-15
cost at iteration 298 : 1.505543463651462
Newton iteration 0
err = 0.02509664661367184
Newton iteration 1
err = 7.435562457103207e-15
cost at iteration 299 : 1.5017957705896965
Newton iteration 0
err = 0.024951684752153428
Newton iteration 1
err = 7.584314153186356e-15
cost at iteration 300 : 1.4980904726491875
Newton iteration 0
err = 0.024807965786468485
Newton iteration 1
err = 7.623211761175861e-15
cost at iteration 301 : 1.4944271448079767
Newton iteration 0
err = 0.024665478276945887
Newton iteration 1
err = 7.534415310744736e-15
cost at iteration 302 : 1.4908053670742984
Newton iteration 0
err = 0.024524210933878448
Newton iteration 1
err = 6.948741615770957e-15
cost at iteration 303 : 1.487224724418807
Newton iteration 0
err = 0.024384152615336496
Newton iteration 1
err = 7.647145393369557e-15
cost at iteration 304 : 1.4836848067078714
Newton iteration 0
err = 0.024245292324699415
Newton iteration 1
err = 7.06546624158858e-15
cost at iteration 305 : 1.4801852086377938
Newton iteration 0
err = 0.02410761920848713
Newton iteration 1
err = 7.53735802647261e-15
cost at iteration 306 : 1.4767255296702095
Newton iteration 0
err = 0.02397112255423721
Newton iteration 1
err = 7.283318303577721e-15
cost at iteration 307 : 1.4733053739682442
Newton iteration 0
err = 0.023835791788056202
Newton iteration 1
err = 7.333658832930132e-15
cost at iteration 308 : 1.4699243503340491
Newton iteration 0
err = 0.02370161647281716
Newton iteration 1
err = 7.160977070905204e-15
cost at iteration 309 : 1.4665820721467613
Newton iteration 0
err = 0.023568586305822124
Newton iteration 1
err = 7.180923128721638e-15
cost at iteration 310 : 1.463278157301896
Newton iteration 0
err = 0.023436691116734555
Newton iteration 1
err = 7.452080377693612e-15
cost at iteration 311 : 1.4600122281512942
Newton iteration 0
err = 0.023305920865706663
Newton iteration 1
err = 7.469386492665296e-15
cost at iteration 312 : 1.456783911444226
Newton iteration 0
err = 0.023176265641149874
Newton iteration 1
err = 7.05563947776917e-15
cost at iteration 313 : 1.453592838269358
Newton iteration 0
err = 0.023047715657897542
Newton iteration 1
err = 6.9702950709158405e-15
cost at iteration 314 : 1.4504386439974586
Newton iteration 0
err = 0.022920261255288895
Newton iteration 1
err = 7.067074135898416e-15
cost at iteration 315 : 1.4473209682250379
Newton iteration 0
err = 0.02279389289509277
Newton iteration 1
err = 7.359202477254275e-15
cost at iteration 316 : 1.4442394547189723
Newton iteration 0
err = 0.022668601159807355
Newton iteration 1
err = 6.9680686706595326e-15
cost at iteration 317 : 1.4411937513617123
Newton iteration 0
err = 0.022544376750721073
Newton iteration 1
err = 7.818626527741644e-15
cost at iteration 318 : 1.4381835100975193
Newton iteration 0
err = 0.02242121048603246
Newton iteration 1
err = 6.868513216054526e-15
cost at iteration 319 : 1.4352083868794991
Newton iteration 0
err = 0.02229909329921659
Newton iteration 1
err = 6.9981385995630444e-15
cost at iteration 320 : 1.4322680416171023
Newton iteration 0
err = 0.02217801623707167
Newton iteration 1
err = 6.809975678551943e-15
cost at iteration 321 : 1.4293621381249835
Newton iteration 0
err = 0.022057970458135755
Newton iteration 1
err = 7.039104936825374e-15
cost at iteration 322 : 1.4264903440721564
Newton iteration 0
err = 0.02193894723085272
Newton iteration 1
err = 6.874515387267434e-15
cost at iteration 323 : 1.4236523309318958
Newton iteration 0
err = 0.021820937931965635
Newton iteration 1
err = 6.8161884695552965e-15
cost at iteration 324 : 1.4208477739328556
Newton iteration 0
err = 0.02170393404488064
Newton iteration 1
err = 6.946049796171718e-15
cost at iteration 325 : 1.4180763520103388
Newton iteration 0
err = 0.021587927157880753
Newton iteration 1
err = 7.020565052423586e-15
cost at iteration 326 : 1.4153377477586393
Newton iteration 0
err = 0.02147290896275567
Newton iteration 1
err = 6.906553854623312e-15
cost at iteration 327 : 1.4126316473839353
Newton iteration 0
err = 0.02135887125293732
Newton iteration 1
err = 7.086819471574583e-15
cost at iteration 328 : 1.4099577406579167
Newton iteration 0
err = 0.021245805922218336
Newton iteration 1
err = 6.789760753033439e-15
cost at iteration 329 : 1.4073157208721752
Newton iteration 0
err = 0.021133704963023296
Newton iteration 1
err = 6.626435611997565e-15
cost at iteration 330 : 1.4047052847930914
Newton iteration 0
err = 0.021022560464935443
Newton iteration 1
err = 7.079174804108771e-15
cost at iteration 331 : 1.402126132617442
Newton iteration 0
err = 0.020912364613301416
Newton iteration 1
err = 7.596051467262828e-15
cost at iteration 332 : 1.3995779679288523
Newton iteration 0
err = 0.02080310968764293
Newton iteration 1
err = 6.485422564876489e-15
cost at iteration 333 : 1.3970604976545442
Newton iteration 0
err = 0.02069478806022344
Newton iteration 1
err = 6.2259718238701865e-15
cost at iteration 334 : 1.394573432022856
Newton iteration 0
err = 0.02058739219473534
Newton iteration 1
err = 6.4775913418277475e-15
cost at iteration 335 : 1.3921164845214715
Newton iteration 0
err = 0.020480914644782924
Newton iteration 1
err = 6.7253109615051366e-15
cost at iteration 336 : 1.3896893718562897
Newton iteration 0
err = 0.0203753480525228
Newton iteration 1
err = 6.418206625657877e-15
cost at iteration 337 : 1.3872918139105188
Newton iteration 0
err = 0.020270685147378193
Newton iteration 1
err = 6.531838168622758e-15
cost at iteration 338 : 1.3849235337048131
Newton iteration 0
err = 0.020166918744588175
Newton iteration 1
err = 6.591120509306196e-15
cost at iteration 339 : 1.3825842573577827
Newton iteration 0
err = 0.020064041743884434
Newton iteration 1
err = 6.7771590611222375e-15
cost at iteration 340 : 1.3802737140468042
Newton iteration 0
err = 0.019962047128375904
Newton iteration 1
err = 6.417868904569248e-15
cost at iteration 341 : 1.3779916359700686
Newton iteration 0
err = 0.019860927963027058
Newton iteration 1
err = 6.6640964111746985e-15
cost at iteration 342 : 1.3757377583084465
Newton iteration 0
err = 0.019760677393508196
Newton iteration 1
err = 6.40250117841675e-15
cost at iteration 343 : 1.3735118191882918
Newton iteration 0
err = 0.019661288644921077
Newton iteration 1
err = 6.814382063579104e-15
cost at iteration 344 : 1.3713135596446613
Newton iteration 0
err = 0.01956275502068379
Newton iteration 1
err = 6.185956616481845e-15
cost at iteration 345 : 1.3691427235852867
Newton iteration 0
err = 0.01946506990107182
Newton iteration 1
err = 6.267419111187413e-15
cost at iteration 346 : 1.3669990577544866
Newton iteration 0
err = 0.01936822674232329
Newton iteration 1
err = 6.534867497082318e-15
cost at iteration 347 : 1.364882311698321
Newton iteration 0
err = 0.01927221907513693
Newton iteration 1
err = 6.756303440668222e-15
cost at iteration 348 : 1.3627922377295518
Newton iteration 0
err = 0.019177040503908484
Newton iteration 1
err = 6.33259031543098e-15
cost at iteration 349 : 1.36072859089365
Newton iteration 0
err = 0.019082684705179116
Newton iteration 1
err = 6.229324677804126e-15
cost at iteration 350 : 1.3586911289349992
Newton iteration 0
err = 0.018989145426779777
Newton iteration 1
err = 6.4555105148366385e-15
cost at iteration 351 : 1.3566796122634308
Newton iteration 0
err = 0.018896416486651316
Newton iteration 1
err = 6.2273718513285985e-15
cost at iteration 352 : 1.3546938039217367
Newton iteration 0
err = 0.018804491771654424
Newton iteration 1
err = 6.3108726943394105e-15
cost at iteration 353 : 1.3527334695529398
Newton iteration 0
err = 0.018713365236733608
Newton iteration 1
err = 7.233120106997587e-15
cost at iteration 354 : 1.3507983773687198
Newton iteration 0
err = 0.018623030903518654
Newton iteration 1
err = 6.404668438865247e-15
cost at iteration 355 : 1.3488882981177963
Newton iteration 0
err = 0.018533482859504503
Newton iteration 1
err = 5.969444456880012e-15
cost at iteration 356 : 1.3470030050548205
Newton iteration 0
err = 0.018444715257035078
Newton iteration 1
err = 6.122444158604868e-15
cost at iteration 357 : 1.3451422739099563
Newton iteration 0
err = 0.018356722312132482
Newton iteration 1
err = 5.892235518465987e-15
cost at iteration 358 : 1.3433058828586582
Newton iteration 0
err = 0.018269498303551086
Newton iteration 1
err = 6.173646213943639e-15
cost at iteration 359 : 1.3414936124918218
Newton iteration 0
err = 0.018183037571826315
Newton iteration 1
err = 5.810765128175706e-15
cost at iteration 360 : 1.339705245786437
Newton iteration 0
err = 0.01809733451835703
Newton iteration 1
err = 5.989111736222585e-15
cost at iteration 361 : 1.337940568076825
Newton iteration 0
err = 0.018012383604240397
Newton iteration 1
err = 5.945514451287301e-15
cost at iteration 362 : 1.3361993670259684
Newton iteration 0
err = 0.017928179349425
Newton iteration 1
err = 6.0940408806413e-15
cost at iteration 363 : 1.3344814325972343
Newton iteration 0
err = 0.01784471633191367
Newton iteration 1
err = 5.851655961315726e-15
cost at iteration 364 : 1.332786557026774
Newton iteration 0
err = 0.017761989186568128
Newton iteration 1
err = 6.001268321584067e-15
cost at iteration 365 : 1.3311145347960986
Newton iteration 0
err = 0.017679992604465528
Newton iteration 1
err = 6.369560087098946e-15
cost at iteration 366 : 1.3294651626050498
Newton iteration 0
err = 0.01759872133172245
Newton iteration 1
err = 6.168926911786663e-15
cost at iteration 367 : 1.3278382393450783
Newton iteration 0
err = 0.017518170168895984
Newton iteration 1
err = 5.90289735754596e-15
cost at iteration 368 : 1.3262335660729636
Newton iteration 0
err = 0.017438333969929377
Newton iteration 1
err = 6.174426911674652e-15
cost at iteration 369 : 1.324650945985037
Newton iteration 0
err = 0.017359207641296143
Newton iteration 1
err = 5.9948213860156785e-15
cost at iteration 370 : 1.3230901843914091
Newton iteration 0
err = 0.017280786141140293
Newton iteration 1
err = 5.79212216110854e-15
cost at iteration 371 : 1.321551088690667
Newton iteration 0
err = 0.017203064478691313
Newton iteration 1
err = 5.874699895966492e-15
cost at iteration 372 : 1.3200334683454054
Newton iteration 0
err = 0.0171260377128738
Newton iteration 1
err = 6.0502263532975194e-15
cost at iteration 373 : 1.3185371348568984
Newton iteration 0
err = 0.017049700952271974
Newton iteration 1
err = 6.223870960019374e-15
cost at iteration 374 : 1.3170619017417504
Newton iteration 0
err = 0.01697404935348647
Newton iteration 1
err = 6.043434207164425e-15
cost at iteration 375 : 1.3156075845072315
Newton iteration 0
err = 0.016899078121071037
Newton iteration 1
err = 5.943256304764971e-15
cost at iteration 376 : 1.314174000628201
Newton iteration 0
err = 0.016824782506145792
Newton iteration 1
err = 6.844394568230098e-15
cost at iteration 377 : 1.312760969523343
Newton iteration 0
err = 0.016751157806269276
Newton iteration 1
err = 6.8495683193937266e-15
cost at iteration 378 : 1.3113683125328477
Newton iteration 0
err = 0.01667819936392945
Newton iteration 1
err = 5.923953893807216e-15
cost at iteration 379 : 1.3099958528950681
Newton iteration 0
err = 0.016605902566455295
Newton iteration 1
err = 6.011100033495453e-15
cost at iteration 380 : 1.3086434157245443
Newton iteration 0
err = 0.01653426284487813
Newton iteration 1
err = 5.8669752003181395e-15
cost at iteration 381 : 1.307310827989854
Newton iteration 0
err = 0.016463275673377385
Newton iteration 1
err = 5.812744711672109e-15
cost at iteration 382 : 1.3059979184918404
Newton iteration 0
err = 0.01639293656853612
Newton iteration 1
err = 6.09228836665353e-15
cost at iteration 383 : 1.3047045178421026
Newton iteration 0
err = 0.016323241088515208
Newton iteration 1
err = 5.837020910162792e-15
cost at iteration 384 : 1.303430458441755
Newton iteration 0
err = 0.01625418483256933
Newton iteration 1
err = 5.9253014304586185e-15
cost at iteration 385 : 1.3021755744606078
Newton iteration 0
err = 0.016185763440107907
Newton iteration 1
err = 5.864652114231764e-15
cost at iteration 386 : 1.3009397018163489
Newton iteration 0
err = 0.016117972590316394
Newton iteration 1
err = 6.177141125724994e-15
cost at iteration 387 : 1.2997226781543867
Newton iteration 0
err = 0.016050808001115413
Newton iteration 1
err = 6.179931120366633e-15
cost at iteration 388 : 1.2985243428274889
Newton iteration 0
err = 0.01598426542889171
Newton iteration 1
err = 5.848886668157695e-15
cost at iteration 389 : 1.2973445368761696
Newton iteration 0
err = 0.015918340667485306
Newton iteration 1
err = 5.981263418983771e-15
cost at iteration 390 : 1.2961831030087758
Newton iteration 0
err = 0.015853029547899314
Newton iteration 1
err = 6.107898601837273e-15
cost at iteration 391 : 1.2950398855825656
Newton iteration 0
err = 0.015788327937311185
Newton iteration 1
err = 6.319386942874652e-15
cost at iteration 392 : 1.2939147305842142
Newton iteration 0
err = 0.01572423173879727
Newton iteration 1
err = 5.707167852131487e-15
cost at iteration 393 : 1.292807485611282
Newton iteration 0
err = 0.01566073689047814
Newton iteration 1
err = 5.822553018982528e-15
cost at iteration 394 : 1.2917179998534898
Newton iteration 0
err = 0.01559783936491258
Newton iteration 1
err = 5.570841773342408e-15
cost at iteration 395 : 1.290646124074247
Newton iteration 0
err = 0.015535535168737662
Newton iteration 1
err = 5.529300688326916e-15
cost at iteration 396 : 1.2895917105927512
Newton iteration 0
err = 0.015473820341799003
Newton iteration 1
err = 5.883189559447527e-15
cost at iteration 397 : 1.288554613265949
Newton iteration 0
err = 0.015412690956789921
Newton iteration 1
err = 5.8028068681042725e-15
cost at iteration 398 : 1.2875346874709503
Newton iteration 0
err = 0.015352143118500092
Newton iteration 1
err = 5.995074323107259e-15
cost at iteration 399 : 1.2865317900875928
Newton iteration 0
err = 0.01529217296345702
Newton iteration 1
err = 5.992024095095047e-15
cost at iteration 400 : 1.2855457794811842
Newton iteration 0
err = 0.01523277665917158
Newton iteration 1
err = 5.690705783984125e-15
cost at iteration 401 : 1.2845765154855913
Newton iteration 0
err = 0.015173950403808362
Newton iteration 1
err = 6.134458827159431e-15
cost at iteration 402 : 1.2836238593863976
Newton iteration 0
err = 0.015115690425328728
Newton iteration 1
err = 6.042845223990916e-15
cost at iteration 403 : 1.282687673904337
Newton iteration 0
err = 0.015057992981438884
Newton iteration 1
err = 5.897063683198683e-15
cost at iteration 404 : 1.2817678231791796
Newton iteration 0
err = 0.015000854358606215
Newton iteration 1
err = 5.966105664034391e-15
cost at iteration 405 : 1.2808641727531263
Newton iteration 0
err = 0.0149442708718234
Newton iteration 1
err = 5.659846629868067e-15
cost at iteration 406 : 1.2799765895554072
Newton iteration 0
err = 0.01488823886407971
Newton iteration 1
err = 5.842821371365607e-15
cost at iteration 407 : 1.2791049418861868
Newton iteration 0
err = 0.014832754705723657
Newton iteration 1
err = 5.5748200987261026e-15
cost at iteration 408 : 1.278249099401141
Newton iteration 0
err = 0.014777814794148999
Newton iteration 1
err = 5.690878976044175e-15
cost at iteration 409 : 1.277408933096219
Newton iteration 0
err = 0.014723415553248423
Newton iteration 1
err = 5.8441405195437625e-15
cost at iteration 410 : 1.2765843152922904
Newton iteration 0
err = 0.01466955343293476
Newton iteration 1
err = 5.45332165413058e-15
cost at iteration 411 : 1.2757751196203786
Newton iteration 0
err = 0.014616224908662075
Newton iteration 1
err = 5.5056489564586646e-15
cost at iteration 412 : 1.2749812210068245
Newton iteration 0
err = 0.014563426481126194
Newton iteration 1
err = 5.722646853019058e-15
cost at iteration 413 : 1.2742024956588198
Newton iteration 0
err = 0.01451115467551712
Newton iteration 1
err = 5.6274982105043215e-15
cost at iteration 414 : 1.2734388210497436
Newton iteration 0
err = 0.014459406041409398
Newton iteration 1
err = 6.202802125666374e-15
cost at iteration 415 : 1.2726900759052038
Newton iteration 0
err = 0.014408177152097254
Newton iteration 1
err = 5.770546981147802e-15
cost at iteration 416 : 1.2719561401888462
Newton iteration 0
err = 0.014357464604424672
Newton iteration 1
err = 5.793726147004763e-15
cost at iteration 417 : 1.2712368950887096
Newton iteration 0
err = 0.014307265017978453
Newton iteration 1
err = 5.948249531102655e-15
cost at iteration 418 : 1.2705322230031904
Newton iteration 0
err = 0.014257575035155875
Newton iteration 1
err = 5.549556453513066e-15
cost at iteration 419 : 1.269842007527918
Newton iteration 0
err = 0.014208391320224072
Newton iteration 1
err = 5.879297618825031e-15
cost at iteration 420 : 1.269166133441925
Newton iteration 0
err = 0.014159710559761859
Newton iteration 1
err = 5.7452734608892574e-15
cost at iteration 421 : 1.268504486695106
Newton iteration 0
err = 0.014111529461078945
Newton iteration 1
err = 6.3896739378319985e-15
cost at iteration 422 : 1.2678569543944442
Newton iteration 0
err = 0.014063844753122319
Newton iteration 1
err = 6.171689722003069e-15
cost at iteration 423 : 1.267223424791868
Newton iteration 0
err = 0.014016653185020015
Newton iteration 1
err = 5.366509694857744e-15
cost at iteration 424 : 1.266603787270908
Newton iteration 0
err = 0.013969951526594223
Newton iteration 1
err = 5.762464790329768e-15
cost at iteration 425 : 1.2659979323346957
Newton iteration 0
err = 0.013923736567252171
Newton iteration 1
err = 6.186034317632932e-15
cost at iteration 426 : 1.2654057515931374
Newton iteration 0
err = 0.013878005116458519
Newton iteration 1
err = 5.8697645829827826e-15
cost at iteration 427 : 1.2648271377511338
Newton iteration 0
err = 0.013832754002271762
Newton iteration 1
err = 5.336173758415826e-15
cost at iteration 428 : 1.264261984595809
Newton iteration 0
err = 0.01378798007244223
Newton iteration 1
err = 5.503303333346419e-15
cost at iteration 429 : 1.2637101869853793
Newton iteration 0
err = 0.013743680192794048
Newton iteration 1
err = 5.5886762406567975e-15
cost at iteration 430 : 1.263171640836775
Newton iteration 0
err = 0.013699851247552672
Newton iteration 1
err = 5.451590079248483e-15
cost at iteration 431 : 1.262646243114162
Newton iteration 0
err = 0.013656490139066701
Newton iteration 1
err = 6.1455988648969106e-15
cost at iteration 432 : 1.2621338918174283
Newton iteration 0
err = 0.013613593787267599
Newton iteration 1
err = 5.738012305029323e-15
cost at iteration 433 : 1.261634485970965
Newton iteration 0
err = 0.013571159129389066
Newton iteration 1
err = 5.778859492840211e-15
cost at iteration 434 : 1.261147925612145
Newton iteration 0
err = 0.013529183119917122
Newton iteration 1
err = 5.923486417032145e-15
cost at iteration 435 : 1.260674111780469
Newton iteration 0
err = 0.013487662730018233
Newton iteration 1
err = 5.559852202221101e-15
cost at iteration 436 : 1.2602129465064973
Newton iteration 0
err = 0.013446594947409438
Newton iteration 1
err = 5.6562147892789476e-15
cost at iteration 437 : 1.259764332801013
Newton iteration 0
err = 0.013405976776103916
Newton iteration 1
err = 5.703763191652325e-15
cost at iteration 438 : 1.259328174644364
Newton iteration 0
err = 0.01336580523616169
Newton iteration 1
err = 5.817800866376618e-15
cost at iteration 439 : 1.258904376975956
Newton iteration 0
err = 0.013326077363218803
Newton iteration 1
err = 6.3299934867846884e-15
cost at iteration 440 : 1.2584928456836433
Newton iteration 0
err = 0.013286790208554361
Newton iteration 1
err = 5.5787674553857555e-15
cost at iteration 441 : 1.2580934875935779
Newton iteration 0
err = 0.013247940838603353
Newton iteration 1
err = 5.875665028954993e-15
cost at iteration 442 : 1.257706210459827
Newton iteration 0
err = 0.013209526335043256
Newton iteration 1
err = 6.308463397658376e-15
cost at iteration 443 : 1.2573309229547052
Newton iteration 0
err = 0.013171543793982148
Newton iteration 1
err = 5.359965503124107e-15
cost at iteration 444 : 1.2569675346583193
Newton iteration 0
err = 0.013133990326247156
Newton iteration 1
err = 5.288284530563216e-15
cost at iteration 445 : 1.2566159560489947
Newton iteration 0
err = 0.013096863057138083
Newton iteration 1
err = 5.414375694774223e-15
cost at iteration 446 : 1.2562760984936177
Newton iteration 0
err = 0.013060159125939493
Newton iteration 1
err = 5.366710623576179e-15
cost at iteration 447 : 1.2559478742379209
Newton iteration 0
err = 0.01302387568592076
Newton iteration 1
err = 5.828671748546536e-15
cost at iteration 448 : 1.2556311963971858
Newton iteration 0
err = 0.012988009903911482
Newton iteration 1
err = 5.6671927280742446e-15
cost at iteration 449 : 1.2553259789464075
Newton iteration 0
err = 0.012952558960567578
Newton iteration 1
err = 5.560721988244578e-15
cost at iteration 450 : 1.2550321367116926
Newton iteration 0
err = 0.012917520049776685
Newton iteration 1
err = 5.534973803749219e-15
cost at iteration 451 : 1.254749585360689
Newton iteration 0
err = 0.012882890378465092
Newton iteration 1
err = 5.1384203831663886e-15
cost at iteration 452 : 1.2544782413936248
Newton iteration 0
err = 0.012848667166694545
Newton iteration 1
err = 5.523569549800466e-15
cost at iteration 453 : 1.2542180221344141
Newton iteration 0
err = 0.012814847647228254
Newton iteration 1
err = 5.545068107144237e-15
cost at iteration 454 : 1.2539688457217764
Newton iteration 0
err = 0.012781429065740818
Newton iteration 1
err = 5.625002827229345e-15
cost at iteration 455 : 1.2537306311006549
Newton iteration 0
err = 0.012748408680023575
Newton iteration 1
err = 5.382164647065205e-15
cost at iteration 456 : 1.2535032980134027
Newton iteration 0
err = 0.012715783760563001
Newton iteration 1
err = 5.626101902185502e-15
cost at iteration 457 : 1.253286766991362
Newton iteration 0
err = 0.012683551589885943
Newton iteration 1
err = 5.297123799678775e-15
cost at iteration 458 : 1.2530809593465972
Newton iteration 0
err = 0.012651709462571898
Newton iteration 1
err = 5.247514321301302e-15
cost at iteration 459 : 1.2528857971632268
Newton iteration 0
err = 0.012620254685198413
Newton iteration 1
err = 5.1864529907606165e-15
cost at iteration 460 : 1.2527012032895979
Newton iteration 0
err = 0.012589184575982121
Newton iteration 1
err = 5.929417633523389e-15
cost at iteration 461 : 1.2525271013297758
Newton iteration 0
err = 0.012558496465086446
Newton iteration 1
err = 5.350340246053722e-15
cost at iteration 462 : 1.2523634156358718
Newton iteration 0
err = 0.012528187693897547
Newton iteration 1
err = 5.184389037328501e-15
cost at iteration 463 : 1.2522100712998747
Newton iteration 0
err = 0.012498255615436701
Newton iteration 1
err = 5.319012884650349e-15
cost at iteration 464 : 1.2520669941458866
Newton iteration 0
err = 0.01246869759407648
Newton iteration 1
err = 5.560518466604486e-15
cost at iteration 465 : 1.2519341107224065
Newton iteration 0
err = 0.012439511005175144
Newton iteration 1
err = 6.059595639270101e-15
cost at iteration 466 : 1.2518113482945947
Newton iteration 0
err = 0.012410693235454653
Newton iteration 1
err = 5.520391388566282e-15
cost at iteration 467 : 1.251698634836755
Newton iteration 0
err = 0.012382241682422003
Newton iteration 1
err = 5.3865678087956365e-15
cost at iteration 468 : 1.2515958990247538
Newton iteration 0
err = 0.01235415375482515
Newton iteration 1
err = 5.2858908212564884e-15
cost at iteration 469 : 1.2515030702288537
Newton iteration 0
err = 0.01232642687193126
Newton iteration 1
err = 6.327437380811941e-15
cost at iteration 470 : 1.2514200785061478
Newton iteration 0
err = 0.01229905846382115
Newton iteration 1
err = 5.499855409901974e-15
cost at iteration 471 : 1.2513468545934314
Newton iteration 0
err = 0.012272045971458015
Newton iteration 1
err = 5.175984348766833e-15
cost at iteration 472 : 1.25128332990008
Newton iteration 0
err = 0.012245386846256773
Newton iteration 1
err = 6.083251342207983e-15
cost at iteration 473 : 1.2512294365010628
Newton iteration 0
err = 0.012219078550117295
Newton iteration 1
err = 5.80007850161522e-15
cost at iteration 474 : 1.2511851071297697
Newton iteration 0
err = 0.012193118555581094
Newton iteration 1
err = 5.45350643293988e-15
cost at iteration 475 : 1.2511502751713044
Newton iteration 0
err = 0.01216750434524448
Newton iteration 1
err = 5.961584381602905e-15
cost at iteration 476 : 1.2511248746554224
Newton iteration 0
err = 0.012142233412599776
Newton iteration 1
err = 5.273484203131633e-15
cost at iteration 477 : 1.2511088402501607
Newton iteration 0
err = 0.012117303260865664
Newton iteration 1
err = 5.1859008661746714e-15
cost at iteration 478 : 1.2511021072547472
Newton iteration 0
err = 0.012092711403952024
Newton iteration 1
err = 5.212613895612797e-15
scale = 5e-05
New cost during line search: 1.24416530401301
cost at iteration 479 : 1.24416530401301
Newton iteration 0
err = 0.006034500361259126
Newton iteration 1
err = 5.1850160002518524e-15
scale = 5e-05
New cost during line search: 1.2441630501393042
cost at iteration 480 : 1.2441630501393042
Newton iteration 0
err = 0.006028510531511427
Newton iteration 1
err = 5.699703000782396e-15
scale = 5e-05
New cost during line search: 1.2441630851978775
scale = 2.5e-05
New cost during line search: 1.2441629892777155
cost at iteration 481 : 1.2441629892777155
Newton iteration 0
err = 0.003011349663225307
Newton iteration 1
err = 5.3918315175923764e-15
scale = 5e-05
New cost during line search: 1.244164165163216
scale = 2.5e-05
New cost during line search: 1.2441634989728694
scale = 1.25e-05
New cost during line search: 1.2441632245647072
scale = 6.25e-06
New cost during line search: 1.244163102031241
scale = 3.125e-06
New cost during line search: 1.2441630444319927
scale = 1.5625e-06
New cost during line search: 1.2441630165492428
scale = 7.8125e-07
New cost during line search: 1.2441630028370567
scale = 3.90625e-07
New cost during line search: 1.244162996038274
scale = 1.953125e-07
New cost during line search: 1.2441629926532174
scale = 9.765625e-08
New cost during line search: 1.2441629909642733
scale = 4.8828125e-08
New cost during line search: 1.2441629901206601
scale = 2.44140625e-08
New cost during line search: 1.24416298969912
scale = 1.220703125e-08
New cost during line search: 1.2441629894884019
scale = 6.103515625e-09
New cost during line search: 1.2441629893830144
scale = 3.0517578125e-09
New cost during line search: 1.2441629893303696
scale = 1.52587890625e-09
New cost during line search: 1.244162989304043
scale = 7.62939453125e-10
New cost during line search: 1.2441629892908286
scale = 3.814697265625e-10
New cost during line search: 1.2441629892842996
scale = 1.9073486328125e-10
New cost during line search: 1.2441629892809771
scale = 9.5367431640625e-11
New cost during line search: 1.2441629892793518
scale = 4.76837158203125e-11
New cost during line search: 1.2441629892785286
scale = 2.384185791015625e-11
New cost during line search: 1.2441629892780874
scale = 1.1920928955078126e-11
New cost during line search: 1.2441629892778983
scale = 5.960464477539063e-12
New cost during line search: 1.2441629892777857
scale = 2.9802322387695314e-12
New cost during line search: 1.2441629892777466
scale = 1.4901161193847657e-12
New cost during line search: 1.2441629892777295
scale = 7.450580596923828e-13
No more descent can be found
=================Mesh Quality Check=================
max quality: 1.441807871875525 min quality: 0.5167710698008765 mean quality: 0.8594784827643153
===================================================
check_mesh_quality()
c = pose_deformation_equation(p=5)
gfX = GridFunction(VEC)
gfX.Set((x *(1-x)*y*(1-y),0))
gfset = GridFunction(VEC)
gfsetOld = GridFunction(VEC)
gfset.Set((0,0))
scene = Draw(gfset,mesh,"gfset")
SetVisualization (deformation=True)
optimize_shape(linear=False)
mesh.SetDeformation(gfset)
check_mesh_quality()
=================Mesh Quality Check=================
max quality: 1.441807871875525 min quality: 0.5167710698008765 mean quality: 0.8594784827643153
===================================================
Using p = 5
Cost at initial design 13.962388102754325
cost at iteration 0 : 13.962388102754325
Newton iteration 0
err = 69.38374121332043
Newton iteration 1
err = 95525973.66348527
Newton iteration 2
err = 46534454.951519154
Newton iteration 3
err = 22668761.32833453
Newton iteration 4
err = 11042844.290232865
Newton iteration 5
err = 5379403.322999564
Newton iteration 6
err = 2620518.71335221
Newton iteration 7
err = 1276557.6244199935
Newton iteration 8
err = 621861.3742348586
Newton iteration 9
err = 302933.1079457248
Newton iteration 10
err = 147570.61879367885
Newton iteration 11
err = 71887.44606638905
Newton iteration 12
err = 35019.19889199272
Newton iteration 13
err = 17059.226994771438
Newton iteration 14
err = 8310.219017450478
Newton iteration 15
err = 4048.233680696883
Newton iteration 16
err = 1972.0533601981629
Newton iteration 17
err = 960.6643252697774
Newton iteration 18
err = 467.97698908025376
Newton iteration 19
err = 227.9697209864447
Newton iteration 20
err = 111.05296347167638
Newton iteration 21
err = 54.098079001127424
Newton iteration 22
err = 26.35274210207727
Newton iteration 23
err = 12.837171975807319
Newton iteration 24
err = 6.256405109294945
Newton iteration 25
err = 3.055553881474144
Newton iteration 26
err = 1.5032775702643142
Newton iteration 27
err = 0.7539379142719076
Newton iteration 28
err = 0.39516650093499783
Newton iteration 29
err = 0.1962113835173666
Newton iteration 30
err = 0.07616340191361548
Newton iteration 31
err = 0.0161952102187628
Newton iteration 32
err = 0.0008527401283532236
Newton iteration 33
err = 2.4418088889929785e-06
Newton iteration 34
err = 2.0062203907033034e-11
cost at iteration 1 : 13.934752179253156
Newton iteration 0
err = 0.1653834014651807
Newton iteration 1
err = 3.979479133693186e-05
Newton iteration 2
err = 7.957871044945422e-10
cost at iteration 2 : 13.908869188124667
Newton iteration 0
err = 0.16007838567877328
Newton iteration 1
err = 3.870229556976139e-05
Newton iteration 2
err = 7.584259436602396e-10
cost at iteration 3 : 13.883178024567632
Newton iteration 0
err = 0.15494615756371
Newton iteration 1
err = 3.766084712112133e-05
Newton iteration 2
err = 7.234984402115556e-10
cost at iteration 4 : 13.857668858567187
Newton iteration 0
err = 0.1499811640330529
Newton iteration 1
err = 3.666793450350137e-05
Newton iteration 2
err = 6.908339011113341e-10
cost at iteration 5 : 13.832332481890504
Newton iteration 0
err = 0.14517803283706313
Newton iteration 1
err = 3.572121220495211e-05
Newton iteration 2
err = 6.60263167729356e-10
cost at iteration 6 : 13.807160268538782
Newton iteration 0
err = 0.14053156671882644
Newton iteration 1
err = 3.4818486783809284e-05
Newton iteration 2
err = 6.316453715935145e-10
cost at iteration 7 : 13.78214413771727
Newton iteration 0
err = 0.13603673776060105
Newton iteration 1
err = 3.395770417769166e-05
Newton iteration 2
err = 6.048447608968938e-10
cost at iteration 8 : 13.757276519160792
Newton iteration 0
err = 0.13168868191436386
Newton iteration 1
err = 3.3136938089904465e-05
Newton iteration 2
err = 5.797341565124422e-10
cost at iteration 9 : 13.732550320666476
Newton iteration 0
err = 0.12748269371029505
Newton iteration 1
err = 3.235437933179608e-05
Newton iteration 2
err = 5.56202405922325e-10
cost at iteration 10 : 13.707958897692258
Newton iteration 0
err = 0.12341422113652271
Newton iteration 1
err = 3.1608326096894855e-05
Newton iteration 2
err = 5.341329099250298e-10
cost at iteration 11 : 13.683496024889928
Newton iteration 0
err = 0.11947886068575622
Newton iteration 1
err = 3.089717506403407e-05
Newton iteration 2
err = 5.134375212418207e-10
cost at iteration 12 : 13.659155869449927
Newton iteration 0
err = 0.11567235256220998
Newton iteration 1
err = 3.0219413283596676e-05
Newton iteration 2
err = 4.940149573669752e-10
cost at iteration 13 : 13.634932966142163
Newton iteration 0
err = 0.11199057604318337
Newton iteration 1
err = 2.9573610736510203e-05
Newton iteration 2
err = 4.757887483157533e-10
cost at iteration 14 : 13.61082219394503
Newton iteration 0
err = 0.1084295449911515
Newton iteration 1
err = 2.895841360895716e-05
Newton iteration 2
err = 4.58674688496704e-10
cost at iteration 15 : 13.586818754161266
Newton iteration 0
err = 0.10498540351029195
Newton iteration 1
err = 2.8372538117719427e-05
Newton iteration 2
err = 4.426029538964038e-10
cost at iteration 16 : 13.562918149926984
Newton iteration 0
err = 0.10165442174235445
Newton iteration 1
err = 2.781476491709504e-05
Newton iteration 2
err = 4.2750353825139356e-10
cost at iteration 17 : 13.539116167023822
Newton iteration 0
err = 0.09843299179885555
Newton iteration 1
err = 2.7283934032305238e-05
Newton iteration 2
err = 4.1331902806538765e-10
cost at iteration 18 : 13.515408855912371
Newton iteration 0
err = 0.09531762382206924
Newton iteration 1
err = 2.6778940249599286e-05
Newton iteration 2
err = 3.9998511048991195e-10
cost at iteration 19 : 13.491792514908903
Newton iteration 0
err = 0.09230494217315402
Newton iteration 1
err = 2.629872896108289e-05
Newton iteration 2
err = 3.8745411602311227e-10
cost at iteration 20 : 13.468263674432137
Newton iteration 0
err = 0.08939168174115018
Newton iteration 1
err = 2.5842292424781e-05
Newton iteration 2
err = 3.756707515319509e-10
cost at iteration 21 : 13.444819082252364
Newton iteration 0
err = 0.0865746843694662
Newton iteration 1
err = 2.540866636047218e-05
Newton iteration 2
err = 3.6459457775500474e-10
cost at iteration 22 : 13.42145568967985
Newton iteration 0
err = 0.08385089539543973
Newton iteration 1
err = 2.499692694405459e-05
Newton iteration 2
err = 3.5418250913762986e-10
cost at iteration 23 : 13.398170638630528
Newton iteration 0
err = 0.08121736029902628
Newton iteration 1
err = 2.4606188084728796e-05
Newton iteration 2
err = 3.443927347695955e-10
cost at iteration 24 : 13.374961249516515
Newton iteration 0
err = 0.07867122145659185
Newton iteration 1
err = 2.423559898811694e-05
Newton iteration 2
err = 3.35190857824499e-10
cost at iteration 25 : 13.351825009905221
Newton iteration 0
err = 0.07620971499641337
Newton iteration 1
err = 2.3884342009685574e-05
Newton iteration 2
err = 3.2654130299985124e-10
cost at iteration 26 : 13.328759563901423
Newton iteration 0
err = 0.07383016775139907
Newton iteration 1
err = 2.355163072972866e-05
Newton iteration 2
err = 3.184129010365336e-10
cost at iteration 27 : 13.305762702204017
Newton iteration 0
err = 0.07152999430656994
Newton iteration 1
err = 2.323670824209527e-05
Newton iteration 2
err = 3.107789920303499e-10
cost at iteration 28 : 13.282832352795541
Newton iteration 0
err = 0.06930669413638757
Newton iteration 1
err = 2.2938845657968435e-05
Newton iteration 2
err = 3.0361176175778246e-10
cost at iteration 29 : 13.259966572224714
Newton iteration 0
err = 0.06715784883000046
Newton iteration 1
err = 2.265734076610297e-05
Newton iteration 2
err = 2.9688810519160065e-10
cost at iteration 30 : 13.237163537442639
Newton iteration 0
err = 0.06508111939966152
Newton iteration 1
err = 2.2391516862928474e-05
Newton iteration 2
err = 2.9058246355070095e-10
cost at iteration 31 : 13.214421538158547
Newton iteration 0
err = 0.06307424367048843
Newton iteration 1
err = 2.2140721702169e-05
Newton iteration 2
err = 2.846745086947417e-10
cost at iteration 32 : 13.191738969682694
Newton iteration 0
err = 0.06113503374691637
Newton iteration 1
err = 2.1904326588968548e-05
Newton iteration 2
err = 2.7914525544891714e-10
cost at iteration 33 : 13.169114326223509
Newton iteration 0
err = 0.059261373554205375
Newton iteration 1
err = 2.168172555486305e-05
Newton iteration 2
err = 2.7397597230784513e-10
cost at iteration 34 : 13.146546194611872
Newton iteration 0
err = 0.057451216450950256
Newton iteration 1
err = 2.147233463200448e-05
Newton iteration 2
err = 2.691545187931426e-10
cost at iteration 35 : 13.124033248424368
Newton iteration 0
err = 0.05570258291042372
Newton iteration 1
err = 2.1275591218018307e-05
Newton iteration 2
err = 2.646617110668048e-10
cost at iteration 36 : 13.10157424248032
Newton iteration 0
err = 0.05401355826715598
Newton iteration 1
err = 2.109095345805131e-05
Newton iteration 2
err = 2.6048540761035893e-10
cost at iteration 37 : 13.079168007689233
Newton iteration 0
err = 0.05238229052713169
Newton iteration 1
err = 2.091789974569797e-05
Newton iteration 2
err = 2.5660795548213085e-10
cost at iteration 38 : 13.056813446225696
Newton iteration 0
err = 0.05080698823785021
Newton iteration 1
err = 2.075592820393548e-05
Newton iteration 2
err = 2.530263506901833e-10
cost at iteration 39 : 13.034509527011236
Newton iteration 0
err = 0.04928591841650187
Newton iteration 1
err = 2.060455623027167e-05
Newton iteration 2
err = 2.497234385687985e-10
cost at iteration 40 : 13.012255281484032
Newton iteration 0
err = 0.04781740453345619
Newton iteration 1
err = 2.0463320062343084e-05
Newton iteration 2
err = 2.466916017938375e-10
cost at iteration 41 : 12.990049799637204
Newton iteration 0
err = 0.04639982454906527
Newton iteration 1
err = 2.033177436745505e-05
Newton iteration 2
err = 2.439211066824284e-10
cost at iteration 42 : 12.967892226309438
Newton iteration 0
err = 0.045031609000942496
Newton iteration 1
err = 2.0209491822957732e-05
Newton iteration 2
err = 2.414044578616109e-10
cost at iteration 43 : 12.94578175771133
Newton iteration 0
err = 0.043711239140922356
Newton iteration 1
err = 2.0096062729260756e-05
Newton iteration 2
err = 2.3913171029803674e-10
cost at iteration 44 : 12.923717638173049
Newton iteration 0
err = 0.042437245117982594
Newton iteration 1
err = 1.9991094615979194e-05
Newton iteration 2
err = 2.370986581939119e-10
cost at iteration 45 : 12.90169915709876
Newton iteration 0
err = 0.04120820420733661
Newton iteration 1
err = 1.9894211849599213e-05
Newton iteration 2
err = 2.352992442704463e-10
cost at iteration 46 : 12.879725646114599
Newton iteration 0
err = 0.040022739082168385
Newton iteration 1
err = 1.980505523952205e-05
Newton iteration 2
err = 2.337256606303884e-10
cost at iteration 47 : 12.857796476398315
Newton iteration 0
err = 0.038879516128408585
Newton iteration 1
err = 1.972328164685503e-05
Newton iteration 2
err = 2.3237397221252913e-10
cost at iteration 48 : 12.835911056178656
Newton iteration 0
err = 0.03777724379966409
Newton iteration 1
err = 1.9648563580626895e-05
Newton iteration 2
err = 2.312400645822846e-10
cost at iteration 49 : 12.814068828393879
Newton iteration 0
err = 0.03671467101220227
Newton iteration 1
err = 1.9580588802687717e-05
Newton iteration 2
err = 2.3031968943383468e-10
cost at iteration 50 : 12.792269268499494
Newton iteration 0
err = 0.035690585579110964
Newton iteration 1
err = 1.951905991753995e-05
Newton iteration 2
err = 2.2960584151286877e-10
cost at iteration 51 : 12.770511882414597
Newton iteration 0
err = 0.034703812682531816
Newton iteration 1
err = 1.946369395871884e-05
Newton iteration 2
err = 2.2909982075148822e-10
cost at iteration 52 : 12.748796204600286
Newton iteration 0
err = 0.03375321338396836
Newton iteration 1
err = 1.9414221986048657e-05
Newton iteration 2
err = 2.2879579784684218e-10
cost at iteration 53 : 12.727121796258643
Newton iteration 0
err = 0.03283768317226119
Newton iteration 1
err = 1.9370388659518247e-05
Newton iteration 2
err = 2.2869239362203937e-10
cost at iteration 54 : 12.705488243647784
Newton iteration 0
err = 0.03195615054965868
Newton iteration 1
err = 1.933195183103412e-05
Newton iteration 2
err = 2.2878511469500165e-10
cost at iteration 55 : 12.68389515650325
Newton iteration 0
err = 0.031107575655835402
Newton iteration 1
err = 1.929868212607858e-05
Newton iteration 2
err = 2.2907324036470654e-10
cost at iteration 56 : 12.662342166560158
Newton iteration 0
err = 0.030290948931083875
Newton iteration 1
err = 1.9270362523484163e-05
Newton iteration 2
err = 2.2955564201463845e-10
cost at iteration 57 : 12.640828926169023
Newton iteration 0
err = 0.029505289818631777
Newton iteration 1
err = 1.9246787951226814e-05
Newton iteration 2
err = 2.3023049095460936e-10
cost at iteration 58 : 12.61935510700028
Newton iteration 0
err = 0.028749645508260254
Newton iteration 1
err = 1.922776486537214e-05
Newton iteration 2
err = 2.3109646983083572e-10
cost at iteration 59 : 12.597920398830784
Newton iteration 0
err = 0.028023089721468292
Newton iteration 1
err = 1.9213110858595057e-05
Newton iteration 2
err = 2.321511859820282e-10
cost at iteration 60 : 12.576524508407468
Newton iteration 0
err = 0.02732472154018269
Newton iteration 1
err = 1.9202654246036426e-05
Newton iteration 2
err = 2.3339566618893463e-10
cost at iteration 61 : 12.555167158383666
Newton iteration 0
err = 0.02665366428086308
Newton iteration 1
err = 1.9196233687256903e-05
Newton iteration 2
err = 2.348285220768461e-10
cost at iteration 62 : 12.533848086322804
Newton iteration 0
err = 0.026009064415096463
Newton iteration 1
err = 1.9193697788336936e-05
Newton iteration 2
err = 2.364489416968931e-10
cost at iteration 63 : 12.51256704376589
Newton iteration 0
err = 0.025390090539169886
Newton iteration 1
err = 1.9194904734163153e-05
Newton iteration 2
err = 2.3825846491152083e-10
cost at iteration 64 : 12.491323795357815
Newton iteration 0
err = 0.024795932394166
Newton iteration 1
err = 1.9199721918646004e-05
Newton iteration 2
err = 2.4025729305966295e-10
cost at iteration 65 : 12.470118118029832
Newton iteration 0
err = 0.02422579993877051
Newton iteration 1
err = 1.9208025591171783e-05
Newton iteration 2
err = 2.424435362353477e-10
cost at iteration 66 : 12.44894980023375
Newton iteration 0
err = 0.0236789224761036
Newton iteration 1
err = 1.9219700503408023e-05
Newton iteration 2
err = 2.448217031206039e-10
cost at iteration 67 : 12.42781864122484
Newton iteration 0
err = 0.02315454783665432
Newton iteration 1
err = 1.923463958612323e-05
Newton iteration 2
err = 2.4738942506964383e-10
cost at iteration 68 : 12.40672445039091
Newton iteration 0
err = 0.022651941618807476
Newton iteration 1
err = 1.925274361801653e-05
Newton iteration 2
err = 2.501497563538057e-10
cost at iteration 69 : 12.385667046623333
Newton iteration 0
err = 0.022170386487614553
Newton iteration 1
err = 1.9273920921481275e-05
Newton iteration 2
err = 2.5310335673416403e-10
cost at iteration 70 : 12.364646257728754
Newton iteration 0
err = 0.021709181533457703
Newton iteration 1
err = 1.929808705961451e-05
Newton iteration 2
err = 2.5625381777137316e-10
cost at iteration 71 : 12.343661919878256
Newton iteration 0
err = 0.0212676416899205
Newton iteration 1
err = 1.9325164555610053e-05
Newton iteration 2
err = 2.5959944349798417e-10
cost at iteration 72 : 12.32271387709122
Newton iteration 0
err = 0.020845097212369715
Newton iteration 1
err = 1.935508262371922e-05
Newton iteration 2
err = 2.631453977900381e-10
cost at iteration 73 : 12.301801980752408
Newton iteration 0
err = 0.020440893215118897
Newton iteration 1
err = 1.9387776896846784e-05
Newton iteration 2
err = 2.6689592041981143e-10
cost at iteration 74 : 12.280926089159916
Newton iteration 0
err = 0.020054389267094404
Newton iteration 1
err = 1.9423189205114426e-05
Newton iteration 2
err = 2.708494746066684e-10
cost at iteration 75 : 12.260086067101575
Newton iteration 0
err = 0.019684959044516368
Newton iteration 1
err = 1.9461267314617638e-05
Newton iteration 2
err = 2.7501388383764224e-10
cost at iteration 76 : 12.23928178545892
Newton iteration 0
err = 0.019331990037441543
Newton iteration 1
err = 1.9501964732748676e-05
Newton iteration 2
err = 2.793915292194201e-10
cost at iteration 77 : 12.21851312083586
Newton iteration 0
err = 0.018994883308412872
Newton iteration 1
err = 1.9545240486624467e-05
Newton iteration 2
err = 2.839843674446678e-10
cost at iteration 78 : 12.197779955211471
Newton iteration 0
err = 0.018673053299673954
Newton iteration 1
err = 1.9591058936352454e-05
Newton iteration 2
err = 2.887981910456866e-10
cost at iteration 79 : 12.17708217561512
Newton iteration 0
err = 0.018365927684582145
Newton iteration 1
err = 1.9639389591711877e-05
Newton iteration 2
err = 2.938406026635228e-10
cost at iteration 80 : 12.15641967382163
Newton iteration 0
err = 0.018072947259868247
Newton iteration 1
err = 1.9690206947247617e-05
Newton iteration 2
err = 2.991119047573253e-10
cost at iteration 81 : 12.135792346066657
Newton iteration 0
err = 0.017793565873736518
Newton iteration 1
err = 1.9743490326403185e-05
Newton iteration 2
err = 3.0462088217597053e-10
cost at iteration 82 : 12.115200092780002
Newton iteration 0
err = 0.017527250384397196
Newton iteration 1
err = 1.9799223733302733e-05
Newton iteration 2
err = 3.1037277445294874e-10
cost at iteration 83 : 12.094642818336116
Newton iteration 0
err = 0.01727348064523154
Newton iteration 1
err = 1.9857395733847835e-05
Newton iteration 2
err = 3.1637348724595416e-10
cost at iteration 84 : 12.074120430820338
Newton iteration 0
err = 0.01703174950953975
Newton iteration 1
err = 1.9917999327007666e-05
Newton iteration 2
err = 3.226325483948901e-10
cost at iteration 85 : 12.053632841810622
Newton iteration 0
err = 0.016801562851512346
Newton iteration 1
err = 1.9981031848223014e-05
Newton iteration 2
err = 3.2915526433688025e-10
cost at iteration 86 : 12.033179966172652
Newton iteration 0
err = 0.016582439596270215
Newton iteration 1
err = 2.004649486746529e-05
Newton iteration 2
err = 3.359501607842966e-10
cost at iteration 87 : 12.012761721868623
Newton iteration 0
err = 0.016373911755880523
Newton iteration 1
err = 2.011439411338072e-05
Newton iteration 2
err = 3.430267176814542e-10
cost at iteration 88 : 11.9923780297781
Newton iteration 0
err = 0.016175524464836742
Newton iteration 1
err = 2.0184739404896924e-05
Newton iteration 2
err = 3.5039366913150407e-10
cost at iteration 89 : 11.972028813530821
Newton iteration 0
err = 0.015986836011493054
Newton iteration 1
err = 2.0257544587512306e-05
Newton iteration 2
err = 3.58061673042036e-10
cost at iteration 90 : 11.951713999349382
Newton iteration 0
err = 0.015807417861197735
Newton iteration 1
err = 2.0332827494493125e-05
Newton iteration 2
err = 3.6604139661135556e-10
cost at iteration 91 : 11.931433515903233
Newton iteration 0
err = 0.015636854667044085
Newton iteration 1
err = 2.0410609908154145e-05
Newton iteration 2
err = 3.7434359164310153e-10
cost at iteration 92 : 11.911187294171327
Newton iteration 0
err = 0.015474744265819422
Newton iteration 1
err = 2.0490917542715592e-05
Newton iteration 2
err = 3.829805289931906e-10
cost at iteration 93 : 11.890975267314026
Newton iteration 0
err = 0.015320697655525444
Newton iteration 1
err = 2.0573780029864628e-05
Newton iteration 2
err = 3.919654526177815e-10
cost at iteration 94 : 11.870797370553184
Newton iteration 0
err = 0.015174338953247718
Newton iteration 1
err = 2.0659230928063215e-05
Newton iteration 2
err = 4.013134672759139e-10
cost at iteration 95 : 11.850653541060415
Newton iteration 0
err = 0.015035305331077577
Newton iteration 1
err = 2.074730772780248e-05
Newton iteration 2
err = 4.110379676218823e-10
cost at iteration 96 : 11.830543717852102
Newton iteration 0
err = 0.01490324692931117
Newton iteration 1
err = 2.0838051892650166e-05
Newton iteration 2
err = 4.21154220623152e-10
cost at iteration 97 : 11.81046784169165
Newton iteration 0
err = 0.014777826746352021
Newton iteration 1
err = 2.0931508883741437e-05
Newton iteration 2
err = 4.316813778380524e-10
cost at iteration 98 : 11.790425854998444
Newton iteration 0
err = 0.014658720504719235
Newton iteration 1
err = 2.1027728228847623e-05
Newton iteration 2
err = 4.4263528697141394e-10
cost at iteration 99 : 11.770417701762034
Newton iteration 0
err = 0.014545616494004859
Newton iteration 1
err = 2.1126763575502378e-05
Newton iteration 2
err = 4.5403539840906605e-10
cost at iteration 100 : 11.750443327463028
Newton iteration 0
err = 0.01443821539114155
Newton iteration 1
err = 2.122867277838644e-05
Newton iteration 2
err = 4.659037570571297e-10
cost at iteration 101 : 11.730502678998207
Newton iteration 0
err = 0.014336230058629842
Newton iteration 1
err = 2.133351799243295e-05
Newton iteration 2
err = 4.782611136332689e-10
cost at iteration 102 : 11.710595704611597
Newton iteration 0
err = 0.01423938532270385
Newton iteration 1
err = 2.144136577921816e-05
Newton iteration 2
err = 4.911307903137802e-10
cost at iteration 103 : 11.69072235382991
Newton iteration 0
err = 0.014147417732197971
Newton iteration 1
err = 2.1552287227326138e-05
Newton iteration 2
err = 5.045399853199991e-10
cost at iteration 104 : 11.670882577401686
Newton iteration 0
err = 0.014060075300214232
Newton iteration 1
err = 2.1666358098721194e-05
Newton iteration 2
err = 5.185109486099702e-10
cost at iteration 105 : 11.651076327241729
Newton iteration 0
err = 0.013977117230139353
Newton iteration 1
err = 2.1783658972485023e-05
Newton iteration 2
err = 5.330748779614502e-10
cost at iteration 106 : 11.631303556378807
Newton iteration 0
err = 0.013898313628141494
Newton iteration 1
err = 2.190427541819263e-05
Newton iteration 2
err = 5.48261561764901e-10
cost at iteration 107 : 11.611564218906677
Newton iteration 0
err = 0.013823445203732573
Newton iteration 1
err = 2.202829818198757e-05
Newton iteration 2
err = 5.64104563503875e-10
cost at iteration 108 : 11.59185826993931
Newton iteration 0
err = 0.013752302960664497
Newton iteration 1
err = 2.2155823386193663e-05
Newton iteration 2
err = 5.806367617766908e-10
cost at iteration 109 : 11.572185665568735
Newton iteration 0
err = 0.01368468787998681
Newton iteration 1
err = 2.2286952753832947e-05
Newton iteration 2
err = 5.9789720255865e-10
cost at iteration 110 : 11.552546362825959
Newton iteration 0
err = 0.01362041059691616
Newton iteration 1
err = 2.2421793839785554e-05
Newton iteration 2
err = 6.159236576046548e-10
cost at iteration 111 : 11.532940319645318
Newton iteration 0
err = 0.01355929107373845
Newton iteration 1
err = 2.256046029099803e-05
Newton iteration 2
err = 6.347601047436843e-10
cost at iteration 112 : 11.513367494830643
Newton iteration 0
err = 0.013501158270124457
Newton iteration 1
err = 2.2703072118057304e-05
Newton iteration 2
err = 6.544522667119729e-10
cost at iteration 113 : 11.49382784802461
Newton iteration 0
err = 0.013445849812561434
Newton iteration 1
err = 2.2849755989266716e-05
Newton iteration 2
err = 6.75048097427665e-10
cost at iteration 114 : 11.474321339680246
Newton iteration 0
err = 0.013393211664539233
Newton iteration 1
err = 2.300064553798688e-05
Newton iteration 2
err = 6.965990842230796e-10
cost at iteration 115 : 11.45484793103473
Newton iteration 0
err = 0.013343097798627279
Newton iteration 1
err = 2.3155881699313018e-05
Newton iteration 2
err = 7.191631784842471e-10
cost at iteration 116 : 11.43540758408521
Newton iteration 0
err = 0.013295369871895071
Newton iteration 1
err = 2.3315613053601007e-05
Newton iteration 2
err = 7.427980683812784e-10
cost at iteration 117 : 11.416000261566968
Newton iteration 0
err = 0.01324989690581398
Newton iteration 1
err = 2.3479996194372114e-05
Newton iteration 2
err = 7.67570677313356e-10
cost at iteration 118 : 11.396625926932856
Newton iteration 0
err = 0.013206554971345267
Newton iteration 1
err = 2.36491961197143e-05
Newton iteration 2
err = 7.935477015028336e-10
cost at iteration 119 : 11.377284544335392
Newton iteration 0
err = 0.013165226880569389
Newton iteration 1
err = 2.3823386628648956e-05
Newton iteration 2
err = 8.20804991172258e-10
cost at iteration 120 : 11.35797607860985
Newton iteration 0
err = 0.013125801885137164
Newton iteration 1
err = 2.4002750746670334e-05
Newton iteration 2
err = 8.494206907645162e-10
cost at iteration 121 : 11.338700495259182
Newton iteration 0
err = 0.013088175382484624
Newton iteration 1
err = 2.4187481157376705e-05
Newton iteration 2
err = 8.794786734006334e-10
cost at iteration 122 : 11.319457760441008
Newton iteration 0
err = 0.013052248630370458
Newton iteration 1
err = 2.4377780654375423e-05
Newton iteration 2
err = 9.110734558766703e-10
cost at iteration 123 : 11.300247840955231
Newton iteration 0
err = 0.01301792846980446
Newton iteration 1
err = 2.457386260109387e-05
Newton iteration 2
err = 9.443018561769567e-10
cost at iteration 124 : 11.281070704233523
Newton iteration 0
err = 0.012985127057254117
Newton iteration 1
err = 2.4775951402067472e-05
Newton iteration 2
err = 9.7926818329193e-10
cost at iteration 125 : 11.261926318330408
Newton iteration 0
err = 0.012953761605895824
Newton iteration 1
err = 2.498428297817992e-05
Newton iteration 2
err = 1.0160863631884724e-09
cost at iteration 126 : 11.24281465191515
Newton iteration 0
err = 0.012923754136357716
Newton iteration 1
err = 2.519910524733234e-05
Newton iteration 2
err = 1.0548784742779523e-09
cost at iteration 127 : 11.223735674265345
Newton iteration 0
err = 0.01289503123703417
Newton iteration 1
err = 2.5420678605761696e-05
Newton iteration 2
err = 1.0957739518233584e-09
cost at iteration 128 : 11.20468935526117
Newton iteration 0
err = 0.012867523833871042
Newton iteration 1
err = 2.564927639418488e-05
Newton iteration 2
err = 1.1389127165181045e-09
cost at iteration 129 : 11.185675665381524
Newton iteration 0
err = 0.012841166969829003
Newton iteration 1
err = 2.5885185364067842e-05
Newton iteration 2
err = 1.1844455041994396e-09
cost at iteration 130 : 11.1666945757007
Newton iteration 0
err = 0.012815899593699812
Newton iteration 1
err = 2.612870611265229e-05
Newton iteration 2
err = 1.232533474787202e-09
cost at iteration 131 : 11.147746057886511
Newton iteration 0
err = 0.012791664358429943
Newton iteration 1
err = 2.6380153505690426e-05
Newton iteration 2
err = 1.2833479819657505e-09
cost at iteration 132 : 11.12883008419928
Newton iteration 0
err = 0.012768407428599802
Newton iteration 1
err = 2.6639857055408187e-05
Newton iteration 2
err = 1.3370762977409443e-09
cost at iteration 133 : 11.10994662749188
Newton iteration 0
err = 0.012746078296952785
Newton iteration 1
err = 2.690816125670828e-05
Newton iteration 2
err = 1.393915233529478e-09
cost at iteration 134 : 11.091095661211272
Newton iteration 0
err = 0.012724629609825467
Newton iteration 1
err = 2.7185425871996256e-05
Newton iteration 2
err = 1.4540795181166342e-09
cost at iteration 135 : 11.072277159400391
Newton iteration 0
err = 0.012704017001136833
Newton iteration 1
err = 2.747202615075324e-05
Newton iteration 2
err = 1.5177930923282437e-09
cost at iteration 136 : 11.053491096701482
Newton iteration 0
err = 0.012684198934655378
Newton iteration 1
err = 2.7768352966197956e-05
Newton iteration 2
err = 1.5853022893425856e-09
cost at iteration 137 : 11.034737448360112
Newton iteration 0
err = 0.012665136554533088
Newton iteration 1
err = 2.8074812864442036e-05
Newton iteration 2
err = 1.6568643119494645e-09
cost at iteration 138 : 11.016016190230676
Newton iteration 0
err = 0.012646793543384932
Newton iteration 1
err = 2.8391828007110677e-05
Newton iteration 2
err = 1.732756715136318e-09
cost at iteration 139 : 10.997327298782277
Newton iteration 0
err = 0.0126291359880176
Newton iteration 1
err = 2.871983599469449e-05
Newton iteration 2
err = 1.8132740508754525e-09
cost at iteration 140 : 10.978670751106096
Newton iteration 0
err = 0.012612132252244815
Newton iteration 1
err = 2.9059289556744532e-05
Newton iteration 2
err = 1.898728731499999e-09
cost at iteration 141 : 10.960046524923301
Newton iteration 0
err = 0.012595752856544335
Newton iteration 1
err = 2.9410656097024283e-05
Newton iteration 2
err = 1.989451688034235e-09
cost at iteration 142 : 10.94145459859431
Newton iteration 0
err = 0.012579970364208846
Newton iteration 1
err = 2.9774417074828568e-05
Newton iteration 2
err = 2.0857948895311195e-09
cost at iteration 143 : 10.922894951128796
Newton iteration 0
err = 0.012564759273698973
Newton iteration 1
err = 3.015106722799608e-05
Newton iteration 2
err = 2.1881262140167715e-09
cost at iteration 144 : 10.904367562196974
Newton iteration 0
err = 0.012550095916704536
Newton iteration 1
err = 3.0541113628607755e-05
Newton iteration 2
err = 2.2968364254121554e-09
cost at iteration 145 : 10.885872412140944
Newton iteration 0
err = 0.01253595836177467
Newton iteration 1
err = 3.094507457994192e-05
Newton iteration 2
err = 2.412329711945419e-09
cost at iteration 146 : 10.867409481988647
Newton iteration 0
err = 0.012522326322899202
Newton iteration 1
err = 3.13634783654517e-05
Newton iteration 2
err = 2.535032036065279e-09
cost at iteration 147 : 10.848978753466945
Newton iteration 0
err = 0.012509181072888996
Newton iteration 1
err = 3.179686190408378e-05
Newton iteration 2
err = 2.665383852604303e-09
cost at iteration 148 : 10.830580209017246
Newton iteration 0
err = 0.012496505361047754
Newton iteration 1
err = 3.224576936246809e-05
Newton iteration 2
err = 2.803838911677908e-09
cost at iteration 149 : 10.812213831810764
Newton iteration 0
err = 0.012484283334820089
Newton iteration 1
err = 3.271075080353895e-05
Newton iteration 2
err = 2.9508666115106892e-09
cost at iteration 150 : 10.79387960576576
Newton iteration 0
err = 0.01247250046494682
Newton iteration 1
err = 3.3192361037198554e-05
Newton iteration 2
err = 3.1069454092879777e-09
cost at iteration 151 : 10.7755775155651
Newton iteration 0
err = 0.012461143473996093
Newton iteration 1
err = 3.369115883843271e-05
Newton iteration 2
err = 3.2725598046699425e-09
cost at iteration 152 : 10.757307546674973
Newton iteration 0
err = 0.012450200267533315
Newton iteration 1
err = 3.420770680369968e-05
Newton iteration 2
err = 3.4482065726890676e-09
cost at iteration 153 : 10.73906968536392
Newton iteration 0
err = 0.0124396598680036
Newton iteration 1
err = 3.47425722001425e-05
Newton iteration 2
err = 3.634385129403419e-09
cost at iteration 154 : 10.720863918723941
Newton iteration 0
err = 0.012429512350786064
Newton iteration 1
err = 3.529632925498369e-05
Newton iteration 2
err = 3.831602352332445e-09
cost at iteration 155 : 10.702690234690634
Newton iteration 0
err = 0.012419748782284053
Newton iteration 1
err = 3.5869563490169264e-05
Newton iteration 2
err = 4.04037551724632e-09
cost at iteration 156 : 10.68454862206549
Newton iteration 0
err = 0.012410361159985889
Newton iteration 1
err = 3.646287880226096e-05
Newton iteration 2
err = 4.261241050538261e-09
cost at iteration 157 : 10.66643907053829
Newton iteration 0
err = 0.012401342354375485
Newton iteration 1
err = 3.707690820468913e-05
Newton iteration 2
err = 4.494764212998257e-09
cost at iteration 158 : 10.648361570709627
Newton iteration 0
err = 0.012392686053052757
Newton iteration 1
err = 3.771232923082586e-05
Newton iteration 2
err = 4.741566314469282e-09
cost at iteration 159 : 10.630316114114846
Newton iteration 0
err = 0.012384386707272339
Newton iteration 1
err = 3.8369885186867116e-05
Newton iteration 2
err = 5.0023479379468866e-09
cost at iteration 160 : 10.612302693247567
Newton iteration 0
err = 0.012376439481937729
Newton iteration 1
err = 3.905041344994022e-05
Newton iteration 2
err = 5.2779387052899975e-09
cost at iteration 161 : 10.594321301583989
Newton iteration 0
err = 0.012368840209780777
Newton iteration 1
err = 3.975488198813815e-05
Newton iteration 2
err = 5.569352490039284e-09
cost at iteration 162 : 10.576371933607524
Newton iteration 0
err = 0.01236158535155964
Newton iteration 1
err = 4.048443505122463e-05
Newton iteration 2
err = 5.877858953409508e-09
cost at iteration 163 : 10.55845458483335
Newton iteration 0
err = 0.012354671964298906
Newton iteration 1
err = 4.124044847075189e-05
Newton iteration 2
err = 6.205074367692843e-09
cost at iteration 164 : 10.54056925183342
Newton iteration 0
err = 0.012348097680111816
Newton iteration 1
err = 4.2024594228211534e-05
Newton iteration 2
err = 6.553066371884739e-09
cost at iteration 165 : 10.522715932262177
Newton iteration 0
err = 0.01234186069892655
Newton iteration 1
err = 4.283891270554718e-05
Newton iteration 2
err = 6.9244586919814264e-09
cost at iteration 166 : 10.504894624882294
Newton iteration 0
err = 0.012335959798989264
Newton iteration 1
err = 4.3685889375221465e-05
Newton iteration 2
err = 7.322552687718486e-09
cost at iteration 167 : 10.48710532959084
Newton iteration 0
err = 0.012330394369104367
Newton iteration 1
err = 4.456853075168658e-05
Newton iteration 2
err = 7.751426945462715e-09
cost at iteration 168 : 10.469348047448188
Newton iteration 0
err = 0.012325164467236461
Newton iteration 1
err = 4.549043224489716e-05
Newton iteration 2
err = 8.216033633436336e-09
cost at iteration 169 : 10.451622780706623
Newton iteration 0
err = 0.012320270909423678
Newton iteration 1
err = 4.64558287751062e-05
Newton iteration 2
err = 8.72225710114321e-09
cost at iteration 170 : 10.433929532842695
Newton iteration 0
err = 0.012315715392350192
Newton iteration 1
err = 4.746961799161584e-05
Newton iteration 2
err = 9.276960180847907e-09
cost at iteration 171 : 10.416268308592462
Newton iteration 0
err = 0.012311500651797438
Newton iteration 1
err = 4.853734658012786e-05
Newton iteration 2
err = 9.887989606364224e-09
cost at iteration 172 : 10.398639113991496
Newton iteration 0
err = 0.012307630656274677
Newton iteration 1
err = 4.9665153013417095e-05
Newton iteration 2
err = 1.0564187102901804e-08
Newton iteration 3
err = 1.8233564040988024e-14
cost at iteration 173 : 10.381041956420185
Newton iteration 0
err = 0.012304110833451495
Newton iteration 1
err = 5.085966562464336e-05
Newton iteration 2
err = 1.1315390989351814e-08
Newton iteration 3
err = 1.7787299707904648e-14
cost at iteration 174 : 10.363476844657052
Newton iteration 0
err = 0.01230094832267872
Newton iteration 1
err = 5.212786276184161e-05
Newton iteration 2
err = 1.2152475052473384e-08
Newton iteration 3
err = 1.7349121556607614e-14
cost at iteration 175 : 10.345943788940067
Newton iteration 0
err = 0.012298152244176751
Newton iteration 1
err = 5.3476910883903163e-05
Newton iteration 2
err = 1.3087412704283389e-08
Newton iteration 3
err = 1.871455396751302e-14
cost at iteration 176 : 10.328442801039486
Newton iteration 0
err = 0.012295733972304414
Newton iteration 1
err = 5.4914004689581795e-05
Newton iteration 2
err = 1.4133413158873195e-08
Newton iteration 3
err = 1.8808958329644537e-14
cost at iteration 177 : 10.310973894341576
Newton iteration 0
err = 0.01229370739825941
Newton iteration 1
err = 5.644623835126511e-05
Newton iteration 2
err = 1.5305099215434312e-08
Newton iteration 3
err = 1.9326667953338592e-14
cost at iteration 178 : 10.293537083947292
Newton iteration 0
err = 0.012292089167024359
Newton iteration 1
err = 5.808053629417801e-05
Newton iteration 2
err = 1.661875188392611e-08
Newton iteration 3
err = 1.8109226494483065e-14
cost at iteration 179 : 10.276132386785072
Newton iteration 0
err = 0.012290898874954062
Newton iteration 1
err = 5.9823664946027574e-05
Newton iteration 2
err = 1.8092613091068748e-08
Newton iteration 3
err = 1.7564372682958154e-14
cost at iteration 180 : 10.258759821739352
Newton iteration 0
err = 0.012290159216890376
Newton iteration 1
err = 6.16823343295594e-05
Newton iteration 2
err = 1.9747261659996388e-08
Newton iteration 3
err = 1.8765461197360116e-14
cost at iteration 181 : 10.241419409796512
Newton iteration 0
err = 0.012289896076918064
Newton iteration 1
err = 6.366338323003658e-05
Newton iteration 2
err = 2.16060791505228e-08
Newton iteration 3
err = 1.7787466749043256e-14
cost at iteration 182 : 10.224111174207465
Newton iteration 0
err = 0.012290138561118276
Newton iteration 1
err = 6.577402774877794e-05
Newton iteration 2
err = 2.3695820893734014e-08
Newton iteration 3
err = 1.801377840474905e-14
cost at iteration 183 : 10.206835140667062
Newton iteration 0
err = 0.012290918976429973
Newton iteration 1
err = 6.802214420423649e-05
Newton iteration 2
err = 2.6047313122456472e-08
Newton iteration 3
err = 1.9448335228159856e-14
cost at iteration 184 : 10.189591337511139
Newton iteration 0
err = 0.012292272762426974
Newton iteration 1
err = 7.0416555552149e-05
Newton iteration 2
err = 2.8696285519762736e-08
Newton iteration 3
err = 1.884096207872271e-14
cost at iteration 185 : 10.172379795929428
Newton iteration 0
err = 0.01229423838576332
Newton iteration 1
err = 7.296729618248695e-05
Newton iteration 2
err = 3.168433802872885e-08
Newton iteration 3
err = 2.078065411271825e-14
cost at iteration 186 : 10.155200550194584
Newton iteration 0
err = 0.012296857206278145
Newton iteration 1
err = 7.56858413437247e-05
Newton iteration 2
err = 3.50600525264954e-08
Newton iteration 3
err = 2.0564839526370023e-14
cost at iteration 187 : 10.138053637906946
Newton iteration 0
err = 0.01230017332275999
Newton iteration 1
err = 7.858530212621059e-05
Newton iteration 2
err = 3.888028032501466e-08
Newton iteration 3
err = 2.2868968974132482e-14
cost at iteration 188 : 10.120939100253366
Newton iteration 0
err = 0.012304233403718912
Newton iteration 1
err = 8.168060273910043e-05
Newton iteration 2
err = 4.321172135047506e-08
Newton iteration 3
err = 2.5096560457613616e-14
cost at iteration 189 : 10.103856982279527
Newton iteration 0
err = 0.01230908650697172
Newton iteration 1
err = 8.498867215398047e-05
Newton iteration 2
err = 4.8133011040960096e-08
Newton iteration 3
err = 2.894508299023985e-14
cost at iteration 190 : 10.086807333176145
Newton iteration 0
err = 0.01231478389176927
Newton iteration 1
err = 8.852869677817834e-05
Newton iteration 2
err = 5.373772185714223e-08
Newton iteration 3
err = 3.384465764475792e-14
cost at iteration 191 : 10.069790206576535
Newton iteration 0
err = 0.0123213788309675
Newton iteration 1
err = 9.232249433997831e-05
Newton iteration 2
err = 6.013895852418721e-08
Newton iteration 3
err = 4.087117127122929e-14
cost at iteration 192 : 10.052805660865943
Newton iteration 0
err = 0.012328926440053645
Newton iteration 1
err = 9.639508034211186e-05
Newton iteration 2
err = 6.747654059348854e-08
Newton iteration 3
err = 5.0028235460533076e-14
cost at iteration 193 : 10.035853759502295
Newton iteration 0
err = 0.012337483556003531
Newton iteration 1
err = 0.0001007755028593423
Newton iteration 2
err = 7.592819852547094e-08
Newton iteration 3
err = 6.165859050605587e-14
cost at iteration 194 : 10.018934571347906
Newton iteration 0
err = 0.012347108721927739
Newton iteration 1
err = 0.00010549800964131747
Newton iteration 2
err = 8.572656006446195e-08
Newton iteration 3
err = 7.766143612890255e-14
cost at iteration 195 : 10.002048171014076
Newton iteration 0
err = 0.012357862361309683
Newton iteration 1
err = 0.0001106035672606993
Newton iteration 2
err = 9.718389217316569e-08
Newton iteration 3
err = 9.941112682673601e-14
cost at iteration 196 : 9.985194639220616
Newton iteration 0
err = 0.012369807247722811
Newton iteration 1
err = 0.00011614165387388719
Newton iteration 2
err = 1.1072622013189192e-07
Newton iteration 3
err = 1.291765036593663e-13
cost at iteration 197 : 9.968374063176789
Newton iteration 0
err = 0.012383009373436452
Newton iteration 1
err = 0.00012217207389195398
Newton iteration 2
err = 1.2693727999214811e-07
Newton iteration 3
err = 1.7158506037343033e-13
cost at iteration 198 : 9.951586536990716
Newton iteration 0
err = 0.012397539255241416
Newton iteration 1
err = 0.00012876628622800467
Newton iteration 2
err = 1.4661056483943413e-07
Newton iteration 3
err = 2.3486002309556344e-13
cost at iteration 199 : 9.934832162119742
Newton iteration 0
err = 0.012413473528612424
Newton iteration 1
err = 0.00013600741548291289
Newton iteration 2
err = 1.7080467304352815e-07
Newton iteration 3
err = 3.318844884080939e-13
cost at iteration 200 : 9.918111047873042
Newton iteration 0
err = 0.012430896288924097
Newton iteration 1
err = 0.00014398776886633165
Newton iteration 2
err = 2.008938217699363e-07
Newton iteration 3
err = 4.834922260035211e-13
cost at iteration 201 : 9.901423311976957
Newton iteration 0
err = 0.012449898936251718
Newton iteration 1
err = 0.0001528023314832948
Newton iteration 2
err = 2.3860148924241467e-07
Newton iteration 3
err = 7.242987490647969e-13
cost at iteration 202 : 9.88476908119633
Newton iteration 0
err = 0.01247057614923254
Newton iteration 1
err = 0.00016253626045895244
Newton iteration 2
err = 2.8599747434497306e-07
Newton iteration 3
err = 1.1084882808312597e-12
cost at iteration 203 : 9.868148491975601
Newton iteration 0
err = 0.012493013915133206
Newton iteration 1
err = 0.00017324352100803514
Newton iteration 2
err = 3.4541831917974555e-07
Newton iteration 3
err = 1.7247213047934312e-12
cost at iteration 204 : 9.85156169099608
Newton iteration 0
err = 0.012517263083462142
Newton iteration 1
err = 0.00018491197370533902
Newton iteration 2
err = 4.192219678381694e-07
Newton iteration 3
err = 2.7070359429132887e-12
cost at iteration 205 : 9.83500883543544
Newton iteration 0
err = 0.0125432883813894
Newton iteration 1
err = 0.00019740701464829262
Newton iteration 2
err = 5.091865962984497e-07
Newton iteration 3
err = 4.251363827784453e-12
cost at iteration 206 : 9.81849009252109
Newton iteration 0
err = 0.012570877816956915
Newton iteration 1
err = 0.00021038189387592853
Newton iteration 2
err = 6.151900786917023e-07
Newton iteration 3
err = 6.609292927535375e-12
cost at iteration 207 : 9.802005637667902
Newton iteration 0
err = 0.012599490802467373
Newton iteration 1
err = 0.00022314209699310388
Newton iteration 2
err = 7.326004776070764e-07
Newton iteration 3
err = 9.988610152632147e-12
cost at iteration 208 : 9.785555650020397
Newton iteration 0
err = 0.01262801715316255
Newton iteration 1
err = 0.00023446765222186924
Newton iteration 2
err = 8.479037836207625e-07
Newton iteration 3
err = 1.4218161692357797e-11
cost at iteration 209 : 9.769140303575066
Newton iteration 0
err = 0.012654422777052572
Newton iteration 1
err = 0.00024247016156277073
Newton iteration 2
err = 9.336824296693441e-07
Newton iteration 3
err = 1.8102876993505744e-11
cost at iteration 210 : 9.752759751352622
Newton iteration 0
err = 0.012675300323488287
Newton iteration 1
err = 0.000244770602460703
Newton iteration 2
err = 9.492701950481514e-07
Newton iteration 3
err = 1.9129926515377813e-11
cost at iteration 211 : 9.736414099903048
Newton iteration 0
err = 0.012685486689362702
Newton iteration 1
err = 0.0002396685808087788
Newton iteration 2
err = 8.605358063796654e-07
Newton iteration 3
err = 1.5352502899024855e-11
cost at iteration 212 : 9.720103373262262
Newton iteration 0
err = 0.012678212691960454
Newton iteration 1
err = 0.00022900847970583473
Newton iteration 2
err = 6.840324732515209e-07
Newton iteration 3
err = 8.805191704451343e-12
cost at iteration 213 : 9.703827471531778
Newton iteration 0
err = 0.012646541206614497
Newton iteration 1
err = 0.0002208957185917427
Newton iteration 2
err = 5.22966689346756e-07
Newton iteration 3
err = 4.1416630326860395e-12
cost at iteration 214 : 9.687586139223368
Newton iteration 0
err = 0.012586392844751335
Newton iteration 1
err = 0.00022450335903146442
Newton iteration 2
err = 5.211695561220227e-07
Newton iteration 3
err = 3.8258944850706094e-12
cost at iteration 215 : 9.671378963930431
Newton iteration 0
err = 0.012499633096018899
Newton iteration 1
err = 0.0002372786082359165
Newton iteration 2
err = 6.614704289140628e-07
Newton iteration 3
err = 7.14451032496481e-12
cost at iteration 216 : 9.655205414130458
Newton iteration 0
err = 0.012394274759087597
Newton iteration 1
err = 0.00024707393778918864
Newton iteration 2
err = 7.887869700017052e-07
Newton iteration 3
err = 1.08018893333042e-11
cost at iteration 217 : 9.639064900436525
Newton iteration 0
err = 0.012280998957710105
Newton iteration 1
err = 0.00024622226405694047
Newton iteration 2
err = 8.147828082753429e-07
Newton iteration 3
err = 1.1504768130178703e-11
cost at iteration 218 : 9.62295683269248
Newton iteration 0
err = 0.012168955736803157
Newton iteration 1
err = 0.00023513476937008805
Newton iteration 2
err = 7.456842229774209e-07
Newton iteration 3
err = 9.34197911162977e-12
cost at iteration 219 : 9.606880656389658
Newton iteration 0
err = 0.012063789570376734
Newton iteration 1
err = 0.00021784170102066448
Newton iteration 2
err = 6.284615085426893e-07
Newton iteration 3
err = 6.371559931310892e-12
cost at iteration 220 : 9.590835868498706
Newton iteration 0
err = 0.01196793615792913
Newton iteration 1
err = 0.00019822397766633868
Newton iteration 2
err = 5.039291793968795e-07
Newton iteration 3
err = 3.926784384183423e-12
cost at iteration 221 : 9.574822020125978
Newton iteration 0
err = 0.01188177468161272
Newton iteration 1
err = 0.0001787905223304271
Newton iteration 2
err = 3.934708505207827e-07
Newton iteration 3
err = 2.294102495784768e-12
cost at iteration 222 : 9.558838712805453
Newton iteration 0
err = 0.011804654631116875
Newton iteration 1
err = 0.00016080852190692008
Newton iteration 2
err = 3.0377791865066424e-07
Newton iteration 3
err = 1.3102969066548512e-12
cost at iteration 223 : 9.542885592522522
Newton iteration 0
err = 0.011735537866697104
Newton iteration 1
err = 0.0001447521705739963
Newton iteration 2
err = 2.341364977884601e-07
Newton iteration 3
err = 7.454128240984134e-13
cost at iteration 224 : 9.526962343458358
Newton iteration 0
err = 0.011673328150561187
Newton iteration 1
err = 0.00013067310191967093
Newton iteration 2
err = 1.8119767720502107e-07
Newton iteration 3
err = 4.274975110605572e-13
cost at iteration 225 : 9.51106868226787
Newton iteration 0
err = 0.011617014009844835
Newton iteration 1
err = 0.00011842963885865633
Newton iteration 2
err = 1.4126355874317424e-07
Newton iteration 3
err = 2.4905222740444714e-13
cost at iteration 226 : 9.495204353139398
Newton iteration 0
err = 0.011565716250459919
Newton iteration 1
err = 0.00010781010052778667
Newton iteration 2
err = 1.1113371406889459e-07
Newton iteration 3
err = 1.4827984752100391e-13
cost at iteration 227 : 9.479369123642961
Newton iteration 0
err = 0.01151869208482935
Newton iteration 1
err = 9.859280053177322e-05
Newton iteration 2
err = 8.829261205829106e-08
Newton iteration 3
err = 9.040825947731074e-14
cost at iteration 228 : 9.463562781282617
Newton iteration 0
err = 0.011475322082215923
Newton iteration 1
err = 9.05721292484673e-05
Newton iteration 2
err = 7.084949711119817e-08
Newton iteration 3
err = 5.723896100731357e-14
cost at iteration 229 : 9.447785130644588
Newton iteration 0
err = 0.011435092112948192
Newton iteration 1
err = 8.356771945999547e-05
Newton iteration 2
err = 5.741230685377768e-08
Newton iteration 3
err = 3.877100067040091e-14
cost at iteration 230 : 9.432035991037212
Newton iteration 0
err = 0.011397575470893687
Newton iteration 1
err = 7.742572790745611e-05
Newton iteration 2
err = 4.6964240445235616e-08
Newton iteration 3
err = 2.772227345326458e-14
cost at iteration 231 : 9.41631519453251
Newton iteration 0
err = 0.011362417083487081
Newton iteration 1
err = 7.201678665348485e-05
Newton iteration 2
err = 3.876312563224858e-08
Newton iteration 3
err = 2.2218246478313102e-14
cost at iteration 232 : 9.40062258433328
Newton iteration 0
err = 0.011329320268357642
Newton iteration 1
err = 6.723281973122311e-05
Newton iteration 2
err = 3.226534187028741e-08
Newton iteration 3
err = 1.9198398906464284e-14
cost at iteration 233 : 9.38495801340529
Newton iteration 0
err = 0.01129803590576554
Newton iteration 1
err = 6.298372621363269e-05
Newton iteration 2
err = 2.707030943748325e-08
Newton iteration 3
err = 1.9041089393383047e-14
cost at iteration 234 : 9.36932134332474
Newton iteration 0
err = 0.01126835369522867
Newton iteration 1
err = 5.919434424469781e-05
Newton iteration 2
err = 2.2880660244994022e-08
Newton iteration 3
err = 1.7097884672059588e-14
cost at iteration 235 : 9.35371244330148
Newton iteration 0
err = 0.01124009513181097
Newton iteration 1
err = 5.580183218689942e-05
Newton iteration 2
err = 1.9473834896623463e-08
Newton iteration 3
err = 1.677509532865484e-14
cost at iteration 236 : 9.338131189347598
Newton iteration 0
err = 0.011213107869280954
Newton iteration 1
err = 5.275347774490113e-05
Newton iteration 2
err = 1.6681841626400464e-08
Newton iteration 3
err = 1.6964672898322622e-14
cost at iteration 237 : 9.322577463564683
Newton iteration 0
err = 0.011187261189272615
Newton iteration 1
err = 5.0004894921011636e-05
Newton iteration 2
err = 1.4376785249343987e-08
Newton iteration 3
err = 1.6444674769673342e-14
cost at iteration 238 : 9.307051153529985
Newton iteration 0
err = 0.011162442348189949
Newton iteration 1
err = 4.7518552680031076e-05
Newton iteration 2
err = 1.2460468455832e-08
Newton iteration 3
err = 1.6081874153370927e-14
cost at iteration 239 : 9.291552151765647
Newton iteration 0
err = 0.011138553619880219
Newton iteration 1
err = 4.52625788364844e-05
Newton iteration 2
err = 1.0856870145464966e-08
Newton iteration 3
err = 1.6909940691669037e-14
cost at iteration 240 : 9.27608035527692
Newton iteration 0
err = 0.011115509890539213
Newton iteration 1
err = 4.320978850616844e-05
Newton iteration 2
err = 9.506662684896346e-09
cost at iteration 241 : 9.260635665148042
Newton iteration 0
err = 0.01109323669328479
Newton iteration 1
err = 4.133689418979085e-05
Newton iteration 2
err = 8.363190255496531e-09
cost at iteration 242 : 9.245217986188868
Newton iteration 0
err = 0.011071668593876668
Newton iteration 1
err = 3.9623862108218774e-05
Newton iteration 2
err = 7.389490675793974e-09
cost at iteration 243 : 9.229827226622453
Newton iteration 0
err = 0.011050747858202132
Newton iteration 1
err = 3.8053386170797945e-05
Newton iteration 2
err = 6.556072433715143e-09
cost at iteration 244 : 9.214463297808825
Newton iteration 0
err = 0.011030423346911684
Newton iteration 1
err = 3.6610456584234745e-05
Newton iteration 2
err = 5.839248769735352e-09
cost at iteration 245 : 9.199126114000364
Newton iteration 0
err = 0.011010649593856916
Newton iteration 1
err = 3.528200476058788e-05
Newton iteration 2
err = 5.219869923687483e-09
cost at iteration 246 : 9.183815592122667
Newton iteration 0
err = 0.010991386034311827
Newton iteration 1
err = 3.405660986798309e-05
Newton iteration 2
err = 4.682362556963818e-09
cost at iteration 247 : 9.168531651579954
Newton iteration 0
err = 0.010972596355459385
Newton iteration 1
err = 3.292425537989784e-05
Newton iteration 2
err = 4.21398565053313e-09
cost at iteration 248 : 9.153274214078976
Newton iteration 0
err = 0.01095424794744804
Newton iteration 1
err = 3.1876126267885305e-05
Newton iteration 2
err = 3.804259953751312e-09
cost at iteration 249 : 9.138043203472176
Newton iteration 0
err = 0.010936311437394968
Newton iteration 1
err = 3.090443940851102e-05
Newton iteration 2
err = 3.4445183710506476e-09
cost at iteration 250 : 9.122838545614764
Newton iteration 0
err = 0.010918760292181086
Newton iteration 1
err = 3.0002301188172677e-05
Newton iteration 2
err = 3.1275580074778088e-09
cost at iteration 251 : 9.107660168236208
Newton iteration 0
err = 0.010901570478566271
Newton iteration 1
err = 2.9163587494163017e-05
Newton iteration 2
err = 2.847365996770924e-09
cost at iteration 252 : 9.092508000823408
Newton iteration 0
err = 0.01088472017117744
Newton iteration 1
err = 2.8382842196245974e-05
Newton iteration 2
err = 2.5988969336474037e-09
cost at iteration 253 : 9.077381974514356
Newton iteration 0
err = 0.010868189500704546
Newton iteration 1
err = 2.7655190947172904e-05
Newton iteration 2
err = 2.3778995973321764e-09
cost at iteration 254 : 9.062282022001302
Newton iteration 0
err = 0.010851960336147954
Newton iteration 1
err = 2.6976267739107517e-05
Newton iteration 2
err = 2.1807802599890927e-09
cost at iteration 255 : 9.047208077442363
Newton iteration 0
err = 0.010836016095607323
Newton iteration 1
err = 2.6342152106898755e-05
Newton iteration 2
err = 2.0044832744903863e-09
cost at iteration 256 : 9.03216007638003
Newton iteration 0
err = 0.010820341581602564
Newton iteration 1
err = 2.5749315261149977e-05
Newton iteration 2
err = 1.8464049019669568e-09
cost at iteration 257 : 9.017137955667316
Newton iteration 0
err = 0.010804922837127553
Newton iteration 1
err = 2.5194573725586567e-05
Newton iteration 2
err = 1.7043189951899303e-09
cost at iteration 258 : 9.002141653399084
Newton iteration 0
err = 0.010789747019568079
Newton iteration 1
err = 2.4675049307022844e-05
Newton iteration 2
err = 1.5763115263483126e-09
cost at iteration 259 : 8.987171108848893
Newton iteration 0
err = 0.010774802289776652
Newton iteration 1
err = 2.4188134439868054e-05
Newton iteration 2
err = 1.4607360117503682e-09
cost at iteration 260 : 8.972226262411315
Newton iteration 0
err = 0.010760077714450814
Newton iteration 1
err = 2.3731462078858236e-05
Newton iteration 2
err = 1.3561691815687223e-09
cost at iteration 261 : 8.957307055547618
Newton iteration 0
err = 0.010745563179758718
Newton iteration 1
err = 2.3302879477351736e-05
Newton iteration 2
err = 1.2613768675618173e-09
cost at iteration 262 : 8.94241343073645
Newton iteration 0
err = 0.010731249314932488
Newton iteration 1
err = 2.2900425291110337e-05
Newton iteration 2
err = 1.1752881634939582e-09
cost at iteration 263 : 8.927545331427662
Newton iteration 0
err = 0.010717127424264446
Newton iteration 1
err = 2.2522309516230512e-05
Newton iteration 2
err = 1.0969672547159895e-09
cost at iteration 264 : 8.912702701999114
Newton iteration 0
err = 0.010703189426732842
Newton iteration 1
err = 2.2166895884663442e-05
Newton iteration 2
err = 1.0255986352976929e-09
cost at iteration 265 : 8.897885487717275
Newton iteration 0
err = 0.010689427802012782
Newton iteration 1
err = 2.1832686362993954e-05
Newton iteration 2
err = 9.604673230285925e-10
cost at iteration 266 : 8.883093634699643
Newton iteration 0
err = 0.010675835542283288
Newton iteration 1
err = 2.151830747443586e-05
Newton iteration 2
err = 9.009460237232743e-10
cost at iteration 267 : 8.868327089880596
Newton iteration 0
err = 0.010662406108988382
Newton iteration 1
err = 2.1222498206759814e-05
Newton iteration 2
err = 8.464810795787929e-10
cost at iteration 268 : 8.853585800978708
Newton iteration 0
err = 0.01064913339404595
Newton iteration 1
err = 2.0944099282605742e-05
Newton iteration 2
err = 7.96586408505783e-10
cost at iteration 269 : 8.838869716466716
Newton iteration 0
err = 0.010636011684949834
Newton iteration 1
err = 2.0682043641563114e-05
Newton iteration 2
err = 7.508321746296746e-10
cost at iteration 270 : 8.824178785543264
Newton iteration 0
err = 0.010623035633293374
Newton iteration 1
err = 2.0435347951697033e-05
Newton iteration 2
err = 7.088359571098894e-10
cost at iteration 271 : 8.809512958106332
Newton iteration 0
err = 0.010610200226403915
Newton iteration 1
err = 2.0203105050455344e-05
Newton iteration 2
err = 6.702615984309886e-10
cost at iteration 272 : 8.794872184728579
Newton iteration 0
err = 0.0105975007615993
Newton iteration 1
err = 1.9984477172279958e-05
Newton iteration 2
err = 6.348078483870471e-10
cost at iteration 273 : 8.78025641663373
Newton iteration 0
err = 0.01058493282296307
Newton iteration 1
err = 1.9778689900939212e-05
Newton iteration 2
err = 6.022078684365181e-10
cost at iteration 274 : 8.765665605675103
Newton iteration 0
err = 0.010572492260166036
Newton iteration 1
err = 1.958502672355081e-05
Newton iteration 2
err = 5.722250135421281e-10
cost at iteration 275 : 8.751099704314528
Newton iteration 0
err = 0.01056017516934909
Newton iteration 1
err = 1.9402824156208822e-05
Newton iteration 2
err = 5.446455358956044e-10
cost at iteration 276 : 8.736558665603429
Newton iteration 0
err = 0.010547977875574177
Newton iteration 1
err = 1.923146733222954e-05
Newton iteration 2
err = 5.192814229874704e-10
cost at iteration 277 : 8.722042443164213
Newton iteration 0
err = 0.010535896916934864
Newton iteration 1
err = 1.907038604805638e-05
Newton iteration 2
err = 4.959608447190959e-10
cost at iteration 278 : 8.707550991173367
Newton iteration 0
err = 0.010523929030007702
Newton iteration 1
err = 1.891905116748715e-05
Newton iteration 2
err = 4.745318490291801e-10
cost at iteration 279 : 8.693084264345016
Newton iteration 0
err = 0.010512071136594888
Newton iteration 1
err = 1.8776971385383153e-05
Newton iteration 2
err = 4.5486010156903315e-10
cost at iteration 280 : 8.678642217915655
Newton iteration 0
err = 0.010500320331534438
Newton iteration 1
err = 1.8643690290093208e-05
Newton iteration 2
err = 4.368194105503116e-10
cost at iteration 281 : 8.66422480762989
Newton iteration 0
err = 0.010488673871558945
Newton iteration 1
err = 1.8518783697116382e-05
Newton iteration 2
err = 4.203002790227481e-10
cost at iteration 282 : 8.649831989726483
Newton iteration 0
err = 0.010477129165134757
Newton iteration 1
err = 1.8401857226762575e-05
Newton iteration 2
err = 4.0520430764477503e-10
cost at iteration 283 : 8.635463720925639
Newton iteration 0
err = 0.010465683763020274
Newton iteration 1
err = 1.8292544107491393e-05
Newton iteration 2
err = 3.9143902275094776e-10
cost at iteration 284 : 8.621119958416713
Newton iteration 0
err = 0.010454335349679222
Newton iteration 1
err = 1.819050316594789e-05
Newton iteration 2
err = 3.78923333266539e-10
cost at iteration 285 : 8.606800659846714
Newton iteration 0
err = 0.010443081735345167
Newton iteration 1
err = 1.8095417004308304e-05
Newton iteration 2
err = 3.6758224293467255e-10
cost at iteration 286 : 8.59250578330947
Newton iteration 0
err = 0.010431920848720777
Newton iteration 1
err = 1.8006990327490742e-05
Newton iteration 2
err = 3.573469164011389e-10
cost at iteration 287 : 8.578235287335051
Newton iteration 0
err = 0.010420850730261727
Newton iteration 1
err = 1.792494842911105e-05
Newton iteration 2
err = 3.481577536261013e-10
cost at iteration 288 : 8.563989130880133
Newton iteration 0
err = 0.010409869525946783
Newton iteration 1
err = 1.7849035787593058e-05
Newton iteration 2
err = 3.399536765081022e-10
cost at iteration 289 : 8.549767273318723
Newton iteration 0
err = 0.010398975481580032
Newton iteration 1
err = 1.7779014803610412e-05
Newton iteration 2
err = 3.326861907013549e-10
cost at iteration 290 : 8.535569674433505
Newton iteration 0
err = 0.010388166937496652
Newton iteration 1
err = 1.771466463184349e-05
Newton iteration 2
err = 3.263041792635424e-10
cost at iteration 291 : 8.521396294407083
Newton iteration 0
err = 0.010377442323663458
Newton iteration 1
err = 1.7655780113992115e-05
Newton iteration 2
err = 3.2076283763241824e-10
cost at iteration 292 : 8.507247093814817
Newton iteration 0
err = 0.01036680015521094
Newton iteration 1
err = 1.7602170797941442e-05
Newton iteration 2
err = 3.1602080798532066e-10
cost at iteration 293 : 8.49312203361659
Newton iteration 0
err = 0.010356239028215645
Newton iteration 1
err = 1.755366004547435e-05
Newton iteration 2
err = 3.120400275841814e-10
cost at iteration 294 : 8.479021075150364
Newton iteration 0
err = 0.010345757615867103
Newton iteration 1
err = 1.751008421189567e-05
Newton iteration 2
err = 3.0878307932244484e-10
cost at iteration 295 : 8.464944180125306
Newton iteration 0
err = 0.010335354664924816
Newton iteration 1
err = 1.747129189104011e-05
Newton iteration 2
err = 3.062159957712191e-10
cost at iteration 296 : 8.450891310615278
Newton iteration 0
err = 0.010325028992358562
Newton iteration 1
err = 1.74371432240774e-05
Newton iteration 2
err = 3.0430846535554714e-10
cost at iteration 297 : 8.436862429053274
Newton iteration 0
err = 0.010314779482390997
Newton iteration 1
err = 1.7407509270039755e-05
Newton iteration 2
err = 3.0303159886840744e-10
cost at iteration 298 : 8.422857498225401
Newton iteration 0
err = 0.010304605083590454
Newton iteration 1
err = 1.7382271427044022e-05
Newton iteration 2
err = 3.0235782817468916e-10
cost at iteration 299 : 8.408876481265745
Newton iteration 0
err = 0.01029450480635507
Newton iteration 1
err = 1.7361320901541812e-05
Newton iteration 2
err = 3.022615727882552e-10
cost at iteration 300 : 8.394919341651374
Newton iteration 0
err = 0.010284477720448798
Newton iteration 1
err = 1.734455821813535e-05
Newton iteration 2
err = 3.027201948375371e-10
cost at iteration 301 : 8.380986043197575
Newton iteration 0
err = 0.010274522952835
Newton iteration 1
err = 1.7331892789270654e-05
Newton iteration 2
err = 3.0371329287697595e-10
cost at iteration 302 : 8.367076550053328
Newton iteration 0
err = 0.010264639685620037
Newton iteration 1
err = 1.7323242502307707e-05
Newton iteration 2
err = 3.0522088788762664e-10
cost at iteration 303 : 8.353190826697203
Newton iteration 0
err = 0.010254827154192117
Newton iteration 1
err = 1.7318533354047955e-05
Newton iteration 2
err = 3.072268301869809e-10
cost at iteration 304 : 8.339328837933415
Newton iteration 0
err = 0.010245084645508936
Newton iteration 1
err = 1.731769912376214e-05
Newton iteration 2
err = 3.097146622120062e-10
cost at iteration 305 : 8.325490548888231
Newton iteration 0
err = 0.010235411496504617
Newton iteration 1
err = 1.73206810647966e-05
Newton iteration 2
err = 3.126729693447604e-10
cost at iteration 306 : 8.311675925006464
Newton iteration 0
err = 0.010225807092709474
Newton iteration 1
err = 1.732742763759152e-05
Newton iteration 2
err = 3.1608944962370677e-10
cost at iteration 307 : 8.29788493204824
Newton iteration 0
err = 0.010216270866902069
Newton iteration 1
err = 1.733789427018766e-05
Newton iteration 2
err = 3.1995462744244705e-10
cost at iteration 308 : 8.284117536086221
Newton iteration 0
err = 0.01020680229796515
Newton iteration 1
err = 1.7352043132509136e-05
Newton iteration 2
err = 3.242600057048941e-10
cost at iteration 309 : 8.270373703502608
Newton iteration 0
err = 0.010197400909811584
Newton iteration 1
err = 1.7369842956466188e-05
Newton iteration 2
err = 3.290010343147126e-10
cost at iteration 310 : 8.256653400986828
Newton iteration 0
err = 0.010188066270470756
Newton iteration 1
err = 1.7391268860904182e-05
Newton iteration 2
err = 3.3417325115878346e-10
cost at iteration 311 : 8.242956595533146
Newton iteration 0
err = 0.01017879799123705
Newton iteration 1
err = 1.7416302210589806e-05
Newton iteration 2
err = 3.397748357500292e-10
cost at iteration 312 : 8.22928325443864
Newton iteration 0
err = 0.010169595725921105
Newton iteration 1
err = 1.7444930497827835e-05
Newton iteration 2
err = 3.4580517083605715e-10
cost at iteration 313 : 8.215633345301129
Newton iteration 0
err = 0.010160459170286292
Newton iteration 1
err = 1.7477147234072907e-05
Newton iteration 2
err = 3.5226334855039417e-10
cost at iteration 314 : 8.202006836017606
Newton iteration 0
err = 0.010151388061481605
Newton iteration 1
err = 1.751295187589034e-05
Newton iteration 2
err = 3.591553442215056e-10
cost at iteration 315 : 8.188403694782872
Newton iteration 0
err = 0.010142382177609063
Newton iteration 1
err = 1.755234976180241e-05
Newton iteration 2
err = 3.664843495301553e-10
cost at iteration 316 : 8.174823890087946
Newton iteration 0
err = 0.010133441337422337
Newton iteration 1
err = 1.7595352071572173e-05
Newton iteration 2
err = 3.7425804862622527e-10
cost at iteration 317 : 8.161267390719201
Newton iteration 0
err = 0.010124565400060985
Newton iteration 1
err = 1.7641975803739507e-05
Newton iteration 2
err = 3.824829122975472e-10
cost at iteration 318 : 8.147734165757226
Newton iteration 0
err = 0.010115754264902079
Newton iteration 1
err = 1.769224377640475e-05
Newton iteration 2
err = 3.9117123453281395e-10
cost at iteration 319 : 8.13422418457646
Newton iteration 0
err = 0.010107007871518296
Newton iteration 1
err = 1.7746184635732924e-05
Newton iteration 2
err = 4.0033125469969907e-10
cost at iteration 320 : 8.12073741684424
Newton iteration 0
err = 0.010098326199735079
Newton iteration 1
err = 1.7803832894577003e-05
Newton iteration 2
err = 4.0998101829613647e-10
cost at iteration 321 : 8.107273832520715
Newton iteration 0
err = 0.010089709269664026
Newton iteration 1
err = 1.7865228983606987e-05
Newton iteration 2
err = 4.2013163852567323e-10
cost at iteration 322 : 8.09383340185873
Newton iteration 0
err = 0.010081157142043629
Newton iteration 1
err = 1.793041932269758e-05
Newton iteration 2
err = 4.308036900990181e-10
cost at iteration 323 : 8.080416095403713
Newton iteration 0
err = 0.010072669918452992
Newton iteration 1
err = 1.7999456406894193e-05
Newton iteration 2
err = 4.420146513294219e-10
cost at iteration 324 : 8.067021883994007
Newton iteration 0
err = 0.010064247741761732
Newton iteration 1
err = 1.8072398920311412e-05
Newton iteration 2
err = 4.5378584053413096e-10
cost at iteration 325 : 8.05365073876125
Newton iteration 0
err = 0.010055890796595975
Newton iteration 1
err = 1.8149311860283086e-05
Newton iteration 2
err = 4.661409257472278e-10
cost at iteration 326 : 8.040302631131096
Newton iteration 0
err = 0.010047599309983956
Newton iteration 1
err = 1.8230266682994816e-05
Newton iteration 2
err = 4.791053903184e-10
cost at iteration 327 : 8.026977532823976
Newton iteration 0
err = 0.010039373551996699
Newton iteration 1
err = 1.8315341474506054e-05
Newton iteration 2
err = 4.927075792914001e-10
cost at iteration 328 : 8.013675415856108
Newton iteration 0
err = 0.01003121383655408
Newton iteration 1
err = 1.8404621134716773e-05
Newton iteration 2
err = 5.069763919625628e-10
cost at iteration 329 : 8.000396252540833
Newton iteration 0
err = 0.010023120522354248
Newton iteration 1
err = 1.8498197590319304e-05
Newton iteration 2
err = 5.219457044743445e-10
cost at iteration 330 : 7.98714001548986
Newton iteration 0
err = 0.010015094013871174
Newton iteration 1
err = 1.8596170023404392e-05
Newton iteration 2
err = 5.376507411111589e-10
cost at iteration 331 : 7.97390667761529
Newton iteration 0
err = 0.010007134762483173
Newton iteration 1
err = 1.8698645123893623e-05
Newton iteration 2
err = 5.54130822123991e-10
cost at iteration 332 : 7.960696212131149
Newton iteration 0
err = 0.009999243267681233
Newton iteration 1
err = 1.8805737368903218e-05
Newton iteration 2
err = 5.714264511013578e-10
cost at iteration 333 : 7.947508592555636
Newton iteration 0
err = 0.00999142007849219
Newton iteration 1
err = 1.8917569326535872e-05
Newton iteration 2
err = 5.895813670581057e-10
cost at iteration 334 : 7.934343792713496
Newton iteration 0
err = 0.009983665794949054
Newton iteration 1
err = 1.9034271976626107e-05
Newton iteration 2
err = 6.086474896437384e-10
cost at iteration 335 : 7.921201786738643
Newton iteration 0
err = 0.009975981069690464
Newton iteration 1
err = 1.9155985070018674e-05
Newton iteration 2
err = 6.28673868835676e-10
cost at iteration 336 : 7.908082549076849
Newton iteration 0
err = 0.009968366609787881
Newton iteration 1
err = 1.9282857508541858e-05
Newton iteration 2
err = 6.497169604145605e-10
cost at iteration 337 : 7.89498605448894
Newton iteration 0
err = 0.009960823178635428
Newton iteration 1
err = 1.9415047753028568e-05
Newton iteration 2
err = 6.718383401351325e-10
cost at iteration 338 : 7.881912278054166
Newton iteration 0
err = 0.009953351598019265
Newton iteration 1
err = 1.955272426425585e-05
Newton iteration 2
err = 6.951029043824951e-10
cost at iteration 339 : 7.8688611951736505
Newton iteration 0
err = 0.009945952750389186
Newton iteration 1
err = 1.969606597618108e-05
Newton iteration 2
err = 7.195811061109158e-10
cost at iteration 340 : 7.855832781574667
Newton iteration 0
err = 0.009938627581249238
Newton iteration 1
err = 1.9845262796823227e-05
Newton iteration 2
err = 7.45348866507181e-10
cost at iteration 341 : 7.842827013314456
Newton iteration 0
err = 0.0099313771018026
Newton iteration 1
err = 2.0000516146452823e-05
Newton iteration 2
err = 7.724873027428906e-10
cost at iteration 342 : 7.829843866785014
Newton iteration 0
err = 0.009924202391663162
Newton iteration 1
err = 2.016203953044025e-05
Newton iteration 2
err = 8.010859720016777e-10
cost at iteration 343 : 7.81688331871779
Newton iteration 0
err = 0.00991710460196803
Newton iteration 1
err = 2.033005914524305e-05
Newton iteration 2
err = 8.312397433966234e-10
cost at iteration 344 : 7.803945346188948
Newton iteration 0
err = 0.009910084958523632
Newton iteration 1
err = 2.0504814526263714e-05
Newton iteration 2
err = 8.630503207097703e-10
cost at iteration 345 : 7.791029926624875
Newton iteration 0
err = 0.009903144765262523
Newton iteration 1
err = 2.0686559229127302e-05
Newton iteration 2
err = 8.96629958404208e-10
cost at iteration 346 : 7.778137037808137
Newton iteration 0
err = 0.009896285407998912
Newton iteration 1
err = 2.0875561553211083e-05
Newton iteration 2
err = 9.320982069230523e-10
cost at iteration 347 : 7.765266657883698
Newton iteration 0
err = 0.009889508358301418
Newton iteration 1
err = 2.107210529907077e-05
Newton iteration 2
err = 9.695859122376951e-10
cost at iteration 348 : 7.7524187653657135
Newton iteration 0
err = 0.009882815177749332
Newton iteration 1
err = 2.1276490575419525e-05
Newton iteration 2
err = 1.0092307867918985e-09
cost at iteration 349 : 7.739593339144463
Newton iteration 0
err = 0.00987620752248644
Newton iteration 1
err = 2.1489034630325355e-05
Newton iteration 2
err = 1.0511865409069397e-09
cost at iteration 350 : 7.726790358494299
Newton iteration 0
err = 0.009869687147961296
Newton iteration 1
err = 2.1710072730444545e-05
Newton iteration 2
err = 1.0956166410079054e-09
cost at iteration 351 : 7.7140098030812645
Newton iteration 0
err = 0.009863255914100777
Newton iteration 1
err = 2.1939959074371744e-05
Newton iteration 2
err = 1.1426997053618905e-09
cost at iteration 352 : 7.701251652972075
Newton iteration 0
err = 0.009856915790792505
Newton iteration 1
err = 2.217906773881653e-05
Newton iteration 2
err = 1.1926271033226805e-09
cost at iteration 353 : 7.688515888642841
Newton iteration 0
err = 0.009850668863638135
Newton iteration 1
err = 2.242779364812339e-05
Newton iteration 2
err = 1.2456082019461151e-09
cost at iteration 354 : 7.6758024909888745
Newton iteration 0
err = 0.00984451734018273
Newton iteration 1
err = 2.2686553576899238e-05
Newton iteration 2
err = 1.3018685072092917e-09
cost at iteration 355 : 7.663111441334938
Newton iteration 0
err = 0.009838463556459646
Newton iteration 1
err = 2.2955787158191435e-05
Newton iteration 2
err = 1.36165406582935e-09
cost at iteration 356 : 7.65044272144599
Newton iteration 0
err = 0.00983250998390936
Newton iteration 1
err = 2.3235957907955807e-05
Newton iteration 2
err = 1.4252301950942225e-09
cost at iteration 357 : 7.6377963135387805
Newton iteration 0
err = 0.009826659236764242
Newton iteration 1
err = 2.352755423391383e-05
Newton iteration 2
err = 1.49288632533243e-09
cost at iteration 358 : 7.625172200293843
Newton iteration 0
err = 0.009820914079738454
Newton iteration 1
err = 2.383109042637388e-05
Newton iteration 2
err = 1.5649353565276393e-09
cost at iteration 359 : 7.6125703648686915
Newton iteration 0
err = 0.009815277436272347
Newton iteration 1
err = 2.4147107604190455e-05
Newton iteration 2
err = 1.6417168268046813e-09
cost at iteration 360 : 7.599990790911121
Newton iteration 0
err = 0.009809752397008572
Newton iteration 1
err = 2.447617458732318e-05
Newton iteration 2
err = 1.7235993217923515e-09
cost at iteration 361 : 7.587433462574027
Newton iteration 0
err = 0.009804342228818135
Newton iteration 1
err = 2.48188886797379e-05
Newton iteration 2
err = 1.8109843092377124e-09
cost at iteration 362 : 7.57489836453053
Newton iteration 0
err = 0.009799050384092726
Newton iteration 1
err = 2.5175876288010577e-05
Newton iteration 2
err = 1.9043048029681677e-09
cost at iteration 363 : 7.562385481990354
Newton iteration 0
err = 0.009793880510463161
Newton iteration 1
err = 2.5547793369221164e-05
Newton iteration 2
err = 2.0040332416835827e-09
cost at iteration 364 : 7.549894800716963
Newton iteration 0
err = 0.009788836460704373
Newton iteration 1
err = 2.5935325613858137e-05
Newton iteration 2
err = 2.1106783873768954e-09
cost at iteration 365 : 7.537426307045587
Newton iteration 0
err = 0.009783922302982055
Newton iteration 1
err = 2.633918829296073e-05
Newton iteration 2
err = 2.2247944438722137e-09
cost at iteration 366 : 7.52497998790282
Newton iteration 0
err = 0.00977914233113634
Newton iteration 1
err = 2.67601256944044e-05
Newton iteration 2
err = 2.34697987725527e-09
cost at iteration 367 : 7.512555830826592
Newton iteration 0
err = 0.009774501075013096
Newton iteration 1
err = 2.7198910006745193e-05
Newton iteration 2
err = 2.477880909865271e-09
cost at iteration 368 : 7.500153823988005
Newton iteration 0
err = 0.009770003310599204
Newton iteration 1
err = 2.7656339532215763e-05
Newton iteration 2
err = 2.6181961669149398e-09
cost at iteration 369 : 7.487773956213732
Newton iteration 0
err = 0.009765654069783313
Newton iteration 1
err = 2.8133236050470487e-05
Newton iteration 2
err = 2.768677745677034e-09
cost at iteration 370 : 7.475416217010188
Newton iteration 0
err = 0.009761458649468441
Newton iteration 1
err = 2.863044111672924e-05
Newton iteration 2
err = 2.9301326690914217e-09
cost at iteration 371 : 7.463080596588835
Newton iteration 0
err = 0.009757422619622782
Newton iteration 1
err = 2.9148811054707186e-05
Newton iteration 2
err = 3.103427406174945e-09
cost at iteration 372 : 7.45076708589259
Newton iteration 0
err = 0.009753551829872264
Newton iteration 1
err = 2.9689210334377875e-05
Newton iteration 2
err = 3.289485296903913e-09
cost at iteration 373 : 7.43847567662398
Newton iteration 0
err = 0.009749852414065207
Newton iteration 1
err = 3.0252502951988608e-05
Newton iteration 2
err = 3.4892889093365995e-09
cost at iteration 374 : 7.426206361274561
Newton iteration 0
err = 0.009746330791994353
Newton iteration 1
err = 3.083954138060381e-05
Newton iteration 2
err = 3.7038746812260253e-09
cost at iteration 375 : 7.413959133155532
Newton iteration 0
err = 0.009742993667551856
Newton iteration 1
err = 3.145115253959402e-05
Newton iteration 2
err = 3.934329998704516e-09
cost at iteration 376 : 7.4017339864302665
Newton iteration 0
err = 0.009739848021939982
Newton iteration 1
err = 3.2088120135433247e-05
Newton iteration 2
err = 4.181787107736253e-09
cost at iteration 377 : 7.389530916147467
Newton iteration 0
err = 0.009736901100909051
Newton iteration 1
err = 3.275116260255677e-05
Newton iteration 2
err = 4.447408347274089e-09
cost at iteration 378 : 7.377349918276473
Newton iteration 0
err = 0.009734160393787737
Newton iteration 1
err = 3.3440905713660976e-05
Newton iteration 2
err = 4.732369029639533e-09
cost at iteration 379 : 7.36519098974291
Newton iteration 0
err = 0.009731633602747201
Newton iteration 1
err = 3.415784876741146e-05
Newton iteration 2
err = 5.037834855387816e-09
cost at iteration 380 : 7.35305412846606
Newton iteration 0
err = 0.009729328599216587
Newton iteration 1
err = 3.490232305345404e-05
Newton iteration 2
err = 5.364930720597321e-09
cost at iteration 381 : 7.340939333396536
Newton iteration 0
err = 0.009727253364437449
Newton iteration 1
err = 3.5674441092446135e-05
Newton iteration 2
err = 5.7146938102749805e-09
cost at iteration 382 : 7.328846604554824
Newton iteration 0
err = 0.009725415910288651
Newton iteration 1
err = 3.647403490225224e-05
Newton iteration 2
err = 6.088018614060089e-09
cost at iteration 383 : 7.316775943069449
Newton iteration 0
err = 0.009723824175419567
Newton iteration 1
err = 3.730058128833815e-05
Newton iteration 2
err = 6.485585172367059e-09
cost at iteration 384 : 7.304727351215091
Newton iteration 0
err = 0.00972248589117931
Newton iteration 1
err = 3.815311192487166e-05
Newton iteration 2
err = 6.9077601301855975e-09
cost at iteration 385 : 7.292700832449385
Newton iteration 0
err = 0.009721408410165031
Newton iteration 1
err = 3.903010576142496e-05
Newton iteration 2
err = 7.354486022898466e-09
cost at iteration 386 : 7.2806963914472345
Newton iteration 0
err = 0.009720598489301785
Newton iteration 1
err = 3.992936114590643e-05
Newton iteration 2
err = 7.825133778556886e-09
cost at iteration 387 : 7.268714034132369
Newton iteration 0
err = 0.009720062017840793
Newton iteration 1
err = 4.084784502995097e-05
Newton iteration 2
err = 8.318326199454976e-09
cost at iteration 388 : 7.256753767704086
Newton iteration 0
err = 0.009719803678429164
Newton iteration 1
err = 4.178151679406261e-05
Newton iteration 2
err = 8.831727192912407e-09
cost at iteration 389 : 7.244815600657051
Newton iteration 0
err = 0.009719826528935277
Newton iteration 1
err = 4.272512477712948e-05
Newton iteration 2
err = 9.361799844870752e-09
cost at iteration 390 : 7.232899542792949
Newton iteration 0
err = 0.009720131489545437
Newton iteration 1
err = 4.367197458877147e-05
Newton iteration 2
err = 9.90352798832348e-09
cost at iteration 391 : 7.221005605220383
Newton iteration 0
err = 0.009720716719626439
Newton iteration 1
err = 4.461367011296824e-05
Newton iteration 2
err = 1.0450111266433136e-08
Newton iteration 3
err = 1.398617220183724e-14
cost at iteration 392 : 7.209133800339945
Newton iteration 0
err = 0.009721576866879317
Newton iteration 1
err = 4.553983090504712e-05
Newton iteration 2
err = 1.0992655892544408e-08
Newton iteration 3
err = 1.3763130590648683e-14
cost at iteration 393 : 7.19728414181033
Newton iteration 0
err = 0.009722702171973602
Newton iteration 1
err = 4.6437794013292374e-05
Newton iteration 2
err = 1.1519873590516903e-08
Newton iteration 3
err = 1.457052684168393e-14
cost at iteration 394 : 7.1854566444915955
Newton iteration 0
err = 0.009724077413231105
Newton iteration 1
err = 4.729231447959694e-05
Newton iteration 2
err = 1.2017844522160257e-08
Newton iteration 3
err = 1.357974096275976e-14
cost at iteration 395 : 7.173651324358728
Newton iteration 0
err = 0.009725680679794829
Newton iteration 1
err = 4.8085287474875344e-05
Newton iteration 2
err = 1.2469883301809588e-08
Newton iteration 3
err = 1.3941319377051174e-14
cost at iteration 396 : 7.161868198379745
Newton iteration 0
err = 0.009727481969459016
Newton iteration 1
err = 4.879552672915221e-05
Newton iteration 2
err = 1.2856606152104496e-08
Newton iteration 3
err = 1.4278081006324533e-14
cost at iteration 397 : 7.150107284351836
Newton iteration 0
err = 0.009729441619304462
Newton iteration 1
err = 4.9398649099038134e-05
Newton iteration 2
err = 1.3156273224041452e-08
Newton iteration 3
err = 1.3837656796887847e-14
cost at iteration 398 : 7.138368600685554
Newton iteration 0
err = 0.009731508595989277
Newton iteration 1
err = 4.986713382851552e-05
Newton iteration 2
err = 1.3345534722105363e-08
Newton iteration 3
err = 1.43752149383691e-14
cost at iteration 399 : 7.126652166130711
Newton iteration 0
err = 0.009733618698777856
Newton iteration 1
err = 5.017064703507137e-05
Newton iteration 2
err = 1.3400667598194664e-08
Newton iteration 3
err = 1.4233773671138964e-14
cost at iteration 400 : 7.114957999433765
Newton iteration 0
err = 0.009735692763459997
Newton iteration 1
err = 5.0276745788699966e-05
Newton iteration 2
err = 1.3299386142870326e-08
Newton iteration 3
err = 1.4041397769829267e-14
cost at iteration 401 : 7.103286118918815
Newton iteration 0
err = 0.009737634999680512
Newton iteration 1
err = 5.015209941563411e-05
Newton iteration 2
err = 1.3023230273862849e-08
Newton iteration 3
err = 1.3582801165669328e-14
cost at iteration 402 : 7.091636541984913
Newton iteration 0
err = 0.009739331645368735
Newton iteration 1
err = 4.9764384222212626e-05
Newton iteration 2
err = 1.2560409613699563e-08
Newton iteration 3
err = 1.369468174446171e-14
cost at iteration 403 : 7.080009284514986
Newton iteration 0
err = 0.00974065017666878
Newton iteration 1
err = 4.908501598486364e-05
Newton iteration 2
err = 1.1908833874310026e-08
Newton iteration 3
err = 1.4185136414341598e-14
cost at iteration 404 : 7.068404360195021
Newton iteration 0
err = 0.00974143936038501
Newton iteration 1
err = 4.8092874893604476e-05
Newton iteration 2
err = 1.1078844602222444e-08
Newton iteration 3
err = 1.3497578061456434e-14
cost at iteration 405 : 7.0568217797465325
Newton iteration 0
err = 0.009741530466938741
Newton iteration 1
err = 4.6779140441579204e-05
Newton iteration 2
err = 1.0095040104838732e-08
Newton iteration 3
err = 1.4647977986138938e-14
cost at iteration 406 : 7.045261550085091
Newton iteration 0
err = 0.009740739956494174
Newton iteration 1
err = 4.515327368926383e-05
Newton iteration 2
err = 8.996510004716513e-09
cost at iteration 407 : 7.033723673420393
Newton iteration 0
err = 0.009738873893465777
Newton iteration 1
err = 4.3250028071159e-05
Newton iteration 2
err = 7.834996131122178e-09
cost at iteration 408 : 7.022208146329081
Newton iteration 0
err = 0.009735734215135955
Newton iteration 1
err = 4.1137054130954686e-05
Newton iteration 2
err = 6.670923693965339e-09
cost at iteration 409 : 7.010714958833902
Newton iteration 0
err = 0.009731126777683596
Newton iteration 1
err = 3.892199379656642e-05
Newton iteration 2
err = 5.567971502616132e-09
cost at iteration 410 : 6.999244093535542
Newton iteration 0
err = 0.009724870834959463
Newton iteration 1
err = 3.6756576150920496e-05
Newton iteration 2
err = 4.5876523144821494e-09
cost at iteration 411 : 6.987795524845162
Newton iteration 0
err = 0.00971680930961486
Newton iteration 1
err = 3.48328003798658e-05
Newton iteration 2
err = 3.785731734856403e-09
cost at iteration 412 : 6.976369218367619
Newton iteration 0
err = 0.00970681894911114
Newton iteration 1
err = 3.3363641807648354e-05
Newton iteration 2
err = 3.2109520410634428e-09
cost at iteration 413 : 6.964965130480298
Newton iteration 0
err = 0.009694819290242209
Newton iteration 1
err = 3.2542166755905926e-05
Newton iteration 2
err = 2.901885747580299e-09
cost at iteration 414 : 6.953583208143355
Newton iteration 0
err = 0.009680779348461156
Newton iteration 1
err = 3.248566661107678e-05
Newton iteration 2
err = 2.8752668463732654e-09
cost at iteration 415 : 6.942223388960291
Newton iteration 0
err = 0.009664721135660708
Newton iteration 1
err = 3.319231825939447e-05
Newton iteration 2
err = 3.1148330705041826e-09
cost at iteration 416 : 6.930885601492659
Newton iteration 0
err = 0.009646719478402193
Newton iteration 1
err = 3.4541097713293145e-05
Newton iteration 2
err = 3.580748325470575e-09
cost at iteration 417 : 6.919569765812022
Newton iteration 0
err = 0.00962689809315406
Newton iteration 1
err = 3.633522985543137e-05
Newton iteration 2
err = 4.228112859701795e-09
cost at iteration 418 : 6.908275794256148
Newton iteration 0
err = 0.009605422375915721
Newton iteration 1
err = 3.835857242420371e-05
Newton iteration 2
err = 5.012828726538204e-09
cost at iteration 419 : 6.8970035923446975
Newton iteration 0
err = 0.00958248977504741
Newton iteration 1
err = 4.041702839506683e-05
Newton iteration 2
err = 5.888626635918416e-09
cost at iteration 420 : 6.885753059801368
Newton iteration 0
err = 0.009558318856725322
Newton iteration 1
err = 4.235808147410718e-05
Newton iteration 2
err = 6.805323327299809e-09
cost at iteration 421 : 6.8745240916279915
Newton iteration 0
err = 0.00953313821159469
Newton iteration 1
err = 4.4074574504422617e-05
Newton iteration 2
err = 7.710939790023761e-09
cost at iteration 422 : 6.863316579181749
Newton iteration 0
err = 0.009507176206416591
Newton iteration 1
err = 4.550041240367954e-05
Newton iteration 2
err = 8.556158587950486e-09
cost at iteration 423 : 6.852130411211048
Newton iteration 0
err = 0.009480652315767274
Newton iteration 1
err = 4.660319849368785e-05
Newton iteration 2
err = 9.299015084230928e-09
cost at iteration 424 : 6.840965474818703
Newton iteration 0
err = 0.009453770448706099
Newton iteration 1
err = 4.737634589286554e-05
Newton iteration 2
err = 9.908360999795208e-09
cost at iteration 425 : 6.82982165632848
Newton iteration 0
err = 0.009426714382113193
Newton iteration 1
err = 4.783179862966938e-05
Newton iteration 2
err = 1.0365446718949274e-08
Newton iteration 3
err = 1.4050596565451031e-14
cost at iteration 426 : 6.818698842043592
Newton iteration 0
err = 0.009399645171044123
Newton iteration 1
err = 4.7993816202270816e-05
Newton iteration 2
err = 1.0663656162068951e-08
Newton iteration 3
err = 1.3178458760927205e-14
cost at iteration 427 : 6.807596918891882
Newton iteration 0
err = 0.009372700248258958
Newton iteration 1
err = 4.7893952232831266e-05
Newton iteration 2
err = 1.080684707040156e-08
Newton iteration 3
err = 1.4111951610100419e-14
cost at iteration 428 : 6.796515774959711
Newton iteration 0
err = 0.009345993848134834
Newton iteration 1
err = 4.756719491187188e-05
Newton iteration 2
err = 1.0806899474350026e-08
Newton iteration 3
err = 1.4007377762386024e-14
cost at iteration 429 : 6.785455299922468
Newton iteration 0
err = 0.00931961837814476
Newton iteration 1
err = 4.704915065642696e-05
Newton iteration 2
err = 1.0681039611110315e-08
Newton iteration 3
err = 1.3763700315135681e-14
cost at iteration 430 : 6.774415385379738
Newton iteration 0
err = 0.009293646395272564
Newton iteration 1
err = 4.637411048976975e-05
Newton iteration 2
err = 1.0449354696564205e-08
Newton iteration 3
err = 1.403056298483153e-14
cost at iteration 431 : 6.763395925107754
Newton iteration 0
err = 0.00926813290172403
Newton iteration 1
err = 4.557382607919982e-05
Newton iteration 2
err = 1.0132746137814827e-08
Newton iteration 3
err = 1.3485393971588039e-14
cost at iteration 432 : 6.75239681523993
Newton iteration 0
err = 0.009243117741110946
Newton iteration 1
err = 4.467682952317102e-05
Newton iteration 2
err = 9.751406700657248e-09
cost at iteration 433 : 6.741417954387316
Newton iteration 0
err = 0.009218627939091718
Newton iteration 1
err = 4.3708149991526546e-05
Newton iteration 2
err = 9.323814346922373e-09
cost at iteration 434 : 6.730459243708893
Newton iteration 0
err = 0.009194679887524499
Newton iteration 1
err = 4.268930491848253e-05
Newton iteration 2
err = 8.866163044919178e-09
cost at iteration 435 : 6.719520586942183
Newton iteration 0
err = 0.009171281314503777
Newton iteration 1
err = 4.1638469036240785e-05
Newton iteration 2
err = 8.392134460874286e-09
cost at iteration 436 : 6.708601890400773
Newton iteration 0
err = 0.009148433014795456
Newton iteration 1
err = 4.0570748282000334e-05
Newton iteration 2
err = 7.912911851157044e-09
cost at iteration 437 : 6.697703062947585
Newton iteration 0
err = 0.009126130337977203
Newton iteration 1
err = 3.949850589184493e-05
Newton iteration 2
err = 7.4373476352486464e-09
cost at iteration 438 : 6.68682401594873
Newton iteration 0
err = 0.009104364445456876
Newton iteration 1
err = 3.843170444874152e-05
Newton iteration 2
err = 6.972215193619555e-09
cost at iteration 439 : 6.675964663213248
Newton iteration 0
err = 0.009083123356465854
Newton iteration 1
err = 3.737824028682545e-05
Newton iteration 2
err = 6.5225037621976025e-09
cost at iteration 440 : 6.665124920922531
Newton iteration 0
err = 0.009062392806559407
Newton iteration 1
err = 3.634425600688165e-05
Newton iteration 2
err = 6.0917085020167406e-09
cost at iteration 441 : 6.654304707552929
Newton iteration 0
err = 0.009042156943816111
Newton iteration 1
err = 3.53344234933217e-05
Newton iteration 2
err = 5.682110899608971e-09
cost at iteration 442 : 6.643503943793962
Newton iteration 0
err = 0.00902239888649081
Newton iteration 1
err = 3.435219434994371e-05
Newton iteration 2
err = 5.2950251633383934e-09
cost at iteration 443 : 6.632722552464204
Newton iteration 0
err = 0.009003101164597987
Newton iteration 1
err = 3.340001759933598e-05
Newton iteration 2
err = 4.931014188476379e-09
cost at iteration 444 : 6.62196045842637
Newton iteration 0
err = 0.008984246065244207
Newton iteration 1
err = 3.2479526286996525e-05
Newton iteration 2
err = 4.590070210266184e-09
cost at iteration 445 : 6.611217588502981
Newton iteration 0
err = 0.008965815898659309
Newton iteration 1
err = 3.159169560556029e-05
Newton iteration 2
err = 4.27176425636873e-09
cost at iteration 446 : 6.600493871393054
Newton iteration 0
err = 0.008947793200256985
Newton iteration 1
err = 3.073697562830319e-05
Newton iteration 2
err = 3.975364624295658e-09
cost at iteration 447 : 6.5897892375914
Newton iteration 0
err = 0.008930160880689299
Newton iteration 1
err = 2.9915401815835417e-05
Newton iteration 2
err = 3.6999362129556096e-09
cost at iteration 448 : 6.579103619309782
Newton iteration 0
err = 0.008912902334462796
Newton iteration 1
err = 2.9126686362901595e-05
Newton iteration 2
err = 3.4444105194300977e-09
cost at iteration 449 : 6.568436950401283
Newton iteration 0
err = 0.008896001515751624
Newton iteration 1
err = 2.837029322317834e-05
Newton iteration 2
err = 3.207646777244422e-09
cost at iteration 450 : 6.5577891662876535
Newton iteration 0
err = 0.008879442988206554
Newton iteration 1
err = 2.764549934571673e-05
Newton iteration 2
err = 2.9884727518469274e-09
cost at iteration 451 : 6.547160203889558
Newton iteration 0
err = 0.008863211954658834
Newton iteration 1
err = 2.695144436831966e-05
Newton iteration 2
err = 2.785716011732653e-09
cost at iteration 452 : 6.536550001560378
Newton iteration 0
err = 0.008847294271154552
Newton iteration 1
err = 2.6287170691975027e-05
Newton iteration 2
err = 2.5982305699336286e-09
cost at iteration 453 : 6.525958499023008
Newton iteration 0
err = 0.008831676449089276
Newton iteration 1
err = 2.565165560392926e-05
Newton iteration 2
err = 2.4249103937957684e-09
cost at iteration 454 : 6.515385637309393
Newton iteration 0
err = 0.008816345648428977
Newton iteration 1
err = 2.5043836842105253e-05
Newton iteration 2
err = 2.264696795784614e-09
cost at iteration 455 : 6.504831358704022
Newton iteration 0
err = 0.00880128966424696
Newton iteration 1
err = 2.446263277779022e-05
Newton iteration 2
err = 2.116591407039426e-09
cost at iteration 456 : 6.494295606689477
Newton iteration 0
err = 0.008786496908508004
Newton iteration 1
err = 2.3906958207483914e-05
Newton iteration 2
err = 1.979656308709132e-09
cost at iteration 457 : 6.483778325895214
Newton iteration 0
err = 0.008771956388557135
Newton iteration 1
err = 2.3375736568348085e-05
Newton iteration 2
err = 1.8530134861426442e-09
cost at iteration 458 : 6.473279462048986
Newton iteration 0
err = 0.008757657683433683
Newton iteration 1
err = 2.286790925154169e-05
Newton iteration 2
err = 1.7358485554588326e-09
cost at iteration 459 : 6.462798961931089
Newton iteration 0
err = 0.008743590918807293
Newton iteration 1
err = 2.238244257825825e-05
Newton iteration 2
err = 1.6274048275969789e-09
cost at iteration 460 : 6.4523367733305275
Newton iteration 0
err = 0.008729746741415662
Newton iteration 1
err = 2.1918332888598976e-05
Newton iteration 2
err = 1.5269863735847087e-09
cost at iteration 461 : 6.4418928450042365
Newton iteration 0
err = 0.008716116293235576
Newton iteration 1
err = 2.147461012897602e-05
Newton iteration 2
err = 1.43394844118948e-09
cost at iteration 462 : 6.431467126637739
Newton iteration 0
err = 0.00870269118607923
Newton iteration 1
err = 2.1050340236654653e-05
Newton iteration 2
err = 1.3476994690272858e-09
cost at iteration 463 : 6.421059568808557
Newton iteration 0
err = 0.008689463476587837
Newton iteration 1
err = 2.0644626577187454e-05
Newton iteration 2
err = 1.2676956135162848e-09
cost at iteration 464 : 6.410670122951312
Newton iteration 0
err = 0.0086764256420704
Newton iteration 1
err = 2.025661063570987e-05
Newton iteration 2
err = 1.1934382258702844e-09
cost at iteration 465 : 6.4002987413243035
Newton iteration 0
err = 0.00866357055724329
Newton iteration 1
err = 1.988547213042698e-05
Newton iteration 2
err = 1.124469396643642e-09
cost at iteration 466 : 6.389945376978945
Newton iteration 0
err = 0.008650891471888219
Newton iteration 1
err = 1.953042867864067e-05
Newton iteration 2
err = 1.0603704807272663e-09
cost at iteration 467 : 6.379609983729496
Newton iteration 0
err = 0.00863838198954823
Newton iteration 1
err = 1.9190735124584782e-05
Newton iteration 2
err = 1.0007571372939796e-09
cost at iteration 468 : 6.369292516125287
Newton iteration 0
err = 0.008626036047407462
Newton iteration 1
err = 1.8865682619561323e-05
Newton iteration 2
err = 9.452783660656187e-10
cost at iteration 469 : 6.35899292942415
Newton iteration 0
err = 0.008613847896995182
Newton iteration 1
err = 1.855459751285212e-05
Newton iteration 2
err = 8.936125954387467e-10
cost at iteration 470 : 6.348711179567033
Newton iteration 0
err = 0.008601812086249364
Newton iteration 1
err = 1.8256840121189622e-05
Newton iteration 2
err = 8.454640896520835e-10
cost at iteration 471 : 6.338447223154226
Newton iteration 0
err = 0.008589923442326298
Newton iteration 1
err = 1.7971803416188874e-05
Newton iteration 2
err = 8.005660369466255e-10
cost at iteration 472 : 6.328201017422577
Newton iteration 0
err = 0.008578177055572738
Newton iteration 1
err = 1.7698911661347036e-05
Newton iteration 2
err = 7.586692259153827e-10
cost at iteration 473 : 6.317972520223959
Newton iteration 0
err = 0.008566568264410546
Newton iteration 1
err = 1.7437619027418947e-05
Newton iteration 2
err = 7.195488910173281e-10
cost at iteration 474 : 6.307761690004833
Newton iteration 0
err = 0.008555092641115317
Newton iteration 1
err = 1.7187408210648528e-05
Newton iteration 2
err = 6.829977444390831e-10
cost at iteration 475 : 6.297568485786663
Newton iteration 0
err = 0.008543745978498404
Newton iteration 1
err = 1.6947789064823022e-05
Newton iteration 2
err = 6.488267847339164e-10
cost at iteration 476 : 6.287392867147724
Newton iteration 0
err = 0.00853252427738159
Newton iteration 1
err = 1.6718297263730075e-05
Newton iteration 2
err = 6.168620957784158e-10
cost at iteration 477 : 6.277234794205239
Newton iteration 0
err = 0.008521423734867509
Newton iteration 1
err = 1.6498492997510444e-05
Newton iteration 2
err = 5.869444113357366e-10
cost at iteration 478 : 6.2670942275986805
Newton iteration 0
err = 0.008510440733386096
Newton iteration 1
err = 1.6287959716786377e-05
Newton iteration 2
err = 5.58928071419119e-10
cost at iteration 479 : 6.256971128474109
Newton iteration 0
err = 0.008499571830323493
Newton iteration 1
err = 1.6086302922318067e-05
Newton iteration 2
err = 5.326785626901691e-10
cost at iteration 480 : 6.246865458468624
Newton iteration 0
err = 0.008488813748430869
Newton iteration 1
err = 1.589314900556803e-05
Newton iteration 2
err = 5.080745088045093e-10
cost at iteration 481 : 6.236777179696367
Newton iteration 0
err = 0.008478163366732912
Newton iteration 1
err = 1.5708144142116658e-05
Newton iteration 2
err = 4.850015411856928e-10
cost at iteration 482 : 6.226706254734282
Newton iteration 0
err = 0.00846761771211579
Newton iteration 1
err = 1.5530953237275977e-05
Newton iteration 2
err = 4.63356813103022e-10
cost at iteration 483 : 6.2166526466096155
Newton iteration 0
err = 0.00845717395130059
Newton iteration 1
err = 1.5361258927872325e-05
Newton iteration 2
err = 4.4304467900862085e-10
cost at iteration 484 : 6.206616318786624
Newton iteration 0
err = 0.008446829383521783
Newton iteration 1
err = 1.51987606308227e-05
Newton iteration 2
err = 4.239788635456096e-10
cost at iteration 485 : 6.196597235155458
Newton iteration 0
err = 0.008436581433458499
Newton iteration 1
err = 1.5043173642220036e-05
Newton iteration 2
err = 4.0607817567795717e-10
cost at iteration 486 : 6.186595360020166
Newton iteration 0
err = 0.008426427644746807
Newton iteration 1
err = 1.4894228295865868e-05
Newton iteration 2
err = 3.8926952467193375e-10
cost at iteration 487 : 6.176610658088096
Newton iteration 0
err = 0.00841636567388891
Newton iteration 1
err = 1.4751669153511225e-05
Newton iteration 2
err = 3.734847279917522e-10
cost at iteration 488 : 6.166643094459568
Newton iteration 0
err = 0.008406393284434764
Newton iteration 1
err = 1.4615254252981814e-05
Newton iteration 2
err = 3.5866158274829005e-10
cost at iteration 489 : 6.156692634617778
Newton iteration 0
err = 0.008396508341696638
Newton iteration 1
err = 1.448475439755905e-05
Newton iteration 2
err = 3.447423737542744e-10
cost at iteration 490 : 6.146759244419429
Newton iteration 0
err = 0.008386708807608539
Newton iteration 1
err = 1.4359952477180007e-05
Newton iteration 2
err = 3.316749170462532e-10
cost at iteration 491 : 6.136842890085618
Newton iteration 0
err = 0.008376992736070019
Newton iteration 1
err = 1.4240642847812438e-05
Newton iteration 2
err = 3.194115629701764e-10
cost at iteration 492 : 6.126943538193289
Newton iteration 0
err = 0.008367358268459696
Newton iteration 1
err = 1.4126630729310577e-05
Newton iteration 2
err = 3.079061333481607e-10
cost at iteration 493 : 6.117061155666889
Newton iteration 0
err = 0.008357803629449125
Newton iteration 1
err = 1.4017731648684683e-05
Newton iteration 2
err = 2.971182664492504e-10
cost at iteration 494 : 6.107195709770312
Newton iteration 0
err = 0.008348327123164852
Newton iteration 1
err = 1.391377091951684e-05
Newton iteration 2
err = 2.870110165981683e-10
cost at iteration 495 : 6.097347168099662
Newton iteration 0
err = 0.008338927129416925
Newton iteration 1
err = 1.3814583144599115e-05
Newton iteration 2
err = 2.775494303438451e-10
cost at iteration 496 : 6.087515498575674
Newton iteration 0
err = 0.008329602100302473
Newton iteration 1
err = 1.3720011752742464e-05
Newton iteration 2
err = 2.687010170108033e-10
cost at iteration 497 : 6.077700669436927
Newton iteration 0
err = 0.008320350556976226
Newton iteration 1
err = 1.3629908567620305e-05
Newton iteration 2
err = 2.6043770035983574e-10
cost at iteration 498 : 6.06790264923323
Newton iteration 0
err = 0.0083111710865116
Newton iteration 1
err = 1.3544133397933361e-05
Newton iteration 2
err = 2.527321151711883e-10
cost at iteration 499 : 6.058121406819125
Newton iteration 0
err = 0.00830206233922868
Newton iteration 1
err = 1.3462553650346539e-05
Newton iteration 2
err = 2.455595922548239e-10
cost at iteration 500 : 6.048356911348104
Newton iteration 0
err = 0.008293023025738679
Newton iteration 1
err = 1.3385043980063087e-05
Newton iteration 2
err = 2.3889797742188397e-10
cost at iteration 501 : 6.038609132266453
Newton iteration 0
err = 0.008284051914680703
Newton iteration 1
err = 1.3311485942530303e-05
Newton iteration 2
err = 2.327258024430791e-10
cost at iteration 502 : 6.028878039307869
Newton iteration 0
err = 0.00827514783016297
Newton iteration 1
err = 1.324176768565226e-05
Newton iteration 2
err = 2.2702396529257512e-10
cost at iteration 503 : 6.019163602488128
Newton iteration 0
err = 0.008266309649628545
Newton iteration 1
err = 1.31757836491911e-05
Newton iteration 2
err = 2.2177551338124054e-10
cost at iteration 504 : 6.009465792099981
Newton iteration 0
err = 0.008257536301712302
Newton iteration 1
err = 1.311343428739471e-05
Newton iteration 2
err = 2.1696341770016532e-10
cost at iteration 505 : 5.999784578708088
Newton iteration 0
err = 0.008248826764298986
Newton iteration 1
err = 1.3054625813108487e-05
Newton iteration 2
err = 2.1257229809326458e-10
cost at iteration 506 : 5.990119933144624
Newton iteration 0
err = 0.00824018006266088
Newton iteration 1
err = 1.2999269947508895e-05
Newton iteration 2
err = 2.0858964754675e-10
cost at iteration 507 : 5.980471826504749
Newton iteration 0
err = 0.008231595267720217
Newton iteration 1
err = 1.2947283698630873e-05
Newton iteration 2
err = 2.0500130161193352e-10
cost at iteration 508 : 5.970840230142147
Newton iteration 0
err = 0.00822307149442724
Newton iteration 1
err = 1.2898589142865036e-05
Newton iteration 2
err = 2.0179570040141952e-10
cost at iteration 509 : 5.961225115665147
Newton iteration 0
err = 0.008214607900230295
Newton iteration 1
err = 1.2853113229292838e-05
Newton iteration 2
err = 1.9896095202508026e-10
cost at iteration 510 : 5.9516264549327325
Newton iteration 0
err = 0.008206203683628625
Newton iteration 1
err = 1.2810787592554716e-05
Newton iteration 2
err = 1.9648758284656608e-10
cost at iteration 511 : 5.942044220050898
Newton iteration 0
err = 0.008197858082843155
Newton iteration 1
err = 1.2771548375813476e-05
Newton iteration 2
err = 1.9436502540704425e-10
cost at iteration 512 : 5.932478383368895
Newton iteration 0
err = 0.00818957037456236
Newton iteration 1
err = 1.2735336074268272e-05
Newton iteration 2
err = 1.9258555058362096e-10
cost at iteration 513 : 5.922928917476177
Newton iteration 0
err = 0.008181339872747042
Newton iteration 1
err = 1.2702095380704025e-05
Newton iteration 2
err = 1.9113932040433875e-10
cost at iteration 514 : 5.913395795198665
Newton iteration 0
err = 0.008173165927550229
Newton iteration 1
err = 1.2671775044791708e-05
Newton iteration 2
err = 1.900179976653705e-10
cost at iteration 515 : 5.903878989596287
Newton iteration 0
err = 0.008165047924262675
Newton iteration 1
err = 1.2644327745162479e-05
Newton iteration 2
err = 1.892138001167686e-10
cost at iteration 516 : 5.894378473959487
Newton iteration 0
err = 0.008156985282420606
Newton iteration 1
err = 1.2619709965349516e-05
Newton iteration 2
err = 1.8872181034163823e-10
cost at iteration 517 : 5.884894221806873
Newton iteration 0
err = 0.008148977454812233
Newton iteration 1
err = 1.2597881882418634e-05
Newton iteration 2
err = 1.8853274024777358e-10
cost at iteration 518 : 5.875426206882051
Newton iteration 0
err = 0.008141023926771343
Newton iteration 1
err = 1.2578807265819067e-05
Newton iteration 2
err = 1.8864247169331393e-10
cost at iteration 519 : 5.86597440315173
Newton iteration 0
err = 0.008133124215322533
Newton iteration 1
err = 1.2562453380909902e-05
Newton iteration 2
err = 1.890441721183539e-10
cost at iteration 520 : 5.856538784802926
Newton iteration 0
err = 0.00812527786853264
Newton iteration 1
err = 1.254879089948866e-05
Newton iteration 2
err = 1.8973303071631942e-10
cost at iteration 521 : 5.847119326240894
Newton iteration 0
err = 0.008117484464797473
Newton iteration 1
err = 1.253779382312144e-05
Newton iteration 2
err = 1.9070385625075613e-10
cost at iteration 522 : 5.83771600208676
Newton iteration 0
err = 0.008109743612354246
Newton iteration 1
err = 1.2529439407198066e-05
Newton iteration 2
err = 1.9195337192986783e-10
cost at iteration 523 : 5.828328787175868
Newton iteration 0
err = 0.008102054948624145
Newton iteration 1
err = 1.2523708093553935e-05
Newton iteration 2
err = 1.9347627094272422e-10
cost at iteration 524 : 5.818957656555715
Newton iteration 0
err = 0.008094418139791814
Newton iteration 1
err = 1.2520583456278075e-05
Newton iteration 2
err = 1.9527082172051233e-10
cost at iteration 525 : 5.809602585484208
Newton iteration 0
err = 0.008086832880349774
Newton iteration 1
err = 1.252005214323651e-05
Newton iteration 2
err = 1.973338657118097e-10
cost at iteration 526 : 5.800263549428046
Newton iteration 0
err = 0.008079298892687463
Newton iteration 1
err = 1.2522103827768975e-05
Newton iteration 2
err = 1.9966296784063297e-10
cost at iteration 527 : 5.790940524061151
Newton iteration 0
err = 0.008071815926739475
Newton iteration 1
err = 1.2526731168988648e-05
Newton iteration 2
err = 2.0225740056897487e-10
cost at iteration 528 : 5.781633485263315
Newton iteration 0
err = 0.008064383759728776
Newton iteration 1
err = 1.2533929774933918e-05
Newton iteration 2
err = 2.051162981861608e-10
cost at iteration 529 : 5.7723424091189
Newton iteration 0
err = 0.008057002195855688
Newton iteration 1
err = 1.2543698169531945e-05
Newton iteration 2
err = 2.0823880407221202e-10
cost at iteration 530 : 5.763067271915612
Newton iteration 0
err = 0.00804967106607208
Newton iteration 1
err = 1.2556037765337035e-05
Newton iteration 2
err = 2.1162763254148652e-10
cost at iteration 531 : 5.753808050143435
Newton iteration 0
err = 0.008042390227981545
Newton iteration 1
err = 1.2570952845141713e-05
Newton iteration 2
err = 2.1528023440611236e-10
cost at iteration 532 : 5.744564720493715
Newton iteration 0
err = 0.008035159565601333
Newton iteration 1
err = 1.2588450540062183e-05
Newton iteration 2
err = 2.19201662556367e-10
cost at iteration 533 : 5.735337259858111
Newton iteration 0
err = 0.008027978989361778
Newton iteration 1
err = 1.2608540819728492e-05
Newton iteration 2
err = 2.2339403972087974e-10
cost at iteration 534 : 5.726125645328157
Newton iteration 0
err = 0.008020848435981729
Newton iteration 1
err = 1.2631236481362172e-05
Newton iteration 2
err = 2.2785938479751314e-10
cost at iteration 535 : 5.716929854194277
Newton iteration 0
err = 0.00801376786849668
Newton iteration 1
err = 1.265655314609669e-05
Newton iteration 2
err = 2.3260337412464485e-10
cost at iteration 536 : 5.707749863945526
Newton iteration 0
err = 0.00800673727625705
Newton iteration 1
err = 1.2684509253905609e-05
Newton iteration 2
err = 2.3763050802503653e-10
cost at iteration 537 : 5.698585652268866
Newton iteration 0
err = 0.007999756675014454
Newton iteration 1
err = 1.2715126070717521e-05
Newton iteration 2
err = 2.4294517762429433e-10
cost at iteration 538 : 5.689437197049102
Newton iteration 0
err = 0.007992826106975102
Newton iteration 1
err = 1.2748427687152409e-05
Newton iteration 2
err = 2.485554630297188e-10
cost at iteration 539 : 5.680304476368519
Newton iteration 0
err = 0.007985945641007066
Newton iteration 1
err = 1.2784441030841315e-05
Newton iteration 2
err = 2.5446898429135946e-10
cost at iteration 540 : 5.671187468506765
Newton iteration 0
err = 0.007979115372743435
Newton iteration 1
err = 1.2823195874343547e-05
Newton iteration 2
err = 2.606926238559937e-10
cost at iteration 541 : 5.662086151940796
Newton iteration 0
err = 0.00797233542483054
Newton iteration 1
err = 1.2864724847016013e-05
Newton iteration 2
err = 2.672353276259297e-10
cost at iteration 542 : 5.65300050534515
Newton iteration 0
err = 0.007965605947178054
Newton iteration 1
err = 1.2909063450312273e-05
Newton iteration 2
err = 2.74107773722818e-10
cost at iteration 543 : 5.643930507591898
Newton iteration 0
err = 0.00795892711719301
Newton iteration 1
err = 1.2956250069902644e-05
Newton iteration 2
err = 2.813205765179934e-10
cost at iteration 544 : 5.634876137751178
Newton iteration 0
err = 0.007952299140198101
Newton iteration 1
err = 1.3006325998573172e-05
Newton iteration 2
err = 2.888849659357409e-10
cost at iteration 545 : 5.625837375091374
Newton iteration 0
err = 0.007945722249692583
Newton iteration 1
err = 1.3059335445839894e-05
Newton iteration 2
err = 2.9681386544520436e-10
cost at iteration 546 : 5.616814199079996
Newton iteration 0
err = 0.007939196707749653
Newton iteration 1
err = 1.3115325558048631e-05
Newton iteration 2
err = 3.0512172066967435e-10
cost at iteration 547 : 5.607806589383843
Newton iteration 0
err = 0.007932722805532439
Newton iteration 1
err = 1.3174346434960832e-05
Newton iteration 2
err = 3.1382175503314044e-10
cost at iteration 548 : 5.598814525870272
Newton iteration 0
err = 0.007926300863622006
Newton iteration 1
err = 1.3236451140286027e-05
Newton iteration 2
err = 3.229308888063699e-10
cost at iteration 549 : 5.589837988607559
Newton iteration 0
err = 0.00791993123262582
Newton iteration 1
err = 1.3301695717980074e-05
Newton iteration 2
err = 3.324651266717776e-10
cost at iteration 550 : 5.580876957866417
Newton iteration 0
err = 0.007913614293626985
Newton iteration 1
err = 1.3370139197894228e-05
Newton iteration 2
err = 3.4244398906814325e-10
cost at iteration 551 : 5.571931414120602
Newton iteration 0
err = 0.007907350458765837
Newton iteration 1
err = 1.3441843600033422e-05
Newton iteration 2
err = 3.5288486788560735e-10
cost at iteration 552 : 5.563001338048367
Newton iteration 0
err = 0.007901140171855614
Newton iteration 1
err = 1.3516873936999256e-05
Newton iteration 2
err = 3.638091528859508e-10
cost at iteration 553 : 5.554086710534076
Newton iteration 0
err = 0.007894983908919864
Newton iteration 1
err = 1.3595298198831245e-05
Newton iteration 2
err = 3.7523779567046767e-10
cost at iteration 554 : 5.545187512669055
Newton iteration 0
err = 0.007888882178938204
Newton iteration 1
err = 1.367718734430061e-05
Newton iteration 2
err = 3.8719343013523494e-10
cost at iteration 555 : 5.536303725753673
Newton iteration 0
err = 0.007882835524427284
Newton iteration 1
err = 1.3762615266295192e-05
Newton iteration 2
err = 3.9970136935238735e-10
cost at iteration 556 : 5.527435331298962
Newton iteration 0
err = 0.007876844522166226
Newton iteration 1
err = 1.3851658759751167e-05
Newton iteration 2
err = 4.1278758662906064e-10
cost at iteration 557 : 5.51858231102845
Newton iteration 0
err = 0.007870909783886219
Newton iteration 1
err = 1.3944397461364093e-05
Newton iteration 2
err = 4.2647712264472146e-10
cost at iteration 558 : 5.509744646879991
Newton iteration 0
err = 0.007865031957033157
Newton iteration 1
err = 1.4040913784978192e-05
Newton iteration 2
err = 4.407999286588591e-10
cost at iteration 559 : 5.500922321008218
Newton iteration 0
err = 0.007859211725439227
Newton iteration 1
err = 1.414129282629913e-05
Newton iteration 2
err = 4.5578684482680946e-10
cost at iteration 560 : 5.492115315786662
Newton iteration 0
err = 0.007853449810081781
Newton iteration 1
err = 1.4245622247571616e-05
Newton iteration 2
err = 4.714686960569678e-10
cost at iteration 561 : 5.483323613810045
Newton iteration 0
err = 0.007847746969824118
Newton iteration 1
err = 1.4353992144917534e-05
Newton iteration 2
err = 4.878792651897743e-10
cost at iteration 562 : 5.474547197896926
Newton iteration 0
err = 0.007842104002147191
Newton iteration 1
err = 1.4466494864044308e-05
Newton iteration 2
err = 5.05053755354408e-10
cost at iteration 563 : 5.465786051092438
Newton iteration 0
err = 0.007836521743773838
Newton iteration 1
err = 1.4583224796828617e-05
Newton iteration 2
err = 5.230283390883619e-10
cost at iteration 564 : 5.457040156671103
Newton iteration 0
err = 0.007831001071461322
Newton iteration 1
err = 1.4704278134402176e-05
Newton iteration 2
err = 5.418417227525245e-10
cost at iteration 565 : 5.448309498139927
Newton iteration 0
err = 0.007825542902524175
Newton iteration 1
err = 1.4829752559794658e-05
Newton iteration 2
err = 5.615346262168241e-10
cost at iteration 566 : 5.439594059241327
Newton iteration 0
err = 0.007820148195435493
Newton iteration 1
err = 1.4959746901534555e-05
Newton iteration 2
err = 5.821485959323693e-10
cost at iteration 567 : 5.430893823956671
Newton iteration 0
err = 0.007814817950382491
Newton iteration 1
err = 1.5094360711366724e-05
Newton iteration 2
err = 6.037274939976459e-10
cost at iteration 568 : 5.422208776509639
Newton iteration 0
err = 0.007809553209530245
Newton iteration 1
err = 1.5233693778120226e-05
Newton iteration 2
err = 6.263164951838951e-10
cost at iteration 569 : 5.413538901369829
Newton iteration 0
err = 0.0078043550574769666
Newton iteration 1
err = 1.537784555891436e-05
Newton iteration 2
err = 6.499626483195095e-10
cost at iteration 570 : 5.404884183256695
Newton iteration 0
err = 0.007799224621228297
Newton iteration 1
err = 1.5526914508877222e-05
Newton iteration 2
err = 6.74715171860191e-10
cost at iteration 571 : 5.3962446071432195
Newton iteration 0
err = 0.00779416307028947
Newton iteration 1
err = 1.5680997316618e-05
Newton iteration 2
err = 7.006223676273083e-10
cost at iteration 572 : 5.387620158260255
Newton iteration 0
err = 0.007789171616364521
Newton iteration 1
err = 1.584018800747182e-05
Newton iteration 2
err = 7.277369282406848e-10
cost at iteration 573 : 5.379010822100767
Newton iteration 0
err = 0.007784251512838803
Newton iteration 1
err = 1.600457691846071e-05
Newton iteration 2
err = 7.561103899674247e-10
cost at iteration 574 : 5.370416584423978
Newton iteration 0
err = 0.007779404054052197
Newton iteration 1
err = 1.6174249511360867e-05
Newton iteration 2
err = 7.857960469032893e-10
cost at iteration 575 : 5.361837431260335
Newton iteration 0
err = 0.007774630574098501
Newton iteration 1
err = 1.6349285014340907e-05
Newton iteration 2
err = 8.168463552946523e-10
cost at iteration 576 : 5.353273348915815
Newton iteration 0
err = 0.007769932445300343
Newton iteration 1
err = 1.6529754863235735e-05
Newton iteration 2
err = 8.493155812188167e-10
cost at iteration 577 : 5.344724323977165
Newton iteration 0
err = 0.007765311076140512
Newton iteration 1
err = 1.671572092019712e-05
Newton iteration 2
err = 8.832562938235196e-10
cost at iteration 578 : 5.336190343316463
Newton iteration 0
err = 0.007760767908698464
Newton iteration 1
err = 1.690723343824332e-05
Newton iteration 2
err = 9.187191244633525e-10
cost at iteration 579 : 5.327671394096573
Newton iteration 0
err = 0.007756304415314915
Newton iteration 1
err = 1.710432874518552e-05
Newton iteration 2
err = 9.557545206308228e-10
cost at iteration 580 : 5.3191674637760435
Newton iteration 0
err = 0.007751922094658117
Newton iteration 1
err = 1.7307026609910238e-05
Newton iteration 2
err = 9.94408760495309e-10
cost at iteration 581 : 5.310678540114642
Newton iteration 0
err = 0.007747622466691406
Newton iteration 1
err = 1.751532725049221e-05
Newton iteration 2
err = 1.0347237832010814e-09
cost at iteration 582 : 5.302204611178359
Newton iteration 0
err = 0.007743407066873933
Newton iteration 1
err = 1.7729207949549903e-05
Newton iteration 2
err = 1.0767381359761941e-09
cost at iteration 583 : 5.293745665345269
Newton iteration 0
err = 0.007739277438887648
Newton iteration 1
err = 1.794861922410038e-05
Newton iteration 2
err = 1.1204816258772703e-09
cost at iteration 584 : 5.285301691310395
Newton iteration 0
err = 0.007735235126307767
Newton iteration 1
err = 1.817348051365307e-05
Newton iteration 2
err = 1.165976730542969e-09
cost at iteration 585 : 5.276872678091341
Newton iteration 0
err = 0.007731281662498428
Newton iteration 1
err = 1.840367531779366e-05
Newton iteration 2
err = 1.213234754305661e-09
cost at iteration 586 : 5.268458615033797
Newton iteration 0
err = 0.007727418558806001
Newton iteration 1
err = 1.8639045741658906e-05
Newton iteration 2
err = 1.2622550170965943e-09
cost at iteration 587 : 5.260059491816233
Newton iteration 0
err = 0.007723647290796462
Newton iteration 1
err = 1.88793863929605e-05
Newton iteration 2
err = 1.3130202873279214e-09
cost at iteration 588 : 5.25167529845525
Newton iteration 0
err = 0.007719969281987918
Newton iteration 1
err = 1.9124437562421617e-05
Newton iteration 2
err = 1.3654949555959318e-09
cost at iteration 589 : 5.243306025310007
Newton iteration 0
err = 0.007716385885252525
Newton iteration 1
err = 1.9373877643156313e-05
Newton iteration 2
err = 1.4196215305525536e-09
cost at iteration 590 : 5.234951663087291
Newton iteration 0
err = 0.007712898360883163
Newton iteration 1
err = 1.962731471352546e-05
Newton iteration 2
err = 1.475317226185749e-09
cost at iteration 591 : 5.226612202844605
Newton iteration 0
err = 0.00770950785168892
Newton iteration 1
err = 1.988427726035643e-05
Newton iteration 2
err = 1.5324675369995589e-09
cost at iteration 592 : 5.218287635994532
Newton iteration 0
err = 0.007706215353939754
Newton iteration 1
err = 2.014420396946258e-05
Newton iteration 2
err = 1.5909243586912085e-09
cost at iteration 593 : 5.209977954307041
Newton iteration 0
err = 0.007703021684454247
Newton iteration 1
err = 2.0406432573351692e-05
Newton iteration 2
err = 1.6504995100852887e-09
cost at iteration 594 : 5.201683149911779
Newton iteration 0
err = 0.007699927442642521
Newton iteration 1
err = 2.0670187738890633e-05
Newton iteration 2
err = 1.710957520083692e-09
cost at iteration 595 : 5.193403215299303
Newton iteration 0
err = 0.0076969329675687145
Newton iteration 1
err = 2.09345680062105e-05
Newton iteration 2
err = 1.772012365723289e-09
cost at iteration 596 : 5.1851381433209305
Newton iteration 0
err = 0.007694038289053742
Newton iteration 1
err = 2.1198531828448528e-05
Newton iteration 2
err = 1.8333207316496864e-09
cost at iteration 597 : 5.176887927187565
Newton iteration 0
err = 0.007691243072542982
Newton iteration 1
err = 2.1460882818521602e-05
Newton iteration 2
err = 1.8944752194511004e-09
cost at iteration 598 : 5.168652560466924
Newton iteration 0
err = 0.007688546557093541
Newton iteration 1
err = 2.1720254346073952e-05
Newton iteration 2
err = 1.9549982557593685e-09
cost at iteration 599 : 5.160432037078894
Newton iteration 0
err = 0.007685947485827505
Newton iteration 1
err = 2.1975093732839245e-05
Newton iteration 2
err = 2.014337234123511e-09
=================Mesh Quality Check=================
max quality: 1.37841365416341 min quality: 0.5553651552984573 mean quality: 0.858318409251656
===================================================